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Abstract 

A unified presentation is given on the use of dispersion relations in the real and virtual 
Compton scattering processes off the nucleon. The way in which dispersion relations for 
Compton scattering amplitudes establish connections between low energy nucleon struc- 
ture quantities, such as polarizabilities or anomalous magnetic moments, and the nucleon 
excitation spectrum is reviewed. We discuss various sum rules for forward real and virtual 
Compton scattering, such as the Gerasimov-Drell-Hearn sum rule and its generalizations, 
the Burkhardt-Cottingham sum rule, as well as sum rules for forward nucleon polarizabili- 
ties, and review their experimental status. Subsequently, we address the general case of real 
Compton scattering (RCS). Various types of dispersion relations for RCS are presented as 
tools for extracting nucleon polarizabilities from the RCS data. The information on nu- 
cleon polarizabilities gained in this way is reviewed and the nucleon structure information 
encoded in these quantities is discussed. The dispersion relation formalism is then extended 
to virtual Compton scattering (VCS). The information on generalized nucleon polarizabil- 
ities extracted from recent VCS experiments is described, along with its interpretation in 
nucleon structure models. As a summary, the physics content of the existing data is dis- 
cussed and some perspectives for future theoretical and experimental activities in this field 
are presented. 
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1 Introduction 



The internal structure of the strongly interacting particles has been an increasingly active 
area of experimental and theoretical research over the past 5 decades. Precision experiments 
at high energy have clearly established Quantum Chromodynamics (QCD) as the underlying 
gauge theory describing the interaction between quarks and gluons, the elementary constituents 
of hadronic matter. However, the running coupling constant of QCD grows at low energies, 
and these constituents are confined to colorless hadrons, the mesons and baryons, which are 
the particles eventually observed by the detection devices. Therefore, we have to live with a 
dichotomy: The small value of the coupling constant at high energies allows for an interpre- 
tation of the experiments in terms of perturbative QCD, while the large value at low energies 
calls for a description in terms of the hadronic degrees of freedom, in particular in the approach 
developed as Chiral Perturbation Theory. 

Between these two regions, at excitation energies between a few hundred MeV and 1-2 GeV, 
lies the interesting region of nucleon resonance structures which is beyond the scope of ei- 
ther perturbation scheme. There is some hope that this regime will eventually be described by 
numerical solutions of QCD through lattice gauge calculations. At present, however, our un- 
derstanding of resonance physics is still mostly based on phenomenology. In the absence of 
a descriptive theory it is essential to extract new and precise hadronic structure information, 
and in this quest electromagnetic probes have played a decisive role. In particular, high preci- 
sion Compton scattering experiments have become possible with the advent of modern electron 
accelerators with high current and duty factor, and of laser backscattering facilities, and in com- 
bination with high precision and large acceptance detectors. This intriguing new window offers, 
among other options, the possibility for precise and detailed investigations of the nucleon po- 
larizability as induced by the applied electromagnetic multipole fields. 

The polarizability of a composite system is an elementary structure constant, just as are its 
size and shape. In a macroscopic medium, the electric and magnetic dipole polarizabilities are 
related to the dielectric constant and the magnetic permeability, and these in turn determine the 
index of refraction. These quantities can be studied by considering an incident electromagnetic 
wave inducing dipole oscillations in the constituent atoms or molecules of a target medium. 
These oscillations then emit dipole radiation leading, by way of interference with the incoming 
wave, to the complex amplitude of the transmitted wave. A general feature of these processes is 
the dispersion relation of Kronig and Kramers [1], which connects the real refraction index as 
function of the frequency with a weighted integral of the extinction coefficient over all frequen- 
cies. 

Dispersion theory in general relies on a few basic principles of physics: relativistic covari- 
ance, causality and unitarity. As a first step a complete set of amplitudes has to be constructed, in 
accordance with relativity and without kinematical singularities. Next, causality requires certain 
analytic properties of the amplitudes, which allow for a continuation of the scattering amplitudes 
into the complex plane and lead to dispersion relations connecting the real and imaginary parts 
of these amplitudes. Finally, the imaginary parts can be replaced by absorption cross sections 
by the use of unitarity, and as a result we can, for example, complete the Compton amplitudes 
from experimental information on photoabsorption and photo-induced reactions. 

In the following Sec. 2 we first discuss the classical theory of dispersion and absorption in 
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a medium, and briefly compare the polarizability of macroscopic matter and microscopic sys- 
tems, atoms and nucleons. This is followed by a review of forward Compton scattering and its 
connection to total absorption cross sections. Combining dispersion relations and low energy 
theorems, we obtain sum rules for certain combinations of the polarizabilities and other ground 
state properties, e.g., the Gerasimov-Drell-Heam sum rule for real photons [2,3], and the much 
debated Burkhardt-Cottingham sum rule for virtual Compton scattering [4] as obtained from 
radiative electron scattering. 

We then address the general case of real Compton scattering in Sec. 3. Besides the elec- 
tric and magnetic (dipole) polarizabilities of a scalar system, the spin of the nucleon leads to 
four additional spin or vector polarizabilities, and higher multipole polarizabilities will appear 
with increasing photon energy. We show how these polarizabilities can be obtained from pho- 
ton scattering and photoexcitation processes through a combined analysis based on dispersion 
theory. The results of such an analysis are then compared in detail with the experimental data 
and predictions from theory. In the Sec. 4 we discuss the more general case of virtual Comp- 
ton scattering, which can be achieved by radiative electron-proton scattering. Such experiments 
have become possible only very recently. The non-zero four-momentum transfer squared of the 
virtual photon allows us to study generalized polarizabilities as function of four-momentum 
transfer squared and therefore, in some sense, to explore the spatial distribution of the polariza- 
tion effects. In the last Section, we summarize the pertinent features of our present knowledge 
on the nucleon polarizability and conclude by outlining some remaining challenges for future 
work. 

This review is largely based on dispersion theory whose development is related to Heisen- 
berg's idea that the interaction of particles can be described by their behavior at large distances, 
i.e., in terms of the S matrix [5]. The practical consequences of this program were worked out 
by Mandelstam and others [6]. An excellent primer for the beginner is the textbook of Nussen- 
zveig [7] . In order to feel comfortable on Mandelstam planes and higher Riemann sheets, the 
review of Hoehler [8] is an absolute must for the practitioner. Concerning the structure aspect of 
our review, we refer the reader to a general treatise of the electromagnetic response of hadronic 
systems by Boffi et al. [9], and to the recent book of Thomas and Weise [10], which is focused 
on the structure aspects of the nucleon. 
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2 Forward dispersion relations and sum rules for real and virtual Compton scattering 



2.1 Classical theory of dispersion and absorption in a medium 



The classical theory of Lorentz describes the dispersion in a medium in terms of electrons 
bound by a harmonic force. In the presence of a monochromatic external field, E^, the equations 
of motion take the form 



+ 27,4 + ^i] = -- e-^-* , (1) 



m 



with — e the chargeQ and m the mass of the electron, and 7^ > and cuj > the damping con- 
stant and oscillator frequency, respectively, of a specific bound state j. The stationary solution 
for the displacement is then given by 

e E^^ e"**^* 

^j(^) = — ^ ' (2) 

and the polarization P is obtained by summing the individual dipole moments dj = — e rj over 
all electrons and oscillator frequencies in the medium, 

m = E , I ; = p. e-^-* , (3) 

J m{ujj - 2i'jj to - co^) 

where A'^ is the number of electrons per unit volume, in the state j. The dielectric susceptibility 
X is defined by 

P^ = x(^)E. , (4) 

with 



LOj — 2i 7j 00 — LO^ 



We observe at this point that xi^) 
(I) is square integrable in the upper half-plane (/+) for any line parallel to the real to axis, and 



^ In Sec. 2.1 we shall use Gaussian units as in most of the literature on theoretical electrodynamics, i.e., 
the fine structure constant takes the form Oem = e'^/ch « 1/137 and the classical electron radius is 
rd = /mc^. In all later sections the Heaviside-Lorentz units will be used in order to concur with the 
standard notation of particle physics. 
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(II) has singularities only in the lower-half plane (/-) in the form of pairs of poles at 



w± = ±V^j - 7j - ilj ■ (6) 

According to Titchmarsh's theorem these observations have the following consequences: 
The Fourier transform 



oo 

Xit) = - I Xiuj) e~'^' du; (7) 



— oo 



is causal, i.e., the dielectric susceptibility and the polarization of the medium build up only after 
the electric field is applied, and the real and imaginary parts of x are Hilbert transforms, 

TT J UJ' — UJ 

— oo 

oo / l\ 

= (8) 

TT J Uj' — U! 
— oo 

where V denotes the principal value integral. 
Applying the convolution theorem for Fourier transforms to Eq. (4), we obtain 



PW= j xit-t')E{t')dt' , (9) 

— oo 

with general time profiles P(t) and E(t) of medium polarization and external field, respectively, 
constructed according to Eq. (7). 

The proof of causality follows from integrating the dielectric susceptibility over a contour C+ 
along the real uj axis, for —R < lo < R, and closed by a large half circle with radius R in the 
upper part of the complex a;-plane. Since no singularities appear within this contour. 



c+ 



'dco^O. (10) 



We make contact with the Fourier transform of Eq. (7) by blowing up the contour (i? — > oo) and 
studying the convergence along the half circle. According to our observation (I) the function x 
itself is square integrable in /+, and therefore the convergence depends on the behavior of the 
exponential function exp(— icjr), which depends on the sign of r. In the case of r < the 
convergence is improved by the exponential, and the contribution of the half-circle vanishes in 
the limit i? — > oo. Combining Eqs. (7) and (10), we then obtain 

X(r) = forrO, (11) 
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which enforces causality, as becomes obvious by inspecting Eq. (9): The electric field E(f') will 
affect the polarizability P(t) only at some later time, r = t — t' > 0. For such time, r > 0, 
the contour integral C+ is of course useless for our purpose, because the exponential overrides 
the convergence of % in /_|_. Therefore, the contour has to be closed in the lower half -plane, 
which picks up the contributions from the singularities in 7_. We note in passing that Eq. (11) 
describes the nonrelativistic causality condition, which has to be sharpened by the postulate 
of relativity that no signal can move faster than the velocity of light. Furthermore, causality is 
found to be a direct consequence of analyticity of the Green function x(a;), which in the Lorentz 
model results from the choice of 7^. For 7^ < 0, the poles of % would have moved to the upper 
half -plane of u, and the result would be an acausal response, x(t) > for t < and xi^) — 
for T > 0. 

Next let us study the symmetry properties of x under the ("crossing") transformation uj — > 
—UJ. The real (xr) imaginary (xi) parts of this function can be read off Eq. (5), 



--VAT, 



2^2 + 47^6^2 



(12) 



xM = -Y.N,r-2 — 1¥tt^' (13) 

m J (a;2 — ujj)^ + 47^ uj^ 
and the crossing relations for real uj values are 

This makes it possible to cast Eq. (8) into the form 

J U}'^ — U}^ IT J 



The crossing relations Eq. (14) can be combined and extended to complex values of uj by 



TT J UJ'^ — UJ^ IT J UJ'^ — UJ^ 



X{-UJ*) = . (16) 

In particular, x is real on the imaginary axis and takes on complex conjugate values at points 
situated mirror-symmetrically to this axis. The dielectric susceptibility can be expressed by the 
dielectric constant e, 

- ^ . (17) 
which in turn is related to the refraction index n and the phase velocity vp'm the medium. 
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where k is the wave number, and the magnetic permeability of the medium. 
In the case of // = 1, it is obvious that also {e — 1) and hence (n^ — 1) obey the dispersion 
relations of Eq. (15). In a gas of low density, the refraction index is close to 1, and we can 
approximate (n^ — 1) by 2(n 1). The result is the Lorentz dispersion formula for the oscillator 
model, to be obtained from Eqs. (5), (17) and (18), 



n 



(u;) = l + 27r-E 



m 



(19) 



Let us now discuss the connection between absorption and dispersion on the microscopic level. 
Suppose that a monochromatic plane wave hits a homogeneous and isotropic medium atx — 
and leaves the slab of matter atx — Ax. The incoming wave is denoted by 



Ei„(x,t) = e^('=^-'^*)£;oeo, (20) 

with the linear dispersion uj = ck and the polarization vector gq. 

Having passed the slab of matter with the dispersion of Eq. (19), the wave function is 



^gi^(n«-l)Aa;g-^n,Aa;E.^^Aa;,i) . (21) 

The imaginary part of n is associated with absorption, which defines an extinction coefficient 
K, such that the intensity drops like |Eoutp = e~'*^^|Einp. On the other hand the extinction 
coefficient is related to the product of the total absorption cross section ctt for an individual 
constituent (e.g., a atom) and the number of constituents per volume N, and therefore 



k{u!) — 2ujni/c — N(7t{oj) . (22) 

Further on the elementary level, the incident light wave excites dipole oscillations of the con- 
stituents with electric dipole moments 



d(t) =aEi„(0,t) , (23) 

with oi — a{uS) the electric dipole polarizability of a constituent. We note that here and in the 

following the dipole approximation has been used such that we can neglect retardation effects 
and evaluate the incoming wave at x = 0. Within the slab of matter, the dipole moments radiate, 
thus giving rise to an induced electric field E^. 
The field due to the individual dipole at r', measured at a point r = x e-^ in beam direction, is 
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gi(fcp-a)t) 

= ak^ Eo {q x eo) x g , (24) 

P 

with ^ = (r' — r)/|r' — r| and p = |r' — r|. 

In particular, forward scattering is obtained in the limit /cx ^ 1. Since the incoming field is 
polarized perpendicularly to this axis, we find 

e, = 0) = £;oeo , (25) 

P 

and by definition the forward scattering amplitude 



f{k,e = 0) = ak"^ . (26) 

The total field due to the dipole oscillations, Eg, is obtained by integrating Eq. (24) over the 
volume of the slab and multiplying with A^, the number of particles per volume. The result for 
small Ax is 



Eont^'Ein + Es^{l + 2mkAxNa)Ei^. (27) 

A comparison of Eqs. (26) and (27) with the macroscopic form, Eq. (21), expanded for small 
Ax, yields the connection between the refractive index and the forward scattering amplitude. 



n{uj)-l = 2TiNa{uj) = ^^f{k,e = 0). (28) 
From Eqs. (22) and (28) we obtain the optical theorem. 



Im/H = ^aTH, (29) 

and since ///c^ is proportional to (n — 1) and there follows a dispersion relation for Re/ 
analogous to Eq. (15), 



Re/M = r ,t 2 2^ d^ = / J ' (30) 

^ ^ 

where we have set c = /i = 1 here and in the following. Historically, Eq. (30) expressed in terms 

of n{{jS) — 1, was first derived by Kronig and Kramers [1]. We also note that without the crossing 
symmetry, Eq. (14), the dispersion integral would also need information about the cross section 
at negative energies, which of course is not available. 
In order to prepare for the specific content of this review, several comments are in order: 
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(I) The derivation of the Kramers-Kronig dispersion relation started from a neutral system, an 
atom like the hydrogen atom. Since the total charge is zero, the electromagnetic field can 
only excite the internal degrees of freedom, while the center of mass remains fixed. As a 
consequence the scattering amplitude /(a;) = 0{uj'^), which leads to a differential cross 
section 

The result is Rayleigh scattering which among other things explains the blue sky. How- 
ever, for charged systems like ions, electrons or protons, also the center of mass will be 
accelerated by the electromagnetic field, and the scattering amplitude takes the general 
form 



Q 



Ref{uj,Q) = --^ + 0{uj') . (31) 



tot 



The additional "Thomson" term due to cm. motion results in a finite scattering amplitude 

for u; = and depends only on the total charge Qtot and the total mass Mtot- 
(II) We have defined the electric dipole polarizability a as a complex function a{uj) whose 
real and imaginary parts can be calculated directly from the total absorption cross section 
ariuj). In the Lorentz model this cross section starts as u;^ for small lo. In reality, however, 
the total absorption cross section has a threshold energy loq. The absorption spectrum of, 
say, hydrogen is given by a series of discrete levels (cuis^ip = 10.2 eV, etc.) followed by 
a continuum for u > 13.6 eV. As a result a^iuj) vanishes in a range < cl^ < uq, and 
therefore a (a; ) = Re a (a;) can be expanded in a Taylor series in the vicinity of the origin, 

"(^) = ^ y + 2^ y (32) 

In the following chapters we shall use the term "polarizability" or more exactly "static 
polarizability" only for the first term of the expansion. Moreover, in the dipole expansion 
used in Eq. (23), this first term is solely determined by electric dipole (El) radiation. 



a 



1 f ariu! 



= 0) = — / -^^div' > . (33) 



0)0 

The terms 0{uj'^) in Eq. (32) are then the first order retardation effects for El radiation, 
and the full function a{uj) will be called the "dynamical polarizability" of the system. 
(Ill) Finally, the Lorentz model discards magnetic effects because of the small velocities in- 
volved in atomic systems. In a general derivation, the first term on the rhs of Eq. (32) 
equals the sum of the electric (a) and magnetic (/?) dipole polarizabilities, while the sec- 
ond term describes the retardation of these dipole polarizabilities and the static quadrupole 
polarizabilities. 

Let us finally discuss the polarizability for some specific cases. The Hamiltonian for an elec- 
tron bound by a harmonic restoring force, as in the Lorentz model of Eq. (1), takes the form 
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//= P^ + ^r^ + er-E, (34) 
2m 2 

where the electric field E is assumed to be static and uniform. Substituting r = r' + Ar and 
p = p', where Ar is the displacement due to the electric field, we may rewrite this equation as 



H^^ + ^r- + AE. (35) 
2m 2 

The displacement Ar leads to an induced dipole moment d and an energy shift AE, 



2 2 

d = -eAr = E , A^ = ^ E^ . (36) 

mwQ 2mwQ 



The induced dipole moment d and the energy shift AE are both proportional to the polarizabil- 
ity, a — e'^/mujQ, which can also be read off Eqs. (2) and (23) in the limit a; — > 0. In fact, the 
relation 

6d 5^AE 

is quite general and even survives in quantum mechanics. As a result we can calculate the 

energy of such a system by second order perturbation theory. The perturbation to first order 
(linear Stark effect) vanishes for a system with good parity, and if the system is also spherically 
symmetric, the second order (quadratic Stark effect) yields, 

AE^-Y, ^ ^ — ^ L E^ , (38) 



n>0 



eo 



where are the energies of the eigenstate \n>. Equations (37) and (38) immediately yield the 
static electric dipole polarizability. 



^^^^ |<nMO>P _ (39) 

As an example for a classical extended object we quote the electric (a) and magnetic {(3) dipole 
polarizabilities of small dielectric or permeable spheres of radius a [11], 

= (40) 

e + 2 + 2 

The same quantities for a perfectly conducting sphere are obtained in the limits e — > cxd and 
// — > 0, respectively. 
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Q; = a^ p^--a^. (41) 

The electric polarizability of the conducting sphere is essentially the volume of the sphere, 
up to a factor 47r/3. Due to the different boundary conditions, the magnetic polarizability is 
negative, which corresponds to diamagnetism (/U < 1). In this case the currents and with them 
the magnetizations are induced against the direction of the applied field according to Lenz's 
law. A permeable sphere can be diamagnetic or paramagnetic (// > 1), in the latter case the 
magnetic moments are already preformed and become aligned in the presence of the external 
field. While the magnetic polarizabilities of atoms and molecules are usually very small because 
of l/i — 1| < lO^'^, electric polarizabilities may be quite large compared to the volume. For 
example, the static dielectric constant of water e = 81 leads to a nearly perfect conductor; in the 
visible range this constant is down to £ = 1.8 with the consequence that the index of refraction 
isn = 1.34. 

A quantum mechanical example is the hydrogen atom in nonrelativistic description. Its ground 
state has good parity and spherical symmetry and therefore Eq. (38 ) applies. In this case it is 
even possible to perform the sum over the excited states and to obtain the closed expression [12] 



a{^E)^-a\, (42) 

where as is the Bohr radius. The rms radius of is < >= 3a^, the radius of an equivalent 
hard sphere is given by i?^ = 5a|, and as a result the hydrogen atom is a pretty good conductor, 
a/volume f=:i 1/10. 

In the following sections we report on the polarizabilities of the nucleon. As compared to 
hydrogen and other atoms, we shall find that the nucleon is a dielectric medium with e ~ 1.001, 
i.e., a very good insulator. Furthermore, magnetic effects are a priori of the same order as 
the electric ones, because the charged constituents, the quarks move with velocities close to 
the velocity of light. However the diamagnetic effects of the pion cloud and the paramagnetic 
effects of the quark core of the nucleon tend to cancel, with the result of a relatively small net 
value of (5. We shall see that "virtual" light allows one to gain information about the spatial 
distribution of the polarization densities, which will be particularly interesting to resolve the 
interfering effects of para- and diamagnetism. Furthermore, the nucleon has a spin and therefore 
appears as an anisotropic object in the case of polarized nucleons. This leads to additional spin 
polarizabilities whose closest parallel in classical physics is the Faraday effect. 



2.2 Real Compton scattering (RCS) : nucleon polarizabilities and the GDH sum rule 



In this section we discuss the forward scattering of a real photon by a nucleon. The incident 
photon is characterized by the Lorentz vectors of momentum, q = (go, q) and polarization. 
Ex = (0, £x), with q - q = (real photon) and Ex - q = (transverse polarization). If the photon 
moves in the direction of the z-axis, q, — qoe^, the two polarization vectors may be taken as 
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(43) 



corresponding to circularly polarized light with helicities A = +1 (right-handed) and A = — 1 
(left-handed). The kinematics of the outgoing photon is then described by the corresponding 
primed quantities. 

For the purpose of dispersion relations we choose the lab frame, and introduce the notation 
g'o* = v for the photon energy in that system. The total cm. energy W is expressed in terms of 
V as : W'^ — + 2Mv, where M is the nucleon mass. The forward Compton amplitude then 
takes the form 



This is the most general expression that is 

(I) constructed from the independent vectors e', e, q' = q (forward scattering!), and cr (the 

proton spin operator), 
(11) linear in e' and e, 

(III) obeying the transverse gauge, • q' = e • q = 0, and 

(IV) invariant under rotational and parity transformations. 

Furthermore, the Compton amplitude has to be invariant under photon crossing, correspond- 
ing to the fact that each graph with emission of the final-state photon after the absorption of 
the incident photon has to be accompanied by a graph with the opposite time order, i.e. absorp- 
tion following emission ("crossed diagram"). This symmetry requires that the amplitude T of 
Eq. (44) be invariant under the transformation and v ^ —v, with the result that / is an 

even and g an odd function. 



These two amplitudes can be determined by scattering circularly polarized photons (e.g., helic- 
ity A = 1) off nucleons polarized along or opposite to the photon momentum q. The former 
situation (Fig. la) leads to an intermediate state with helicity 3/2. Since this requires a total 
spin S > 3/2, the transition can only take place on a correlated 3-quark system. The transi- 
tion of Fig. lb, on the other hand, is helicity conserving and possible for an individual quark, 
and therefore should dominate in the realm of deep inelastic scattering. Denoting the Comp- 
ton scattering amplitudes for the two experiments indicated in Fig. 1 by T3/2 and T1/2, we find 
f{u) — (Ti/2 + T'3/2)/2 and g{i') — {Tif2 — Ts/2)/2. In a similar way we define the total 
absorption cross section as the spin average over the two helicity cross sections. 



= - (0-3/2 + 0-1/2) , (46) 
and the transverse-transverse interference term by the helicity difference. 



T{u, ^ = 0) = e'* ■£f{iy)+icr- (e'* x e) g{u) . 



(44) 




(45) 
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(a) 



(b) 



Fig. 1. Spin and helicity of a double polarization experiment. The arrows denote the spin projections 
on the photon momentum, the arrows — > the momenta of the particles. The spin projection and hehcity 

of the photon is assumed to be A = 1. The spin projection and helicity of the target nucleon N are 
denoted by Sz and h, respectively, and the eigenvalues of the excited system N* by the corresponding 
primed quantities. 

a) Helicity 3/2: Transition N{Sz = 1/2, h = -1/2) ^ N*{Sz = h = 3/2), which changes the helicity 
by 2 units. 

b) Hehcity 1/2: Transition N{Sz = -1/2, h = +1/2) ^ N*{Sz = h = +1/2), which conserves the 
hehcity. 

^TT = 2 ^'^^/^ ~ '^^2) ■ (47) 

The optical theorem expresses the unitarity of the scattering matrix by relating the absorption 
cross sections to the imaginary part of the respective forward scattering amplitude. 



Im f{v) = ^((7i/2(i/) + (73/2(1/)) = ^ ariiy) , 



(48) 



Due to the smallness of the fine structure constant we may neglect all purely electromag- 
netic processes in this context, such as photon scattering to finite angles or electron-positron 
pair production in the Coulomb field of the proton. Instead, we shall consider only the coupling 
of the photon to the hadronic chaimels, which start at the threshold for pion production, i.e., at 
a photon lah energy z/q = m^(l -|- mT^/2M) ^ 150 MeV. We shall return to this point later in 
the context of the GDH integral. 

The total photoabsorption cross section ax is shown in Fig. 2. It clearly exhibits 3 resonance 
structures on top of a strong background. These structures correspond, in order, to concen- 
trations of magnetic dipole strength (Ml) in the region of the A(1232) resonance, electric 
dipole strength (El) near the resonances A^*(1520) and A^*(1535), and electric quadrupole 
{E2) strength near the A^*(1675). Since the absorption cross sections are the input for the 
dispersion integrals, we have to discuss the convergence for large v. For energies above the 
resonance region (u > 1.66 GeV which is equivalent to total cm. energy W > 2 GeV), ax is 
very slowly decreasing and reaches a minimum of about 115 //b around W — 10 GeV. At the 
highest energies, W ~ 200 GeV (corresponding with z/ ^ 2 ■ 10^ GeV), experiments at DESY 
[14] have measured an increase with energy of the form ax ~ VT^^, in accordance with Regge 
parametrizations through a soft pomeron exchange mechanism [15]. Therefore, it can not be 
expected that the unweighted integral over ax converges. 
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Recently, also the helicity difference has been measured. The first measurement was carried 
out at MAMI (Mainz) for photon energies in the range 200 MeV< u < 800 MeV [16,17]. As 
shown in Fig. 3, this difference fluctuates much more strongly than the total cross section gt- 
The threshold region is dominated by S-wave pion production, i.e., intermediate states with spin 
1/2 and, therefore, mostly contributes to the cross section (Ti/2- In the region of the A (1232) 
with spin J = 3/2, both helicity cross sections contribute, but since the transition is essentially 
Ml, we find (Jzi2l<^\i2 ~ 3, and cr^j. becomes large and positive. Figure 3 also shows that a3/2 
dominates the proton photoabsorption cross section in the second and third resonance regions. 
It was in fact one of the early successes of the quark model to predict this fact by a cancellation 
of the convection and spin currents in the case of (J1/2 [23,24]. 

The GDH collaboration has now extended the measurement into the energy range up to 
3 GeV at ELS A (Bonn) [22]. These preliminary data show a small positive value of cr^j. up to 
i/ 2 GeV, with some indication of a cross-over to negative values, as has been predicted from 
an extrapolation of DIS data [25]. This is consistent with the fact that the helicity-conserving 
cross section (T1/2 should dominate in DIS, because an individual quark cannot contribute to 
(T3/2 due to its spin. However, the extrapolation from DIS to real photons should be taken with 
a grain of salt. 

Having studied the behavior of the absorption cross sections, we are now in a position to set 
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Fig. 3. The helicity difference 0-3/2 (?^) — (^i/2{^) for the proton. The calculations include the contribution 
of ttN intermediate states (dashed curve) [18], rjN intermediate state (dotted curve) [19], and the ttttN 
intermediate states (dashed-dotted curve) [20,21]. The total sum of these contributions is shown by the 
full curves. The MAMI data are from Ref. [16,17] and the (preUminary) ELSA data from Ref. [22]. 

up dispersion relations. A generic form starts from a Cauchy integral with contour C shown in 
Fig. 4, 



c 



V' — V — IS 



du , 



(49) 



where u > Q and £ > 0, i.e., in the limit £ — > the singularity approaches a physical point 
dii v' — V > {). The contour is closed in the upper half-plane by a large circle of radius R 
that eventually goes to infinity. Since we want to neglect this contribution eventually, the cross 

sections have to converge for z/ ^ cxo sufficiently well. As we have seen before, this requirement 
is certainly not fulfilled by (Jt(«^), and for this reason we have to subtract the dispersion relation 
for /. If we subtract at = 0, i.e., consider f{v) — /(O), we also remove the nucleon pole terms 
at = 0. The remaining contribution comes from the cuts along the real axis, which may be 
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Im(v) 




Re(v) 



Fig. 4. The contour C for the dispersion integral Eq. (49). The physical point Ues at v+ie and approaches 
the Re u axis in the limit e +0. Singularities lie on the real axis, poles in the s- and n-channel 
contributions with intermediate nucleon states at = 0, and cuts for |z^| > uq due to production of a 
pion or heavier systems (27r, K, etc.). In addition, there occur resonances in the lower half-plane on the 
second Riemann sheet. 

expressed in terms of the discontinuity of Im / across the cut for a contour as shown in Fig. 4 
or simply by an integral over Im / as we approach the axis from above. By use of the crossing 
relation and the optical theorem, the subtracted dispersion integral can then be expressed in 
terms of the cross section. 



Re/(.) = /(0) + |^P / ^f^,d.'. (50) 

I'D 

Though the dispersion integral is clearly dominated by hadronic reactions, the subtraction is 
also necessary for a charged lepton, because the integral over ctt also diverges (logarithmically) 
for a purely electromagnetic process. We note that in a hypothetical world where this integral 
would converge, the charge could be predicted from the absorption cross section. 
For the odd function g{u)'we. may expect the existence of an unsubtracted dispersion relation. 



If the integrals exist, the relations Eq. (50) and (51) can be expanded into a Taylor series at 
the origin, which should converge up to the lowest threshold, u — uq : 
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The expansion coefficients in brackets parametrize the electromagnetic response of the medium, 
e.g., the nucleon. These Taylor series may be compared to the predictions of the low energy the- 
orem (LET) of Low [26], and Gell-Mann and Goldberger [27] who showed that the leading and 
next-to-leading terms of the expansions are fixed by the global properties of the system. These 
properties are the mass M, the charge eejy, and the anomalous magnetic moment {e/2M)KN 
for a particle with spin 1/2 like the nucleon (i.e., Cp = 1, e„ = 0, Kp = 1.79, /«„ = —1.91). The 
predictions of the LET start from the observation that the leading term for z/ ^ is described by 
the Born terms, because these have a pole structure in that limit. If constructed from a Lorentz, 
gauge invariant and crossing symmetrical theory, the leading and next-to-leading order terms 
are completely determined by the Bom terms. 



f{u) = -^ + {a + (3)u' + 0{u'), (54) 

^^'^ = -^^^ + ^'^' + '^^^'^- ^^^^ 

The leading term of the no spin-flip amplitude, /(O), is the Thomson term already familiar from 
nonrelativistic theory^. The term 0{u) vanishes because of crossing symmetry, and only the 
term 0(z^^) contains information on the internal structure (spectrum and excitation strengths) of 
the complex system. In the forward direction this information appears as the sum of the electric 
and magnetic dipole polarizabilities. The higher order terms O^u^) contain contributions of 
dipole retardation and higher multipoles, as will be discussed in Sec. 3.10. By comparing with 
Eq. (52), we can construct all higher coefficients of the low energy expansion (LEX), Eq. (54), 
from moments of the total cross section. In particular we obtain Baldin's sum rule [28,29], 



and from the next term of the expansion a relation for dipole retardation and quadrupole polar- 
izability. In the case of the spin-flip amplitude g, the comparison of Eqs. (53) and (55) yields 
the sum rule of Gerasimov [2], Drell and Heam [3], 



ire K, 



TV _ 

1^0 



^^iisn^^iiisnd,' ^i, (57) 



2M2 

and a value for the forward spin polarizability [27,30], 



oo 



1^0 



^ By comparing with Eq. (31) we see that we have now converted to Heaviside-Lorentz units, i.e., 
Oiem = e^/47r = 1/137 and = c'^/AttM, here and in all following sections. 
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Baldin's sum rule was recently reevaluated in Ref. [13]. These authors determined the integral 
by use of multipole expansions of pion photoproduction in the threshold region, old and new 
total photoabsorption cross sections in the resonance region (200 MeV< v < 2 GeV), and a 
parametrization of the high energy tail containing a logarithmical divergence of ut- The result 
is 



Q;f + /3f = (13.69 ± 0.14) • 10-^ fm^ , 

q;" + /3" = (14.40±0.66) • 10-^fm^ (59) 

for proton and neutron, respectively. 

Due to the v"^ weighting of the integral, the forward spin polarizability of the proton can be 
reasonably well determined by the GDH experiment at MAMI. The contribution of the range 
200 MeV< z/ < 800 MeV is 70 = -[1.87 ± 0.08 (stat) ±0.10 (syst)] ■ 10"^ fm\ the threshold 
region is estimated to yield 0.90- 10^"^ fm"^, and only —0.04-10"'' fm'^ are expected from energies 
above 800 MeV [31]. The total result is 

7^ = [-1.01 ± 0.08 (stat) ± 0.10 (syst)] • 10"^ fm^ . (60) 



We postpone a more detailed discussion of the nucleon's polarizability to Sec. 3.9 where exper- 
imental findings and theoretical predictions are compared to the results of dispersion relations. 

As we have seen above, the GDH sum rule is based on very general principles, Lorentz and 
gauge invariance, unitarity, and on one weak assumption: the convergence of an unsubtracted 
dispersion relation (DR). It is of course impossible to ever prove the existence of such a sum 
rule by experiment. However, the question is legitimate whether or not the anomalous magnetic 
moment on the Ihs of Eq. (57) is approximately obtained by integrating the rhs of that equation 
up to some reasonable energy, say 3 or 50 GeV. The comparison will tell us whether the anoma- 
lous magnetic moment measured as a ground state expectation value, is related to the degrees 
of freedom visible to that energy, or whether it is produced by very short distance and possibly 
still unknown phenomena. Concerning the convergence problem, it is interesting to note that the 
GDH sum rule was recently evaluated in QED [32] for the electron at order and shown to 
agree with the Schwinger correction to the anomalous magnetic moment, i.e., k^, = deml (2vr). 
This also gives the electromagnetic correction to the sum rule for the proton, which is of relative 
order {Ke/ i^nY ~ 10~^. 

The GDH sum rule predicts that the integral on the rhs of Eq. (57) should be Ip — 205 jib for 
the proton. The energy range of the MAMI experiment [17] contributes /p(200 — 800 MeV) = 
[226 ± 5 (stat) ± 12 (syst)] yub. The preliminary results of the GDH experiment at ELS A [22] 
shows a positive contribution in the range of the 3rd resonance region, with a maximum value 
of C3/2 — Ci/2 ~ 100//b, but only very small contributions at the higher energies with a possible 
cross-over to negative values at > 1.8 GeV. At high v, above the resonance region, one usu- 
ally invokes Regge phenomenology to argue that the integral converges [33,34]. In particular, 
for the isovector channel (T1/2 — o"3/2 u"^^^ at large u, with —0.5 ^ cti ^ being the inter- 
cept of the ai(1260) meson Regge trajectory. For the isoscalar channel, Regge theory predicts 
a behavior corresponding to cci ~ —0.5, which is the intercept of the isoscalar /i(1285) and 
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/i(1420) Regge trajectories. However, these assumptions should be tested experimentally. The 
approved experiment SLAC E-159 [35] will measure the helicity difference absorption cross 
section (T3/2 — (J\i2 for protons and neutrons in the photon energy range 5 GeV < < 40 GeV. 
This will be the first measurement of (T3/2 — cri/2 above the resonance region, to test the con- 
vergence of the GDH sum rule and to provide a baseline for our understanding of soft Regge 
physics in the spin-dependent forward Compton amplitude. 

According to the latest MAID analysis [31] the threshold region yields /p(thr - 200 MeV) = 
-27.5 yub, with a sign opposite to the resonance region, because pion S-wave production con- 
tributes to (Ji/2 only. Combining this threshold contribution with the MAMI value (between 
200 and 800 MeV), the MAID analysis from 800 MeV to 1.66 GeV, and including model esti- 
mates for the TTTT, T] and K production channels, one obtains an integral value from threshold to 
1.66 GeV of [31] : 

I^{yV <2 GeV) = [241 ± 5 (stat) ± 12 (syst) ± 7 (model)] [ih . (61) 

The quoted model error is essentially due to uncertainties in the helicity structure of the tttt 
and K channels. Based on Regge extrapolations and fits to DIS, the asymptotic contribution 

{v > 1.66 GeV) has been estimated to be (—26 ±7) /ib in Ref. [25], whereas Ref. [36] estimated 
this to be (— 13±2) fib. We take the average of both estimates to be (— 20±9) fih as a range which 
covers the theoretical uncertainty in the evaluations of this asymptotic contribution. Putting all 
contributions together, the result for the integral / of Eq. (57) is 

Ip = [221 ± 5 (stat) ± 12 (syst) ± 11 (model)] /x6 ^ /p(sum rule) = 204.8 fih , (62) 

where the systematical and model errors of different contributions have been added in quadra- 
ture. Assuming that the size of the high-energy contribution for the estimate of Eq. (62) is 
confirmed by the SLAC E-159 experiment in the near future, one can conclude that the GDH 
sum rule seems to work for the proton. Unfortunately, the experimental situation is much less 
clear in the case of the neutron, for which the sum rule predicts 

4 (sum rule) = 233.2 //b . (63) 

From present knowledge of the pion photoproduction multipoles and models of heavier mass in- 
termediate state, one obtains the estimate /„ = [147 (tt) -I- 55 (tttt) — 6 (77)] /jb^ 196 //b [31], 
from the contributions of the tt, tttt and i] production channels, thereby assuming the same two- 
pion contribution as in the case of the proton. This estimate for /„ falls short of the sum rule 
value by about 15 %. Given the model assumptions and the uncertainties in the present data, 
one can certainly not conclude that the neutron sum rule is violated. Possible sources of the dis- 
crepancy may be a neglect of final state interaction for pion production off the "neutron target" 
deuteron, the helicity structure of two-pion production, or the asymptotic contribution, which 
still remain to be investigated ? We shall return to this point in the following Section when 
discussing so-called generalized GDH integrals for virtual photon scattering. In any case, the 
outcome of the planned experiments of the GDH collaboration [37] for the neutron will be of 
extreme interest. 
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2.3 Forward dispersion relations in doubly virtual Compton scattering (WCS) 



In this section we consider the forward scattering of a virtual photon with space-Uke four- 
momentum q, i.e., — Qo — (l^ — —Q^ < 0. The first stage of this process, the absorption 

of the virtual photon, is related to inclusive electroproduction, e + iV — > e' + A^'+ anything, 
where e(e ) and N{N ) are electrons and nucleons, respectively, in the initial (final) state. The 
kinematics of the electron is traditionally described in the lab frame (rest frame of A^), with E 
and E' the initial and final energy of the electron, respectively, and 9 the scattering angle. This 
defines the kinematical values of the emitted photon in terms of four-momentum transfer Q and 
energy transfer u. 



g2 = 4^^' sin2 y = E - e\ (64) 
2 

and the lab photon momentum |qiab| — \/Q^ + In the cm. frame of the hadronic interme- 
diate state, the four-momentum of the virtual photon is g = (a;, qcm) with 



Mv-Q'' M 
^ = ^ , Qcm = — qiab , (65) 

where W is the total energy in the hadronic cm. frame. We further introduce the Mandelstam 
variable s and the Bjorken variable x. 



02 

s = 2Mv ^M^ -Q^ ^W'^ , X = -f — . (66) 

The virtual photon spectrum is normalized according to Hand's definition [38] by the "equiva- 
lent photon energy". 



K^Ku^v{\-x)^ . (67) 

An alternate choice would be to use Oilman's definition [39], Kq = |qiab|- 

The inclusive inelastic cross section may be written in terms of a virtual photon flux factor 
Fy and four partial cross sections [19], 



rva{u,Q'), (68) 



dQdE' 



(69) 
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with the photon polarization 



1 + 2(1 + z/Vg2) tan2 ^/2 ' ^'^^^ 



and the flux factor 

a,™, E' K 1 



27r2 E g2 i-e 



(71) 



In addition to the transverse cross section ax and a'rpj, of Eqs. (46) and (47), the longitudinal 
polarization of the virtual photon gives rise to a longitudinal cross section ai and a longitudinal- 
transverse interference cr^^. The two spin-flip (interference) cross sections can only be mea- 
sured by a double-polarization experiment, with h = ±1 referring to the two helicity states 
of the (relativistic) electron, and Pz and the components of the target polarization in the 
direction of the virtual photon momentum qiab and perpendicular to that direction in the scat- 
tering plane of the electron. In the following we shall change the sign of the two spin-flip cross 
sections in comparison with Ref. [19], i.e., introduce the sign convention used in DIS, 

(jrprp = —a'rprp aud a LT = —<^'lt • C^^) 
The partial cross sections are related to the quark structure functions as follows [19] : 



"6771 7-, 

(Tr = r 1 

^ MK ^ 



K 



I+YF2 F, 



72 1/ M 



(^LT = 7 {gi + 92) , (73) 

with the ratio 7 = Q/v. The helicity cross sections are then given by 

as/2 = (Fi -91 + l'92) . (74) 

Due to the longitudinal degree of freedom, the virtual photon has a third polarization vector Sq 
in addition to the transverse polarization vectors e± defined in Eq. (43). A convenient definition 
of this four- vector is 



We note at this point that the factor 7 in the denominator of o"l is missing in Ref. [19]. 
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£o=^(|q|,0, 0, go), (75) 
where we have chosen the z-axis in the direction of the photon propagation, 

g = (go, 0, 0, |q|). (76) 
All 3 polarization vectors and the photon momentum are orthogonal (in the Lorentz metrics!), 

em-q = , e*^-e^, = (-1)™ Smm' , for m, m' = 0, ±1 . (77) 
The invariant matrix element for the absorption of a photon with helicity m is 

Mm^Sm- (J) , (78) 

where J is the hadronic transition current, which is gauge invariant, 

q-{J) = qo{p)-ci-{j) = Q. (79) 

Being Lorentz invariant, the matrix element Mm can be evaluated in any system of reference, 
e.g., in the lab frame and by use of Eq. (79), 

Mo ^ ^(|q,ab|(p) - = T^ip) = . (80) 

V iQlabl 

The VVCS amplitude for forward scattering takes the form (as a 2 x 2 matrix in nucleon spinor 
space) : 

+ zo- • (£'* X e) gTTiy. g') - %<T ■ [(£'* - e) x q\gLT{y. Q^) , (81) 

where we have generalized the notation of Eq. (44) to the VVCS case. The optical theorem 
relates the imaginary parts of the 4 amplitudes in Eq. (81) to the 4 partial cross sections of 
inclusive scattering, 

K 

Im gTT{i^, g2) = -^(^tt{i^, g^) , 

Im gLT{y, g') = ^^LT{y. g') ■ (82) 
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We note that products K ax etc. are independent of the choice of K, because they are directly 
proportional to the measured cross section (see Eqs. (68) and (71)). Of course, the natural choice 
at this point would be = Kq = |qiab|> because we expect the photon three-momentum on the 
rhs of Eq. (82). However, we shall later evaluate the cross sections by a multipole decomposi- 
tion in the cm. frame for which K — Kh is the standard choice. 

The imaginary parts of the scattering amplitudes, Eqs. (82), get contributions from both elas- 
tic scattering at z/^ = Q'^ /2M and inelastic processes above pion threshold, for v > vq = m^r + 
(m^ + Q'^)/2M. The elastic contributions can be calculated from the direct and crossed Bom 
diagrams of Fig. 5, where the electromagnetic vertex for the transition 7* (g)+A^(p) N{p+q) 
is given by 

= Fn{Q')Y + MQ')^^'"^ , (83) 
with Fd and Fp the nucleon Dirac and Pauli form factors, respectively. The choice of the elec- 





Fig. 5. Bom diagrams for the doubly virtual Compton scattering (VVCS) process. 

tromagnetic vertex according to Eq. (83) ensures gauge invariance when calculating the Bom 
contribution to the VVCS amplitude, and yields : 



M ) 



/Bom/,, _ (^eniQ I p2 




9lT{-, Q')-^[ FuFp - -r^^^ G,Gm ) . (84) 



The electric {Ge) and magnetic (Gm) Sachs form factors are related to the Dirac (Fd) and 
Pauli (Fp) form factors by 

Ge (Q') = Fd (Q^) - r Fp{Q^) , Gm (Q^) = Fn (Q^) + Fp{Q^) , (85) 
with T = Q'^/AM'^, and are normalized to 
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Ge{0) = Bat , Gm{0) ^ Cn + Kn ^ fJ'N , 



(86) 



where e^, kn, and /j.^ are the charge (in units of e), the anomalous and the total magnetic 
moments (in units of e/2M) of the respective nucleon. We have split the elastic contributions 
of Eq. (84) into a real contribution (terms in Fd and Fp) and a complex contribution (terms 
in G^; and Gm)- The latter terms have a structure like the susceptibility of Eq. (5) and fulfill a 
dispersion relation by themselves. By use of Eqs. (73), (82), and (84), the imaginary parts of 
the Bom amplitudes can be related to the elastic contributions of the quark structure functions 
and to the form factors. 



^Im/^- = Ff^ = lGi,5(l-x), 

~ 2g2 ^2 - -^1 - -Q^G^d^l- X) , 

^ Im g-T -9f-^9t = lGl5il-x), 

^ Im glT = ^ (gf + gt) = ^ GeGm - x) . (87) 
Q Q 

These equations describe the imaginary parts of the scattering amplitudes in the physical region 
at X = 1 or i/ = i/s. The continuation of the amplitudes to negative or complex arguments 
follows from crossing symmetry (see, e.g., Eq. (45)) and analyticity (see Eq. (16)). 

According to Eqs. (54) and (55), the low energy theorem for real photons asserts that the 
leading and next-to-leading order terms in an expansion in u are completely determined by the 
pole singularities of the Bom terms. However, in the case of virtual photons the limit i/ — > has 
to be performed with care [40], because 



lim lim /(i/, Q^) 4^ lim lim fiv, Q^) . (88) 

If we choose Q'^ — Q right away, we reproduce the results of real Compton scattering, Eqs. (54) 
and (55), for f{i/) — frii^, — 0) and g{v) — grrii^: — 0), while fz and gir vanish 
because of the longitudinal currents involved. On the other hand, if we choose Q'^ finite and let 
u go to zero, the result is quite different. In particular 

jr%u,Q' = 0) = -^e% + O{u') (89) 



while 



/Born(^ ^ 0, Q') = ^ tiN{2eN + «iv) + 0{Q') . (90) 



The surprising result is that a long-wave real photon couples to a Dirac (point) particle, while a 
long-wave virtual photon couples only to a particle with an anomalous magnetic moment, i.e.. 
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a particle with internal structure. The inelastic contributions, on the other hand, are independent 
of the order of the limits. 

It is now straightforward to construct the full VVCS amplitudes by dispersion relations in u 
at = const. For the amplitude /r (which is even in v), we shall need a subtracted DR as in 
the case of Eq. (50), 



TT J — I/^ 



(91) 



The integral in Eq. (91) gets contributions from both the elastic cross section (nucleon pole) at 
v' — vb and from the inelastic continuum for v' > vq : 



,2 



+ 



27r2 



V 



Re /t(0, - Re ff\Q, Q^) 
dv' . 



i/'(i/'2 - 1/2) 



(92) 



In the case of i^^ = Knii^, Q"^) = 1^(1 — x), the dispersion integral is of the same form as in 
Eq. (50) except for a factor {1 — x) typical for that choice of K. The pole contribution which 
enters in Eq. (92) can be read off Eq. (84), 



Re ff%u, Q') 



M - 1/2 



(93) 



The function frii^, Q^) — fr^^ij^, Q^), i-C, excluding the nucleon pole term, is continuous in 
u. Therefore, one may perform a low energy expansion in u. 



Re/T(0, Q') -Re/r^O, Q'] 



a 



(94) 



where the term in (^(z/^) generalizes the definition of the sum of electric and magnetic polariz- 
abilities at finite Comparing Eqs. (92) and (94), one obtains the generalization of Baldin's 
sum rule to virtual photons. 



a{Q')+P{Q' 



1 °^ 
2^ I 



K{u, g2) ariiy, Q^) 



V 



dv , 



Xo 

J 2xFi{x, Q^)dx, 



(95) 
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where in the last line we have expressed the integral in terms of the nucleon structure function Fi 
using Eq. (73). The Callan-Gross relation [41] implies that in the limit of large Q"^ the integrand 
2x Fi{x, Q"^) F2{x, Q"^), i.e., the generalized Baldin sum rule measures the second moment 
of Fi and, asymptotically, the first moment of F2. We can also define the resonance contribution 
to a + (3 through the integral 



e IVI f 

aresiQ"") + PresiQ^) = ^ / 2x F,{x, Q^) dx , (96) 

TT J 

Xres 

where Xres corresponds with W = 2 GeV. 

In Fig. 6, we show the dependence of a + P and compare a resonance estimate with the 
evaluation for Q'^ > 1 GeV^ obtained from the DIS structure function Fi, using the MRSTOl 
parametrization [42]. For the resonance estimate we use the MAID model [18] for the one-pion 
channel and include an estimate for the r] and nn channels according to Ref. [19]. One sees that 
at = 0, the one-pion channel alone gives about 85 % of Baldin's sum rule. Including the 
estimate for the i] and vrvr channels, one nearly saturates Baldin's sum rule. Going to larger 
than 1 GeV^, we also show the sum rule estimate of Eq. (96) obtained from DIS by including 
only the range W < 2 GeV. The comparison of this result with the resonance estimate of MAID 
shows that the MAID model nicely reproduces the dependence of ctt for M/^ < 2 GeV. By 
comparing the full DIS estimate with the contribution from W < 2 GeV, one notices that the 
sum rule value for a + P at Q"^ < I GeV^ is mainly saturated by the resonance contribution, 
whereas for > 2 GeV^, the non-resonance contribution (W > 2 GeV ) dominates the 
sum rule. Therefore, around (5^ — 1 — 2 GeV^, a transition occurs from a resonance dominated 
description to a partonic description. Such a transition was already noticed in Refs. [43,44] 
where a resonance estimate for a + P was compared with the DIS estimate, giving qualitatively 
similar results as shown here 

As in the case of /t, also the longitudinal amplitude /l (which is also even in u), should obey 
a subtracted DR : 



Re ^(0, Q') - Re /^(O, Q' 



+ V 

27r2 



V' fz/'2 - Z/2) 



■dv' 



(97) 



where the pole part to /l can be read off Eq. (84), 



Re ft\v. Q') = - , , GliQ') . (98) 

iVi ly 

Analogously to Eq. (94), one may again perform a low energy expansion for the non-pole (or 
inelastic) contribution to the function /l(z/, Q"^), defining a longitudinal polarizability ^^(Q^) 

^ Note however that in Refs. [43,44], the virtual photon flux differs from Eq. (96). 
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Fig. 6. dependence of the polarizability a + P (left) and (a + /3) • Q'^/ (2M) (right) for the proton, 
as given by Eq. (95). The dashed (dashed-dotted) curves represent the MAID estimate [18,19] for the tt 
(tt + r/ + vrvr) channels. The upper solid curve is the evaluation using the DIS structure function Fi [42]. 
The lower solid curve is the evaluation for the resonance region (W < 2 GeV) using the same DIS 
structure function. The solid circle at = corresponds to the Baldin sum rule [13]. 

as the coefficient of the u"^ dependent term. Equation (97) then yields a sum rule for this polar- 
izability : 



1 °° 

2^ / 



K{u, g") aL {ly, Q' 



V 



dv , 



Trge J \ 4M2 





F2 {x, _ 2a; Fi {x, Q^) + x^ F2 {x, Q^) \ .(99) 



In the last line we used Eq. (73) to express in terms of the first moment of Fl = F2 — 
2xFi and the third moment of F2. Comparing Eqs. (95) and (99), one sees that at large Q^, 
where Fi, F2 and Q^F^ are independent (modulo logarithmic scaling violations), the ratio 
+ ~ The quantity 0;^ is therefore a measure of higher twist (i.e. twist-4) matrix 

elements. 

In Fig. 7, we show the Q"^ dependence for and compare the IVIAID model (for the one-pion 
channel) with the DIS evaluation of Eq. (99), using the IvlRSTOl parametrization [42] for F2 
and Fx,. By confronting the full DIS estimate with the contribution from the range W < 2 GeV, 
one first notices that around Q'^ c:^ 1 — 2 GeV^, a transition occurs from a resonance dominated 
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towards a partonic description, as is also seen for a + [3 in Fig. 6. Furthermore, by comparing 
the MAID model with the DIS evaluation of Eq. (99) in the range W < 2 GeV, one notices that, 
in contrast to the case of a + P, the MAID model clearly underestimates q;^. This points to a 
lack of longitudinal strength in the phenomenological model, which is to be addressed in future 
analyses. 

Similar as in the case of real photons (Baldin sum rule !), the generalized polarizabilities 

can be, in principle, constructed directly from the experimental data. However, this requires a 
longitudinal-transverse separation of the cross sections at constant over a large energy range. 



a, (10'^ fin ) 



a, • Q^/(2 Mf (10'"^) 



1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 




0.25 0.5 



J \ \ T~~7'r-r- 

0.75 1 
(GeV^) 



6 ^ 

5 \ 

4 ^ 

3 ^ 

2 ^ 

1 ^ 

^ 



•T-r-l-T 







(GeV^) 



Fig. 7. dependence of the polarizability (left) and ul • / (2M)^ (right) for the proton, as given 
by Eq. (99). The dashed curve represents the MAID estimate [18,19] for the one-pion channel. The upper 
solid curve is the evaluation using the DIS structure functions F2 and Fl [42]. The lower solid curve is 
the evaluation for the resonance region {W < 2 GeV) using the same DIS structure functions. 

We next turn to the sum rules for the spin dependent VVCS amplitudes (see also Ref. [45] 
where a nice review of generalized sum rules for spin dependent nucleon structure functions has 
been given). Assuming an appropriate high-energy behavior, the spin-flip amplitude qtt (which 
is odd in p) satisfies an unsubtracted DR as in Eq. (51), 



Re gTT{y,Q ) = — V / du . (100) 



Assuming that the integral of Eq. (100) converges, one can separate the contributions from the 
elastic cross section at v' — vb and the inelastic continuum for v' > uq, and by use of Eq. (82), 
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one obtains : 



D ^ D pole/ , ^ ^ f K{v' , Q"^) (7ttW , Q"^) , , 

Re^TT(i^, Q ) = grri^, Q ) + J _ ^2 ' (^^1) 



I'D 



where the pole part is given by Eq. (84) as : 

Re g^^iu, Q^) —j^^^ G^M') ■ (102) 

Performing next a low energy expansion (LEX) for the non-pole contribution to grrii^, Q^), we 
obtain : 

Re QTTi^, Q') - Re g^^{v, Q') = Ia{Q') v + 7o(Q') + -(103) 

For the 0{y) term, Eq. (101) yields a generalization of the GDH sum rule : 



' - - ^x2^2(x,g^)^ (104) 



^0 / 
^ 



where the integral 1a (Q^) has been introduced in Ref. [19]. At = 0, one recovers the GDH 
sum rule of Eq. (57) as /yi(0) = — However, it has to be realized that several definitions 
have been given how to generalize the integral to finite [19]. The definition I a of Eq. (104) 
has the advantage that the (arbitrary) factor K in the photon flux disappears (see the discussion 
after Eq. (82)). In other definitions the factor Kjv in Eq. (104) is simply replaced by 1, which 
formally makes the integral look like the GDH integral for real photons, Eq. (57). Unfortunately, 
these integrals now depend on the definition of K (see Eq. (67). In the following we call these 
integrals 1b (Gilman's definition) and Ic (Hand's definition), and refer the reader to Ref. [19] 
for the expressions analogous to Eq. (104) and further details. We will discuss the 0{v) term in 
Eq. (103) and the first moment of g\ in detail further on, and turn first to the 0{y^) term. 

From the 0{ifi') term of Eq. (103), one obtains a generalization of the forward spin polariz- 
ability. 



Xo 



p^AM^ ? f 4M2 1 
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At large Q^, the term proportional to g2 in Eq. (105) can be dropped and 70 is then proportional 
to the third moment of gi. 

In Fig. 8, we show the Q'^ dependence of 70 and compare the resonance estimate from MAID 
to the evaluation with the DIS structure function gi for Q'^ > 1 GeV^. For the structure function 
gi, we use the recent fit performed in [46], which also provides la error bands for this distri- 
bution, allowing us to determine the experimental error on 70, as shown by the shaded bands 
in Fig. 8. At low Q^, one sees that the estimate for the one-pion channel completely dominates 
7o and reproduces well its measured value at Q'^ - 0. Ai Q"^ > 2 GeV^, the MAID model 
(tt + + TTTT channels) is also in good agreement with the DIS evaluation of the 1^ < 2 GeV 
range in the integral Eq. (105) for 70. Furthermore, comparing the full DIS estimate with the 
contribution from the range W < 2 GeV, we once more observe the gradual transition from 
the resonance dominated to the partonic region. Around = 4 GeV^, the W < 2 GeV region 
contributes about 30 % to 70, whereas for a + /3 at the same Q'^, this contribution is below 10 %. 
This difference can be understood by comparing the sum rule Eq. (105) for 70 with Eq. (95) for 
a + p. From this comparison, one notices that the sum rule for 70 invokes one additional power 
of u in the denominator, giving higher weight to the resonance region as compared with a + p. 
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Fig. 8. dependence of the polarizability 70 (left) and 70 • Q^/ (2M)^ (right) for the proton, as given by 
Eq. (105). The dashed (dashed-dotted) curves represent the MAID estimate [18,19] for the vr (Tr + ry+Trvr) 
channels. The upper soUd curve is the evaluation using the DIS structure function [46]. The lower 
solid curve is the evaluation for the resonance region {W < 2 GeV) using the DIS structure function. 
The shaded bands represent the corresponding error estimates as given by Ref. [46]. The soUd circle at 
= corresponds to the evaluation of Eq. (60). 

We next turn to the amplitude gixii^: Q^), which is even in u. Assuming an unsubtracted DR 
exists for the amplitude gzr, it takes the form 
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p ( n2\ -D pole/ n2^ , 1 -D 7 ^'Kju' ,Q^)aLT{jy' ,Q'^) , 

RegLxii^, Q ) = Re ^f^^. (i/, Q ) + J / ;2 _ ^2) ' ^^^^^^ 

where the pole part is given by Eq. (84) as : 



Re ^i^«(^, Q^) = - G^(Q^) Gm(Q^) . (107) 

One sees that for the unsubtracted dispersion integral of Eq. (106) to converge, the cross section 
clt{i^: Q^) should drop/ai'ter than I / u at large v. One can then perform a low energy expansion 
for the non-pole contribution to gLri^-, Q^), as : 

(108) 

where h{Q^) has been introduced in Ref. [19] as 



2M' 



Q 



^0 

/ c^a;{^i(x, g2)+^2(x, g')}. (109) 

For the term of Eq. (108), one obtains a generalized longitudinal-transverse polariz- 

ability. 



e2 4M2 



I dii^ {<,, (x, Q^) + (x, Q^)} . (110) 



This function is finite in the limit Q"^ 0, and can be evaluated safely on the basis of dispersion 
relations. We note that in Ref. [19], the quantity 5q differs by the factor (1 — in the integrand. 
At large g^, 5lt is proportional to the third moment of the transverse spin structure function 
Qj- = Qi + g2.ln this limit, Wandzura and Wilczek [47] have shown that when neglecting 
dynamical (twist-3) quark-gluon correlations, the transverse spin structure function qt can be 
expressed in terms of the twist-2 spin structure function gi as : 

^a^, g') + ^2 (x, g^) = / dy^^^^. (HI) 

i y 
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Recent experimental data from SLAC [48,49] for the spin structure function g2 show that the 
measured value of g2ix, Q^) (in the range 0.02 < x < 0.8 and 1 GeV^ < < 39 QeV^) is 
consistent with the Wandzura-Wilczek (WW) relation of Eq. (111). One can therefore evaluate 
the rhs of Eq. (110), to good approximation, by calculating the third moment of both sides of 
Eq. (111). By changing the integration variables (x, y) {z, y) with x = z ■ y, one obtains : 



] dxx' I dy^-1^^^ = dyy'g,iy,Q'). (112) 

X ^ 

Combining Eqs. (105) and (1 10) with Eq. (112) and using the WW relation, we may relate the 
generalized spin polarizabilities 5lt{Q'^) and 7o(g^), at large : 



5Lr(Q') ^ - 7o(Q') , (113) 

In Fig. 9, the dependence of the polarizability 5lt is shown both for the MAID model (for the 
one-pion channel) and for the DIS evaluation of Eq. (113). Comparing the MAID model with 
the DIS evaluation for the range W < 2 GeV, one notices that the MAID model underestimates 
6lt, similarly as was seen in Fig. 7 for a^^. As the polarizability Slt involves a longitudinal 
amplitude, this may again point to a lack of longitudinal strength in the MAID model. 

In order to construct the VVCS amplitudes which are in one-to-one correspondence with the 
quark structure functions, it is useful to cast Eq. (81) into a co variant form. 



T{i,, Q\ = 0) = e';e, I i-r + Uv, Q-' 



+—U-'-iA{p--^-iAu..Q-) 

p-q \ J \ q^ J 



M3 



qa{p-qsf}-s-qpf3)S2{i^,Q^)^, (114) 



where £0123 = +1, and s° is the nucleon covariant spin vector satisfying s ■ p = 0, = -1. 
With the definition of Eq. (114), all four VVCS amplitudes Ti,T2, 5*1 and ^2 have the same 
dimension of mass. Furthermore, the 4 new structure functions are related to the previously 
introduced VVCS amplitudes of Eq. (81) as follows : 



T,{u,Q^) = fT{u,Q^), (115) 

T2(z/, Q') = ^ {M^, Q') + h{^, Q')) , (116) 

^ Note that for — > 00, one can again neglect the elastic contribution and make the replacement 
C^/o inEq.(llO). 
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Fig. 9. dependence of the polarizability 6lt{Q^) (left) and dLT{Q^)-Q^ /C^M)^ (right) for the proton, 
as given by Eq. (110). The dashed curve represents the MAID estimate for the one-pion chaimel [18]. 
The upper solid curve is the evaluation using the DIS structure function gr = 9i + 92 [46]. The lower 
solid curve is the evaluation for the resonance region {W < 2 GeV) using the DIS structure function. 
The shaded bands represent the corresponding error estimates as given by Ref. [46]. 



(117) 
(118) 



The Born contributions to these functions can be expressed in terms of the form factors by use 
of Eq. (84) as follows : 



T, (''.«)--ir(^^'' + „2_^2+,,G„j , 

sr>. Q') = r— - MFd + Fp) , (119) 

One sees that the pole singularities appearing atu = ±i Q, due to the denominators in Eqs. (116)- 
(118), are actually canceled by a corresponding zero in the numerator of the Bom terms. The 
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imaginary parts of the inelastic contributions follow from Eq. (82), 



^ ^ vM K ( Q \ M ^ ^^^^^ 

In order to cancel the singularities atu = ±i Q, the following relations should be fulfilled if the 
partial cross sections are continued into the complex //-plane: 



<7T {iQ: Q') 

<JTT (^g, g') 



-^L {iQ, g') 

i<yLT {iQ, g^) 



(124) 



These relations can be verified by realizing that the singularities at g^ + z/^ = correspond to 
the Siegert limit, qiab 0, which also implies qcm — 0. Furthermore, all multipoles vanish in 
that limit, except for the (unretarded) dipole amplitudes. In the case of one-pion production as 
presented in Ref. [19], these are the amplitudes £^0+ and Lo+ of the transverse and longitudinal 
currents, respectively, which become equal in the Siegert limit. The relations of Eq. (124) then 
follow straightforwardly. 

We next discuss dispersion relations for the spin dependent amplitudes Si and 5*2. The spin- 
dependent VVCS amplitude Si is even in u, and an unsubtracted DR reads 



ReSi(u, g') = Re^r'' + - V [ 

TT J 



u'lmSiiu', g2) , 
72 ' 



where the pole part Re Sf"^'^ is obtained from Eq. (1 19) as 



(125) 



Re Srii^, Q') ^ -^^^^^ MQ') {Fnm + Fp{Q' 
We can next perform a low-energy expansion for Si{v^ g^) — Si°^^{v^ g^) as 



(126) 



Re>Si(i/, Q'')-^QSf\v, g2) = 



2 Qign 



HQ') + 



2 CXprn \ 1 



M J Q 



+ o{u^)im) 



where the leading term in follows from Eq. (125) as 
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2M2 r 

hiQ') ^ -Q^J 9iix, Q')dx 



which reduces to the GDH sum rule at = 0, as/i(0) = -K^i/A.By using Eqs. (104), (1 10) 
and (128), one can verify that the term in z/^ in Eq. (127) can be expressed in terms of Ji 
and 5lt- 
At large Q'^, hiQ^) has the limit 



with □ : 



2M2 



oo 



(129) 



1 

ri(Q2) = J dxgiix, Q2). (130) 



For the first moment Fi, a next-to-leading order (NLO) QCD fit to all available DIS data for 
Qiigi) has been performed in Ref. [50], yielding the values at = 5 GeV^ : 



F? = 0.118 ±0.004 ±0.007, 
F^ = -0.058 ±0.005 ±0.008 , 
Ff - F^ = 0.176 ± 0.003 ± 0.007 . (131) 

For the isovector combination Ff — F", the Bjorken sum rule [51] predicts : 



Tl-TI ^ -gA = 0.211 ± 0.001 , ^ oo , (132) 

6 

where gA is the axial-vector weak coupling constant. The inclusion of QCD corrections up to 
order yields [52] : 



F? - F^ = 1,, |l - - 3.5833 (^)' - 20.2153 (^)'| • (133) 

When evaluating Eq. (133) using 3 light quark flavors in and fixing as(M|) at 0.114, one 
obtains [50] : 

^ At — > oo , one can replace J^ " f^, because the elastic contribution to Pi vanishes like 
and is therefore highly suppressed. 
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rf-r^* = 0.182 ±0.005, at g^ = 5GeV^ 



(134) 



One sees that the experimental value of Eq. (131) is in good agreement with the Bjorken sum 
rule value of Eq. (134). 

In Fig. 10, we show the dependence of Ji for the proton and compare the MAID esti- 
mate with the DIS evaluation for > 1 GeV^, using the parametrization of Ref. [46] for gi. 
One immediately sees that the integral /f has to undergo a sign change from the large nega- 
tive GDH sum rule value at = to the positive value at large as extracted from DIS. 
Recent data from SLAC [55] and JLab/CLAS [57] cover the intermediate range. In par- 
ticular, the JLab/CLAS data, which extend downwards to ~ 0.15 GeV^, clearly confirm 
this sign change in the sum rule, which occurs around ~ 0.25 GeV^. The resonance esti- 
mate of MAID, including n + rj + nn channels also displays such a sign change. Given some 
uncertainties in the evaluation of the tttt channels and higher continua, the calculation qualita- 
tively reproduces the trend of the data for the 1^ < 2 GeV contribution to /f . At larger Q^, 
one again notices the gradual transition from a resonance dominated to a partonic description. 
For example, at = 2 GeV^, the 11^ < 2 GeV region amounts to only to 20 % of the total 
sum rule value for Jf . To gain an understanding of this gradual transition in the integral If, 
it was proposed in [58] to parametrize the Q'^ dependence through a vector meson dominance 
type model. This model was refined in Refs. [59,60] by adding explicit resonance contributions 
which are important at low as discussed above and lead to the following phenomenological 
parametrization : 



as 



(135) 



where If'^res iQ'^) is the resonance contribution to /f '", Ti^'as the asymptotic values for the first 
moments of gi, and the scale // was assumed to be the vector meson mass [58], i.e., fj, — rup. 
Furthermore, the parameter c**'" in Eq. (135) was chosen as 



1 + 



1 



2M2 



.p,n 
as 



p,n 
res 



(0) 



(136) 



so as to reproduce the sum rule at Q"^ = 0. In Fig. 10, we use the parametrization of Eq. (135), 
but take the recent experimental value of Eq. (131) for Ff Furthermore, we use as input for 
the resonance contribution at the real photon point Ii^resi^) (corresponding with W < 2 GeV) 
the experimental value from Ref. [17] : If ^gg(O) = —0.95 (open diamond in Fig. 10). For the 
Q"^ dependence of the resonance contribution, we take the MAID estimate [18,19] rescaled to 
the experimental value If^resi^) ^t the real photon point. It is seen from Fig. 10 that the resulting 
calculation (shown by the thin solid curve) gives a rather good description of the sign change 
occurring in 7f at ~ 0.25 GeV^. 

The following Fig. 1 1 displays the results for the generalized GDH integral for the neutron, 
as derived from the "^He data of the Hall A Collaboration at JLab [61] and corrected for nuclear 
effects according to the procedure of Ref. [64]. Recently, also the generalized GDH integral for 
the deuteron has been measured by the Clas Collaboration at JLab [62,63], and will provide a 
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Fig. 10. dependence of the integral If (upper panel) and If -Q^/ (2M^) (lower panel) for the proton, 
as given by Eq. (128). The dashed (dashed-dotted) curve represent the MAID estimate [18,19] for the 
TT (tt + ?7 + tttt) channels. The thin solid curve, covering the whole range, is the parametrization of 
Eq. (135) evaluated as described in the text. The upper thin dotted curves are the C(p^) HBChPT results 
of Ref. [53], whereas the lower thick dotted curves are the corresponding 0{p^) relativistic BChPT 
results of Ref. [54]. The upper thick solid curve at > 1.25 GeV^ is the evaluation using the DIS 
structure function gf of Ref. [46]. The lower thick sohd curve in the same range is the evaluation for the 
resonance region {W < 2 GeV) using the DIS structure function. The shaded bands around the thick 
solid curves represent the corresponding error estimates as given by Ref. [46]. The open star at = 
corresponds with the MAMI data [16] combined with the estimate for the non-measured contribution in 
the range W <2 GeV, as given by Eq. (61). The solid star is the total value of Eq. (62), which includes 
the estimate f or > 2 GeV. The SLAC data are from Ref. [55], the HERJVIES data are from Ref. [56], 
and the preliminary JLab/CLAS data are from Ref. [57] (iimer error bars are statistical errors only, outer 
error bars include systematical errors). 



cross-check for the extraction of the generalized GDH integral for the neutron. The comparison 
of the existing neutron data with the MAID results in Fig. 1 1 shows the same problem as already 
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discussed for real photons: The helicity difference in the low-energy region is not properly de- 
scribed by the existing phase shift analyses. However, the strong curvature at ^0.1 GeV^ 
agrees nicely with the predictions. In Fig. 11, we also show the corresponding parametrization 
of Eq. (135) for Ia, by using the value of Eq. (131) for T^^as MAID estimate for the 

resonance contribution. It is seen that the resulting calculation gives a rather good description 
for the generalized GDH integral I a for the neutron. 
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Fig. 11. The generalized GDH integral 7^(Q^) vs. for the neutron. The dashed curve represents the 
MAID estimate [18,19] for the vr channel. The thin solid curve is the parametrization of Eq. (136), using 
MAID to calculate the resonance contribution. The lower thin dotted curve is the Oij)^) HBChPT results 
of Ref. [53], whereas the upper thick dotted curve is the corresponding Oijy') relativistic BChPT results 
of Ref. [54]. The thick solid curve for > 1.25 GeV^ is the evaluation using the DIS structure function 
g", and the shaded band represents the corresponding error estimate as given by Ref. [46]. The ^He data 
of the Hall A Collaboration at JLab [61] are corrected for nuclear effects according to Ref. [64]. The 
HERMES data are from Ref. [56]. For both data sets : inner error bars are statistical errors only, outer 
error bars include systematical errors. The GDH sum rule value is indicated by the star. 

In Figs. 10 and 11, we also present the heavy baryon chiral perturbation theory (HBChPT) 
calculation to of Ref. [53], as well as the relativistic baryon ChPT (relativistic BChPT) 

calculation to C(p^) of Ref. [54]. From the comparison of both the HBChPT and relativistic 
BChPT calculations to the individual proton and neutron generalized GDH integrals, one sees 
that the chiral expansion may only be applied in a very limited range of < 0.05 GeV^. 
This can be understood from the phenomenological calculations discussed above, where it be- 
came obvious that the GDH integrals for proton and neutron at small are dominated by the 
A(1232) resonance contribution. However in the p - n difference, the A(1232) contribution and 
other isospin 3/2 resonances drop out. Therefore, it was noted in Ref. [65] that the HBChPT 
expansion may be applied in a larger range for the difference If — In Fig. 12, we dis- 
play the dependence of the proton - neutron difference If — It is indeed seen that the 
dependence of the ChPT calculations, in particular the HBChPT calculation, is much less 
steep for the p - n difference and follows the phenomenological estimate over a larger range. 
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Therefore this opens up the possibility, as discussed in Ref. [65], to extend the range of 
the ChPT calculation upwards in Q^. On the other hand, the extension of the operator product 
expansion for T{ — to a value around Q'^ ~ 0.5 GeV^ requires the control of higher twist 
terms, which lattice QCD estimates show to be rather small [66]. This may open the possibility 
to bridge the gap between the low and high regimes, at least for this particular observable. 
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Fig. 12. Q2 dependence of the integral if - 7f (left) and (if - ij") • Q'^/{2M'^) (right) for the proton - 
neutron difference, as given by Eq. (128). The dashed curve is the MAID estimate [18] for the one-pion 
channel. The thin soUd curve is the parametrization of Eq. (136), using MAID to calculate the resonance 
contribution. The lower thin dotted curves are the 0{p^) HBChPT results of Ref [53], whereas the upper 
thick dotted curves are the corresponding 0{p^) relativistic BChPT results of Ref. [54]. The upper thick 
solid curve for > 1.25 GeV^ is the evaluation using the DIS structure function — [46], whereas 
the lower thick solid curve is the evaluation for the resonance region (W < 2 GeV) using the DIS 
structure function. The shaded bands represent the corresponding error estimates as given by Ref [46]. 
The SLAC data are from Ref. [55], and the HERMES data are from Ref. [56] (for both data sets : inner 
error bars are statistical errors only, outer error bars include systematical errors). The GDH sum rule 
value is indicated by the star. 

The second spin-dependent VVCS amplitude 5*2 is odd in u, which leads to the unsubtracted 
DR 



ReS, 



TT J v'^ — 


TT J v'^ — 

where the pole part Re 82°^^ is obtained from Eq. (1 19) as : 



diy' 



(137) 
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Re Sf'{u, Q') = Fp{Q') (Fn{Q') + Fp{Q')) , (138) 

Assuming further that the high-energy behavior of 5*2 is given by 

S2ii^, Q^) for 1/ ^ oo , with as < -1 , (139) 

one can also write down an unsubtracted dispersion relation for the amplitude v S2 (which is 
even in u). 



dv 

1^ 



Re t/^s i/, Q )) = - V / 



2^ Ju'HmS2{i^, Q^) , , 



= Re (. 5.)P- + ^ P / ^"f-^^;^ ^ d.' , (140) 

I'D 

where the pole part is obtained from Eq. (1 19) as : 



Re {yS2{y, Q^W~^^Fp{Q^) (F^(g^) + Fp(g^)) . (141) 

If we subtract Eq. (140) from Eq. (137) multiplied by v, we obtain the "superconvergence 
relation" (for any value of Q^), 

00 

= y^ Im52(t/, g^)c^i/, (142) 


i.e., the pole contribution and the inelastic contribution to that integral should cancel. Equa- 
tion (142) is known as the Burkhardt-Cottingham (BC) sum rule [4]. When Eq. (142) is ex- 
pressed in terms of the nucleon structure function 5^2 (2^, Q^), the BC sum rule implies the van- 
ishing of the first moment of g2, i.e.. 



= y cix52(x, g'), (143) 



and the convergence condition of Eq. (139) leads to : 



g2{x, g^) ^ for X ^ , with ^2 > -1 ■ (144) 

Separating the elastic and inelastic contributions in Eq. (143) and using Eq. (87), we may ex- 
press the BC sum rule for any value of g^ as : 
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2/M-2 y 1 



(145) 



Alternatively the BC sum rule can be written in terms of the Sachs form factors and the absorp- 
tion cross sections, i.e., by 



M2 fK{u,Q^ 



z^2 + g2 



V 



1 Gm{Q'')[Gm[Q'')-Ge{Q'')) 
4 1 + r 



(146) 



Performing a low energy expansion for [yS-i) — [yS-if^'^'^, we obtain from Eq. (140) Q : 



Re z/52(z/, g') - Re {vS^iv. Q^W' = 



^2\\ 2 



1 



z/^ + 0(i/^),(147) 



in terms of the integrals /2, /i, /a and spin polarizabilities 70 and 5lt introduced before. For the 
Born contribution, we obtain from Eqs. (119) and (141) that 



Re {US2{U, g2))Bom _ j^^ ^^^^(^^ g2))pole 



Up 



Fp{Q^) (Fd{Q^)+Fp{Q^)){U%) 



yielding 



/Bom(g2) ^ IppiQ') {FniQ') + Fp{Q' 



(149) 



It is interesting to note that the Born contribution to {^'82) leads exactly to the BC sum rule 
value of Eq. (145). Furthermore, the Born contribution also leads to /f°''"(g^) = /5°™(g^) = 
-F|(Q2)/4. 

By comparing Eqs. (109), (128) and (146) one obtains that J3, defined by Eq. (108), can be 
expressed as : 



HQ') = h{Q') + HQ'')- (150) 



^ Note that the relation of Ref. [19], i.e., /;^(0) - /^(O) = MV(2aem) • (7o(0) - 5lt{^)) ensures that 
the v'^ term in VS2 has no singularity at = 0. 
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If the BC sum rule holds at = 0, one obtains 13(0) = eNKN/4:. 

The BC sum rule has been shown to be satisfied in the case of quantum electrodynamics by 
a calculation in lowest order of ctem [67]. In perturbative QCD, the BC sum rule was calculated 
for a quark target to first order in ct^ and also shown to hold [68]. Furthermore, it is interesting 
to note that the validity of the Wandzura-Wilczek relation of Eq. (Ill) for the transverse spin 
structure function gi + g2 implies that the BC sum rule is satisfied. Indeed one directly obtains 
Eq. (143) by integrating Eq. (111). 

For a nucleon target, the BC sum rule has recently been evaluated at small in HBChPT at 
order 0{p'^) [69], and has also been shown to hold to this order. 

In Fig. 13, we show the MAID model prediction for hiQ"^) of the proton and compare it 
with the BC sum rule value. It is obvious from Fig. 13 that at small Q^, the one-pion channel 
nearly saturates the BC sum rule prediction. At intermediate values of Q^, the MAID calcula- 
tion starts to fall short of the sum rule, because the tttt channels and higher continua become 
increasingly important. One also notices that the HBChPT result at order 0{p'^) [69] for the first 
moment of g2 remains close to the phenomenological sum rule evaluation, in the range up to 
~ 0.3 GeV^ . For the higher moments of g2, it was shown in Ref. [69] that the A -contribution 
is very small, so that the moments of g2 seem to be a promising observable to bridge the gap 
between the HBChPT description at the lower Q'^ and the perturbative QCD result at the larger 
Q'^. In the large region, the first moment of g2 was recently evaluated by the E155 Collab- 
oration [49] at = 5 GeV^, and the integral of g2 over the range 0.02 < x < 0.8 was found 
to be —0.044 ± 0.008 ± 0.003. Although this value differs significantly from zero, it does not 
represent a conclusive test of the BC sum rule, because the behavior of g2 is still unknown in 
the small x region, which remains to be explored by future experiments. 

r 0.06 




0.25 0.5 0.75 1 2 4 

(GeV^) (GeV^) 

Fig. 13. dependence of the integral I2 (left) and I2 ■ (5^/(2M^) (right) for the proton, as given by 
Eq. (145). The dashed curve represents the MAID estimate [18] for the one-pion channel. The dotted 
curve IS the HBChPT result at order 0{p^) [69]. The solid curve is the Burkhardt-Cottingham sum 
rule {rhs of Eq. (145)), using the dipole parametrization for and the parametrization for G^^/G^j 
following from the recent JLab data [70,71]. The shaded band represents the evaluation using the recent 
SLAC E155 data for 52 integrated over the range 0.02 < x < 0.8 [49]. 
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3 Dispersion relations in real Compton scattering (RCS) 

3.1 Introduction 

As we have seen in the previous section, forward photon scattering is closely related to prop- 
erties of the excitation spectrum of the probed system. By use of dispersion relations (DRs) it 
becomes possible to set up sum rules on the basis of general principles and to determine certain 
combinations of polarizabilities from the knowledge of the absorption cross sections alone. In 
the following we shall discuss the general case of RCS and set up dispersion relations valid for 
all angles. 




(a) (b) (c) (d) 




(e) (f) (g) 

Fig. 14. Some typical intermediate states contributing to Compton scattering off the nucleon. Upper row: 
The direct (a) and crossed (b) Bom diagrams with intermediate nucleons, a typical resonance excitation 
in the s-channel (c) and its crossed version (d). Lower row: Typical mesonic contributions with photon 
scattering off an intermediate pion (e), the pion pole diagram (f) and a correlated two-pion exchange such 
as the "cr meson" (g). 

Some typical processes contributing to RCS are shown in Fig. 14. When the nucleon is taken 

as a structureless Dirac particle, only the nucleon pole terms contribute. These are diagrams (a) 
and (b) for the s and u channels, respectively. The differential cross section for this situation, 
first obtained by Klein and Nishina [72] in 1929, is shown in Fig. 15. The inclusion of the 
anomalous magnetic moment leads to a far more complicated result corresponding to the Powell 
cross section [73] in Fig. 15. If we add the pion pole term, Fig. 14 (f), the cross section drops 
one third to the original result of Klein and Nishina. This term is, of course, due to the decay 
7r° — > 7 + 7, and therefore directly related to the axial anomaly, derived on general grounds 
as Wess-Zumino-Witten term [74]. The pion pole term is often referred to as triangle anomaly, 
because the vertex 7r77 can be resolved into a triangular quark loop, a diagram not allowed in a 
classical theory and only appearing due to the renormalization process of quantum field theory. 
As we see from the figure, the pion pole term gives a considerable contribution for backward 
scattering, its effect is sometimes included in the backward spin polarizability 7v (the index tt 
stands for 6 = 180° !), though from the standpoint of dispersion relations it should be considered 
as a pole term like the nucleon pole terms. Except for the diagrams (a), (b), and (f), all other and 
higher diagrams in Fig. 14 have no pole structure, but correspond to excited states in s-, u- or 
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t-channel processes. As such they lead to dispersive contributions whose lowest terms are given 
by the six leading polarizabilities of RCS on the nucleon. The result of a calculation taking 
account only of the electric and magnetic dipole polarizabilities is labeled LEX in Fig. 15. It is 
obvious that the low energy expansion is correct only up to about 80 MeV, in a region where 
the "world data" scatter and give only limited information on the polarizabilities. Therefore, the 
analysis of the modem data has been based on dispersion theory whose results are labeled by DR 
in the figure. Clearly the higher order terms become more and more important with increasing 
photon energies, particularly after crossing the pion threshold (seen as a kink at about 150 MeV) 
with a sharp rise if the energy increases further towards the A(1232) resonance. 




5 



20 40 60 80 100 120 140 160 

E^(MeV) 

Fig. 15. Differential cross section for Compton scattering off the proton as a function of the lab photon 
energy and at fixed scattering angle ^lab = 135°. The curves show the full cross section from fixed-i 
subtracted dispersion relations (solid), the Klein-Nishina cross section (small dots), the Powell cross 
section (dashed), the Powell plus vr'' pole cross section (large dots), and the low energy expansion (LEX) 
including also the leading order contributions from the scalar polarizabilities (dashed-dotted). 



3.2 Kinematics 



Assuming invariance under parity, charge conjugation and time reversal symmetry, the gen- 
eral amplitude for Compton scattering can be expressed by 6 independent structure functions 
Ai(i',t),i = 1, 6 [75]. These structure functions depend on two Lorentz invariant variables, 
e.g., v and t as defined in the following. Denoting the momenta of the initial state photon and 
proton by q andp respectively, and with corresponding final state momenta q' andp', the familiar 
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Mandelstam variables are 

s = {q + py, t=iq-qf, u = {q - p'f , (151) 



with the constraint s + t + u — 2M^. The variable u is defined by 

s — u 



The orthogonal coordinates of the Mandelstam plane, u and t, are related to the initial (E^) and 
final (E' ) photon lab energies, and to the lab scattering angle ^lab by 



t = -4:E^E'^ sin^ = -2M{E^ - E'^) , 

^ = ^7+4^ = ^(^7 + ^;)- (153) 

The physical regions of the Mandelstam plane are shown in Fig. 16 by the horizontally hatched 
areas. The vertically hatched areas are the spectral regions discussed in detail in App. A of 
Ref. [76]. The boundaries of the physical regions in the s, u and t channels are determined by 
the zeros of the Kibble function 

t, u) = t{us - M^) = . (154) 



In particular the RCS experiment takes place in the s-channel region, limited by the line t = 
(forward scattering, 9 = 0°) and the lower right part of the hyperbola us = (backward 
scattering, 9 = 180°). The w-channel region is obtained by crossing (z/ — — z/), and the t- 
channel region in the upper part of Fig. 16 corresponds to the process 7 + 7 N + N and 
requires a value oft> 4M^. 



3.3 Invariant amplitudes and nucleon polarizabilities 

The invariant Compton tensor can be constructed as 

Tf, = e': u ip', X^)H^' u{p, Xn) , (155) 



where e and e' are the polarization vectors of the incoming and outgoing photon, respectively, 
as defined in Eq. (43), u and u are the nucleon spinors, and (A'^) are the nucleon helicities 
in the initial (final) states respectively. The Compton tensor if can be built from the four- 
momentum vectors and Dirac matrices as follows [77]: 
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Fig. 16. The Mandelstam plane for real Compton scattering. The physical regions are horizontally 
hatched. The spectral regions (with boundaries bj and bu) are vertically hatched. 



+i 75T5 + i 75^T6, (156) 



where 



QM = 1 (g' _ , P'M = _ i^M ^ (157) 



with £° ^ = —1. The six tensorial objects in Eq. (156) form a complete basis, and the am- 
plitudes Ti, . . .Tq of Prange are scalar functions of u and t containing the nucleon dynamics. 
Unfortunately, the Prange amplitudes have singularities in the forward and backward directions 
leading to linear dependencies at these points (kinematical constraints). L'vov [78] has therefore 
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proposed a different tensor basis, resulting in the set of amplitudes 



= 1 [Ti + T3 + u{T2 + T,)] , A2 = \ [2T, + z/(T2 + T4)] , 

A^ = lJ^u\T.-T,-i{T,-m , A, = j£^j2MT,-i,{T,-T,)], (158) 



These L'vov amplitudes have no kinematical constraints and are symmetrical under crossing, 

A{-u,t) = A{u,t) , t = l,...6. (159) 

In the spirit of dispersion relations we build the invariant amplitudes by adding the pole con- 
tributions of Fig. 14 (a), (b) and (f), and an integral over the spectrum of excited intermediate 
states. Furthermore, we define the polarizabilities by subtracting the nucleon pole contributions 
Af from the amplitudes and introduce the quantities 

Af^(i/,t) = A,(z/,t) - Af{u,t). (160) 



The polarizabilities are related to these functions and their derivatives at the origin of the Man- 
delstam plane, u = t = 0, 

a.^Ar(0,0), a..= (^) , "^'*=(^) ■ ^^^^^ 

For the spin-independent (scalar) polarizabilities a = asi and /3 = (3 mi, one finds the two 
combinations 

ttEi + Z^Mi = -77- (% + fle) , (162) 

o^Ei - Pmi = «i ) (163) 
Zn 

related to forward and backward Compton scattering respectively. The 4 spin-dependent (vec- 
tor) polarizabilities 71 to 74 of Ragusa [79], and the multipole spin polarizabilities 'Jeiei, 
7A/1A/1, 1M1E2, 1E1M2 of Rcf . [80] (scc Scctiou 3.10), are defined by : 



70 = 7l - 72 - 274 = -7B1E1 - 7M1A/1 - 1MIE2 - 1E1M2 = ^4 , (164) 

7i3 = 71 + 273 = -1E1E1 + 1E1M2 = - -r~T7 ('^5 + ^e) , (165) 



^ Alternatively, the polarizabilities can be defined by also subtracting the vr*^ pole contribution in the 
case of the amplitude A2. 



49 



7l4 = 7l - = -7E1B1 - 27M1M1 - 7siiW2 = 



1 



(2 04 + 05 - oe) , 



(166) 



1 



(a2 + 05), (167) 



= 7l + 72 + 274 = -IfElEl + IMIMI + 7Mlf;2 - 7B1M2 = " 



27rM 



where 70 and 7^ are the spin polarizabilities in the forward and backward directions respectively. 
Since the 7r° pole (see Fig. 14 (f)) contributes to A2 only, the combinations 70, 713 and 714 of 
Eqs. (164)-(166) are independent of the pole term, and only the backward spin polarizability 
is affected by this term. 

3.4 RCS data for the proton and extraction of proton polarizabilities 

A pioneering experiment in Compton scattering off the proton was performed by Gol'danski 
et al. [81] in 1960. Their result for the electric polarizability was asi = 9 ± 2, with a large 
uncertainty in the normalization of the cross section giving rise to an additional systematical 
error of ±5. We note that here and in the following all scalar polarizabilities are given in units of 
10^^ fm^. The next effort to determine the polarizabilities is due to the group of Baranov [82]. 
The data were taken with a bremsstrahlung beam with photon energies up to 100 MeV, and 
the polarizabilities were obtained by a fit to the low-energy expansion (LEX). However, such 
energies are outside the range of the LEX. A later reevaluation by use of dispersion relations [83] 
lead to center values of a^i ^ 12 and /3mi ~ —6, far outside the range of Baldin's sum rule and 
more recent results for the magnetic polarizability (3 mi- In any case these findings were much 
to the surprise of everybody, because the spin flip transition from the nucleon to the dominant 
A(1232) resonance was expected to provide a large paramagnetic contribution of order /3para a; 
10. The first modem experiment was performed at Illinois in 1991 [84]. It was done with a 
tagged photon beam, thereby improving the capability to measure absolute cross sections, and 
in the region of energies between 32 and 72 MeV where the LEX was applicable. Unfortunately, 
by the same token the cross sections were small with the consequence of large error bars. The 
experiment was repeated by the Saskatoon-Illinois group at higher energies above [85] and 
below [83] the pion threshold, and evaluated in the framework of dispersion relations with 
much improved results on the polarizabilities. These results were confirmed, within the error 
bars, by the Brookhaven group working with photons produced by laser backscattering from 
a high-energy electron beam [86]. Even more precise data were recently obtained by the A2 
collaboration at MAMl, using the TAPS setup at energies below pion threshold [87]. The results 
of these modem experiments are compiled in Table 1. 

A fit to all modem low-energy data constrained by the sum rule relation uei + Pmi — 13.8 ± 
0.4 leads to the results [87]: 



the errors denoting the statistical, systematical and model-dependent errors, in order. This new 
global average confirms, beyond any doubt, the dominance of the electric polarizability uei and 



aEi = 12.1 ± 0.3(stat) t 0.4(syst) ± 0.3(mod) , 
/3mi = 1-6 ± 0.4(stat) ± 0.4(syst) ± 0.4(mod) , 



(168) 
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Table 1 

Values for the scalar polarizabilities of the proton as obtained from the modern experiments. 



Data set 


Energies 


Angles 


OiEl + PmI 


OiEl — PmI 




(MeV) 


(degree) 


(10-^ fm3) 


(10-^ fm3) 


niinois 1991 [84] 


32-72 


60, 135 


15.8 ±4.5 ±0.1 


11.9 ±5.3 ±0.2 


Saskatoon 1993 [85] 


149-286 


24-135 


12.1 ± 1.7 ±0.9 


7.9 ± 1.4 ±2.0 


Saskatoon 1995 [83] 


70-148 


90, 135 


15.0 ±3.1 ±0.4 


10.8 ± 1.8 ± 1.0 


LEGS 1998 [86] 


33-309 


70-130 


13.23 ± 0.86l° j° 


10.11 ±1.74lJ:i 


MAMI/TAPS 2001 [87] 


55-165 


59-155 


13.1 ±0.6 ±0.8 


10.7 ±0.6 ±0.8 



the tiny value of the magnetic polarizability f^Mi, which has to come about by a cancellation of 
the large paramagnetic contribution of the A^A spin-flip transition with a nearly equally strong 
diamagnetic term. 

Much less is known about the spin polarizabilities of the proton, except for the forward spin 
polarizability 70 = 7^ - 72 - 274 = [-1.01 ± 0.08(stat) ± O.l(syst)] • 10"'' fm'^, which is 
determined by the GDH experiment at MAMI and dispersion relations according to Eq. (60). 
However, the only other combination for which there exists experimental information is the 
backward spin polarizability = 71 + 72 + 274. Dispersive contributions from the s-chaimel 
integral have been found to be positive and in the range of 5 < 7^ (disp) < 10 (here and 
in the following in units of 10~^ fm^). In addition to this dispersive part, a large contribution 
comes from the t-channel 7r° exchange, 77r(7r°-pole) ~ —46.7 (see Eqs. (177)-(179)), giving a 
total result of —42 < 7,^ < —37. These theoretical predictions have been challenged by a first 
experimental value presented by the LEGS group [86] who found from a combined analysis of 
pion photoproduction and Compton scattering: 

-f^ = -27.1 ± 2.2 (stat + syst)l|^ (mod) , (169) 



where the first error combines statistical and systematical uncertainties, and the second one 
represents the model error. However, there is now contradicting evidence from recent MAMI 
data obtained both at low energies [87] and in the region of the A resonance ([88], [89], [90]). 
Though these new results vary somewhat depending on the subsets of data in different energy 
regions and on the input of the underlying dispersion analysis, they are well in the range of the 
expectations from both dispersion theory and chiral perturbation theory. Typical values are: 

f -36.1 ± 2.1 (stat + syst) ± 0.8(mod) [87] 
77r = < (170) 
[ -37.9 ± 0.6(stat + syst) ± 3.5(mod) [88] , 

where the statistical plus systematical error dominates in the low energy region while the model 
dependency gives rise to a large uncertainty for the experiments in the A region. 

The predictions for the other spin and higher order polarizabilities from dispersion analysis 
and ChPT will be compared in the following chapters. Unfortunately, these polarizabilities are 
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all small and can hardly be deduced without dedicated polarization studies. This will require a 
new generation of experiments with polarized beams, polarized targets, and recoil polarimetry. 



3.5 Extraction of neutron polarizabilities 

The experimental situation concerning the polarizabilities of the neutron is still quite unsat- 
isfactory. The electric polarizability can in principle be measured by scattering low energy 
neutrons off the Coulomb field of a heavy nucleus, while the magnetic polarizability (3^^ re- 
mains essentially unconstrained. This technique seemed to be very promising until the begin- 
ning of the 90's, when Schmiedmayer et al. [91] published a value of 0,"%-^ = 12.6 ±1.5 (stat) ± 
2.0 (syst), obtained by the scattering neutrons with energies 50 eV< < 50 keV off a ^osp^, 
target. Shortly later Nikolenko and Popov [92] argued that the errors were underestimated by a 
factor of 5. These findings were confirmed by a similar experiment [93] resulting in = 0±5, 
and by a further analysis of the systematical errors [94] leading to the estimate 7 < < 19. 

The two remaining methods to measure a'^^ are quasi-free Compton scattering off a bound 
neutron, or elastic scattering from the deuteron. The first experiment on quasi-free Compton 
scattering by a neutron bound in the deuteron was performed by Rose et al. [95]. Interpreted in 
conjunction with Baldin's sum rule, the result is < a^^-^ < 14 with a mean value a^^ ^ 10.7. 
The small sensitivity of the experiment follows from the fact that Thomson scattering vanishes 
for the neutron, and therefore also the important interference between the Thomson term and 
the leading non-Bom amplitude (present in the LEX of the proton!) is absent for the neutron. 
It was therefore proposed to measure a'%i — jS^ji at photon energies in the A region and at 
backward angles. Of course, the analysis will strongly depend on final-state interactions and 
two-body currents. The quality of the analysis can be tested, to some extent, by also measuring 
the polarizabilities of the bound proton. Such results obtained by the TAPS Collaboration at 
MAMI were quite promising [96], o^ei - Pmi = 10.3 ± 1.7 (stat + syst) ± 1.1 (mod). The 
experiment was then extended to the neutron by the CATS/SENECA Collaboration [97]. Data 
were collected with both a deuterium and a hydrogen target and analyzed within the framework 
of Levchuk et al. [98] by use of different parametrizations of pion photoproduction multipoles 
and nucleon-nucleon interactions. The agreement between the polarizabilities of free and bound 
protons was again quite satisfactory, and the final result for the (bound) neutron was 

- Pmi = 9.8 ± 3.6 (stat)l^:i (syst) ± 2.2 (mod) . (171) 



The quasi-free scattering cross section obtained at MAMI is in good agreement with an earlier 
datum of a Saskatoon group [99] measured at 247 MeV and ^lab = 135°. From the ratio between 
the neutron and the proton results this group derived a most probable value of a^i — Pmi = 12, 
however with a very large error bar. The comparison between proton and neutron shows that 
there is no significant isovector contribution in the scalar polarizabilities of the nucleon. 

The second type of experiment, 70? — > 70?, has been performed at SAL [100] and at MAX-lab 
[101]. An analysis with the formalism of Ref. [98] gave the results : 
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«bi-/3mi = -4.8 ±3.9 [100], 

= +3.2 ±3.1 [101], (172) 

leading to values compatible with zero. By comparing the two methods to extract neutron polar- 
izabilities from deuteron experiments, we observe a clear tendancy that elastic Compton scat- 
tering leads to smaller values than those extracted from quasi-free scattering, which remains to 
be studied by future investigations. 



3.6 Unsubtracted fixed-t dispersion relations 

The invariant amplitudes Ai are free of kinematical singularities and constraints, and obey 
the crossing symmetry Eq. (159). Assuming further analyticity and an appropriate high-energy 
behavior, these amplitudes fulfill unsubtracted DRs at fixed t, 

ReA,{u,t) = Af{u,t) + -V 7dv' ^"^--^^(^J'^) , (173) 

where Af are the Bom (nucleon pole) contributions as in App. A of Ref. [75], Im^Aj the 
discontinuities across the s-channel cuts of the Compton process and vq — txIt^ + (m^ + 
t/2)/{2M). However, such unsubtracted DRs require that at high energies {v — > oo) the am- 
plitudes \TCLsAi{u,t) drop fast enough so that the integral of Eq. (173) is convergent and the 
contribution from the semi-circle at infinity can be neglected. For real Compton scattering, 
Regge theory predicts the following high-energy behavior for i/ — > oo and fixed t [75]: 

~ iz-^W , (Aa + Ae) ~ ^^^^^ 

where aM{t) < 0.5 (for t < 0) is a meson Regge trajectory, and where ap{t) is the Pomeron 
trajectory which has an intercept ap(0) ?a 1.08. Note that the Pomeron dominates the high 
energy behavior of the combination of ^43 + ^46. From the asymptotic behavior of Eq. (174), it 
follows that for RCS unsubtracted dispersion relations do not exist for the amplitudes Ai and 
A2. The reason for the divergence of the unsubtracted integrals is essentially given by fixed 
poles in the t-channel, notably the exchange of the neutral pion (for A2) and of a somewhat 
fictitious cr-meson (for Ai) with a mass of about 600 MeV and a large width, which models the 
two-pion continuum with the quantum numbers 7 = J = 0. In order to obtain useful results for 
these two amplitudes, L'vov et al. [75] proposed to close the contour of the integral in Eq. (173) 
by a semi-circle of finite radius i/max (instead of the usually assumed infinite radius!) in the 
complex plane, i.e. the real parts of Ai and A2 are calculated from the decomposition 

ReA{u,t) = Af{u,t) + Ar\i^,t) + Ar{i^,t) , (175) 
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with A-"* the s-channel integral from pion threshold uq to a finite upper limit Umax, and an 
'asymptotic contribution' representing the contribution along the finite semi-circle of ra- 
dius Vjnax in the complex plane. In the actual calculations, the s-channel integral is typically 
evaluated up to a maximum photon energy = Vmax — t/{AM) ^ 1.5 GeV, for which the 
imaginary part of the amplitudes can be expressed through unitarity by the meson photoproduc- 
tion amplitudes (mainly Itt and 27r photoproduction) taken from experiment. 

All contributions from higher energies are then absorbed in the asymptotic term, which is 
replaced by a finite number of energy independent poles in the t channel. In particular the 
asymptotic part of Ai is parametrized by the exchange of a scalar particle in the t channel, i.e. 
an effective "cr meson" [75], 

ATM^Alit) = E^TlS^^ (176) 

t — 



where is the a mass, and QaNN and F^^^ are the couplings of the a to nucleons and photons 
respectively. In Ref. [75], the product of the a couplings in the numerator of Eq. (176) is used 
as a fit parameter, which determines the value of a — P through Eq. (163). 
In a similar way, the asymptotic part of A2 is described by the i-channel pole : 

Af{0,t) ^ P-'^^3.NN _ ^^^^^ 
t — mi 



The coupling determined through the 7r° 77 decay as 

r(^°-77) = ^^^oF^^o,,. (178) 

Using r (7r° -f-f) = 1.1 A eV [102], one obtains 

^^°77 = ^ -0.0252 GeV-\ (179) 



where /tt = 92.4 MeV is the pion-decay constant and the sign is in accordance with the 7r°77 
coupling in the chiral limit, given by the Wess-Zumino-Witten effective chiral Lagrangian [74]. 
With the ttNN coupling constant taken from Ref. [103], gl^j^/Aix = 13.73, the product of the 
couplings in Eq. (177) takes the value F^o^^ Q-knn ~ -0.331 GeV^\ leading to a value of -46.7 
for the pion pole contribution to 77r. On the other hand, the effective chiral Lagrangian yields 
the value -43.5 [104]. 

This procedure is relatively safe for A2 because of the dominance of the 7r° pole or trian- 
gle anomaly, which is well established both experimentally and on general grounds as Wess- 
Zumino-Witten term. However, it introduces a considerable model- dependence in the case of 
Ai. Though a mesons have been repeatedly reported in the past, their properties were never 
clearly established. Therefore, this particle should be interpreted as a parametrization of the 
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I = J = part of the two-pion spectrum, which shows up differently in different experiments 
and hence has been reported with varying masses and widths. 



3. 7 Subtracted fixed-t dispersion relations 



As has been stated in the previous section, unsubtracted DRs do not converge for the ampli- 
tudes Ai and A2. Moreover, the amplitude converges only slowly, and in practice has to be 
fixed by Baldin's sum rule. In order to avoid the convergence problems and the phenomenology 
necessary to determine the asymptotic contributions, it was suggested to consider DRs at fixed 
t that are once subtracted at = [76], 



ReA{iy,t) = Af{u,t) + A{Q,t) - Af {Q,t) 



2 2 

TT 



+00 



V / du' 



Im,A(^',0 



(180) 



These subtracted DRs should converge for all 6 invariant amplitudes due to the two additional 
powers of v' in the denominator, and they are essentially saturated by the ttA^ intermediate states 
as will be shown later. In other words, the lesser known contributions of two and more pions as 
well as higher continua are small and may be treated reliably by simple models. 

The price to pay for this alternative is the appearance of the subtraction functions Ai{v = 
0, t), which have to be determined at some small (negative) value of t. We do this by setting up 
once-subtracted DRs, this time in the variable t. 



Ai{0,t) - Af{0,t) = [A{0,0) - Af(0,0) 

+00 



t—pole 



(0,0) 



t 

+ - 

TT 



dt 
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,lmtAi{0,t') t 7 lmtAi{0,t') 
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/ 



dt' 



t' {f - t) 



,(181) 



where A*^^"'^(0, t) represents the contribution of poles in the t channel, in particular of the 7r° 
pole in the case of A2, which is given by Eq. (177). 

To evaluate the dispersion integrals, the imaginary part due to s-channel cuts in Eq. (180) 
is determined, through unitarity relation, from the scattering amplitudes of photoproduction on 
the nucleon. Due to the energy denominator l/z/'(z/'^ — z/^) in the subtracted dispersion inte- 
grals, the most important contribution is from the ttN intermediate states, while mechanisms 
involving more pions or heavier mesons in the intermediate states are largely suppressed. In our 
calculation, we evaluate the nN contribution using the multipole amplitudes from the analy- 
sis of Hanstein, Drechsel and Tiator (HDT) [105] at energies u < 500 MeV and at the higher 
energies up to u ~ 1.5 GeV we take as input the SAID multipoles (SP02K solution) [106]. 
The expansion of Im^Aj into this set of multipoles is truncated at a maximum angular momen- 
tum jmax = ^± 1/2 = 7/2, with the exception of the energy range in the unphysical region 
where we use j^ax = 3/2. The higher partial waves with j > jj^ax + 1 are evaluated analyti- 
cally in the one-pion exchange (OPE) approximation. The relevant formulas to implement the 
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calculation are reported in App. B and C of Ref. [75]. The multipion intermediate states are 
approximated by the inelastic decay channels of the irN resonances. Since a multipole analysis 
is not yet available for the two-pion channel, we assume that this inelastic contribution follows 
the helicity structure of the one-pion photoproduction amplitudes. In this approximation, we 
first calculate the resonant part of the pion photoproduction multipoles using the Breit-Wigner 
parametrization of Ref. [106], which is then scaled by a suitable factor to include the inelastic 
decays of the resonances. It was found, however, that in the subtracted dispersion relation for- 
malism, the sensitivity to the multipion channels is very small and that subtracted dispersion 
relations are essentially saturated at ~ 0.4 GeV. 

The imaginary part in the ^-channel integral from 4m^ — > +00 in Eq. (181) is saturated 
by the possible intermediate states for the t-channel process (see, for example, Fig. 14 (e) and 
(g)), which lead to cuts along the positive t-axis. For values of t below the KK threshold, the 
t-channel discontinuity is essentially saturated by tttt intermediate states. As a consequence, 
the dependence of the subtraction functions on the momentum transfer t can be calculated by 
including the experimental information on the i-channel process through tttt intermediate states 
as 77 — > TTTT — s> NN. In Ref [76], a unitarized amplitude for the 77 tttt subprocess was 
constructed, and a good description of the available data was found. This information is then 
combined with the tttt — > NN amplitudes determined from dispersion theory by analytical con- 
tinuation of ttN scattering amplitudes [8]. In practice, the upper limit of integration along the 
positive-^ cut ist — 0.78 GeV^, corresponding to the highest t value for which the tttt — > NN 
amplitudes are given in Ref. [8]. In App. A, we show in detail how the discontinuities luitAi 
of the invariant amplitudes Ai{i = 1, 6) in the t-channel (77 NN) can be expressed in 
terms of the corresponding 77 tttt and tttt — > NN amplitudes. 

The second integral in Eq. (181) extends from —00 to a = — 2 (m^-|-2Mm^) pa —0.56 GeV^. 
As we are interested in evaluating Eq. (181) for small (negative) values of t (|t| <S |a|), the inte- 
gral from —00 to a will be highly suppressed by the denominator of the subtracted DRs, result- 
ing in a small contribution. This contribution is estimated by saturation with the A-resonance 
and non-resonant vrA^ intermediate states. In particular, we calculate the non-resonant irN con- 
tribution to the Compton amplitudes through unitarity relation from the OPE and nucleon-pole 
pion-photoproduction amplitudes, while we consider the amplitudes corresponding to diagrams 
(c) and (d) of Fig. 14 for the A-resonance excitation . Finally, the corresponding contributions 
to the discontinuities of the invariant amplitudes Ai at z/ = and negative t are obtained by 
analytical continuation in the unphysical region. We estimated the total uncertainty resulting 
from the negative t integral and the two-pion contributions, by calculating the cross sections 
with and without the negative t integral and the two-pion contributions. The total difference 
can be estimated of the order of 3-5 %, as long as we restrict ourselves to the calculation of 
observables up to the A-resonance region. 

Once the t dependence of the subtraction functions Ai{0,t) is determined, the subtraction 
constants /lj(0, 0) remain to be fixed. Although all 6 subtraction constants ai to could be 
used as fit parameters, we shall restrict the fit to the parameters oi and 02, or equivalently to 
CK^i — Pmi and 7^. The subtraction constants 04, 05 and will be calculated through an unsub- 
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traded sum rule, as derived from Eq. (173), 



2 Y , Im,^4,5,6(z^',i = 0) 

04,5,6 = — I av ; 

1^0 



(182) 



The remaining subtraction constant 03, which is related to uei + Pmi through Eq. (162), will 
be fixed through Baldin's sum rule. 



3.8 Hyperbolic (fixed-angle) dispersion relations 



As we have seen in the previous sections, DRs at constant t have the shortcoming that the 
dispersion integrals get contributions from the unphysical region between the boundaries of the 
physical s and u channel regions. Though in principle the integrand in this region can be con- 
structed by extrapolating a partial wave expansion of the Compton amplitudes, the calculation 
is limited in practice to low partial waves. In order to improve the convergence for larger values 
of t, fixed-angle DRs have been proposed [107] and applied to Compton scattering [108,109]. 
In particular for ^lab = 180°, the path of integration runs along the lower boundary of the s- 
channel region (see Fig. 17) from the origin of the Mandelstam plane to infinity ("s-channel 
contribution"), complemented by a path of integration in the upper half -plane ("i-channel con- 
tribution"). The s — u crossing symmetric hyperbolic integration paths are given by : 



a){u 



(183) 



where a is in one-to-one correspondence with the lab and cm. scattering angles as 

2 1 + cos 6'iab 1 + COS 6'cm 

a — —M — , a — —s 



1 — cos 9 



lab 



1 — cos 9c 



(184) 



A few contours corresponding with fixed values of a are shown in Fig. 17. Along such a path at 
fixed a, one can write down a dispersion integral as 



ReAi{s, t, a) = Af{s, t, a) + A* t, a 
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+ - / ds' lmsAi{s',i,a) 
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Fig. 17. Integration paths in the s-channel region of the Mandelstam plane for RCS at fixed lab angle. 

where the discontinuity in the s channel lmsAi{s', i, a) is evaluated along the hyperbola given 
by 

{s' - a){u' - a) = h, s' + i + u' = 2M^ , (186) 

and lmtAi{s, t', a) runs along the path defined by the hyperbola 

{s - a){u - a) ^ b, s + t' + u^2M'^. (187) 

The integrals in Eq. (185) have a similar form as in the case of subtracted DRs (Sec. 3.7) except 
that the individual partial waves are multiplied with different kinematical factors depending on 
the angle. Though the problems of the partial wave expansion are now cured in the lower half 
plane (for t < 0), the integration in the upper half plane (for t > 0) still runs through the un- 
physical region 4m^ < t < 4M^. In the latter region the t-channel partial wave expansion of 
luitAi convergences if -0.594 GeV^ < a < 0, corresponding to 101° < 6'iab < 180°. In fact, 
a simultaneous investigation of the s— and chaimel Lehmaim ellipses leads to the result that 
the convergence of the partial wave expansion is limited at positive z/^ by the spectral function 
bjj{s, t), and at negative z/^ by the left thick line shown in Fig. 18 which follows from the semi- 
major axis of the ellipse of convergence [8]. 

Summing up we find that DRs att — const are perfect for t — 0, i.e., ^lab = 0, and run into 
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Fig. 18. Fixed lab angle integration paths for RCS in the z^^ — t plane. The right thick line corresponds to 
the spectral function t) which determines the semi-minor axis of the elhpse of convergence, while 
the left thick hne follows from the semi-major axis of the ellipse of convergence and gives the boundary 
of convergence at u'^ < 0. 

problems with increasingly negative t values, particularly at backwards angles, while DRs at 
^lab = const are best at 6'iab = 180° and loose accuracy with decreasing angle. Therefore, the 
two techniques nicely complement each other. 

Holstein and Nathan [108] investigated backward DRs (^^lab = 180°) in order to get rigorous 
bounds for the backward scalar polarizability, aEi — /3mi, and to connect the polarizabilities 
of nucleons and pions. From the s-channel integral they found for the proton (ckei — Pmi)s — 
4.8 — 10.8 ± 3.0, where the numbers on the rhs refer to one-pion intermediate states with parity 
change, one-pion states without parity change, and the error due to the unknown multipole struc- 
ture of heavier intermediate states. For the t-channel the result was (a ei — PMi)t = 10.3 — 1.7, 
the first contribution due to S-wave tttt states, the much smaller second number from D-waves. 
Correlations are found to play an important role for the tttt S -waves. While the pion polariz- 
ability increases the Born contribution from 16.1 to 19.1, the hadronic interaction decreases 
this value to 10.3 as shown above. Their total sum of the s- and t-channel contributions is still 
considerably smaller than the nowadays accepted experimental value, but the interesting find- 
ing of this investigation was the partial resolution of the phenomenological a meson by the tttt 
continuum. 

The same technique was later applied by L'vov and Nathan [109] to study the puzzle of the 

large dispersive contributions to the backward spin polarizability as deduced by the LEGS ex- 
periment [86]. Their final result was 7^ = —39.5 ± 2.4, obtained by adding the contributions 
of the tt" pole (—45.0 ± 1.6), the s-channel ttN states (7.31 ± 1.8) and small contributions of 
nnN states as well as rj and rj' in the i-chaimel. 
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Our results for the proton polarizabilities from unsubtracted fixed-angle DRs are presented 
in the following three tables. We choose 6 convenient combinations of the polarizabilities as 
seen in Table 2 where the total results are given. For the reason outlined before, the fixed-angle 
results deteriorate with decreasing angle, and therefore only the range 100° < ^lab < 180° is 
shown. The contributions of the s- and i-channel paths are listed separately in Tables 3 and 4, 
respectively. 

Table 2 

The proton polarizabilities evaluated by unsubtracted fixed-angle dispersion relations at different lab 
scattering angles. See Eqs. (162)-(167) for definitions. The values are given in units of 10~^ fm^ for the 
scalar polarizabilities, and 10~^ fm^ for the spin polarizabilities. 





OlEl — PmI 


OlEl + PmI 


disp 
In 


713 


714 


70 


180° 


10.89 


10.80 


1.19 


4.32 


-2.36 


-1.07 


140° 


10.58 


10.93 


7.62 


4.33 


-2.35 


-1.09 


100° 


9.36 


11.41 


6.93 


4.40 


-2.28 


-1.14 



In general we expect that the forward polarizabilities, a^i + (^mi and 70, are described best 
by forward DRs, i.e., the values at ^lab = 0° (last line of Table 3). In comparing with the 
Baldin sum rule, Eq. (59), it becomes obvious that a^i -I- (3 mi is not yet saturated at the upper 
limit of our integration, z/ = 1.5 GeV. However, 70 is in good agreement with the experimental 
analysis, Eq. (60), because of the better convergence of the integral Eq. (58). The backward 
polarizabilities, on the other hand, should be evaluated by paths along backward angles. From 
Table 3 we find indeed that the s-channel contribution for a^i — (3 mi and 7**^ is pretty stable 
fore,ab>100°. 

Table 3 

The s-channel contribution to the proton polarizabilities of Table 2, first integral on the rhs of Eq. (185). 
Note that the values for 6'iab = 0° are identical with the results of fixed-i dispersion relations at i = 0. 
Units as in Table 2. 



^lab 


OiEl — PmI 


OlEl + PmI 


disp 
In 


713 


714 


70 


180° 


-5.56 


1.52 


1.1\ 


2.75 


-2.70 


-1.07 


140° 


-5.63 


7.65 


I.IQ 


2.75 


-2.70 


-1.09 


100° 


-5.76 


8.13 


7.71 


2.82 


-2.62 


-1.14 


60° 


-5.76 


9.25 


7.94 


3.19 


-2.31 


-1.07 


20° 


-5.49 


11.29 


8.89 


4.01 


-1.68 


-0.82 


0° 


-5.30 


11.94 


9.29 


4.28 


-1.54 


-0.75 



However, the astounding fact is the large t-channel contribution for cxei — (3mi (see Table 4). 
From the last line in Table 3, one sees that the s-channel integral up to Uj^ax = 1.5 GeV yields 
only a small contribution of about -1-3.3 to the electric polarizability. It is remarkable to observe 
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Table 4 

The t-channel contribution to the proton polarizabilities of Table 2, last integral on the rhs of Eq. (185). 
Units as in Table 2. 



^lab 


OlEl — PmI 


"El + PmI 


disp 


713 


714 


70 


180° 


16.46 


3.28 


0.08 


1.57 


0.34 





140° 


16.20 


3.28 


-0.08 


1.57 


0.34 





100° 


15.11 


3.28 


-0.78 


1.57 


0.34 






that the bulk contribution resides at energies beyond 1.5 GeV. The bad convergence of the s- 
channel integral is related to a strong concentration of the spectral strength in the t-channel, 
close to two-pion threshold. This effect is clearly reflected by the large t-channel contribution 
of about +9.9 to uei (see Table 2, first line). The integrand for the i-channel integral is shown 
in Fig. 19 for the tt'^tt" and 7r°7r° channel and for the sum of both channels. The maximum 
of the integral is at t ^ 0.09 GeV^, and displays a long tail reaching out to higher values of 
t. It is obvious that this contribution contains the phenomenological a meson, which has to be 
introduced to describe the data in the framework of unsubtracted DRs at constant t. Our best 
value for aEi — I3mi comes from backward angles with an error estimated from the stability in 
the region 140° < ^lab < 180°, 

CKEi -/9mi = 10.7 ±0.2 . (188) 



In Tables 2-4 only the dispersive contribution to the backward spin polarizability has been listed. 
If we add the large 7r° pole contribution (see Eqs. (177)-(179)), we obtain 

7^ = -38.8 ± 1.8 , (189) 



the largest error being due to the value of the tt" pole contribution. 

According to Table 4, the ^-channel contributions for the remaining combinations 713 and 714 
are very stable while the s-channel results depend on the path of integration. In the case of the 
polarizability 713, the backward value and the forward value agree within 1% (see Table 2 and 
the last line of Table 3). However, the forward value of 714 = —1.54 differs from its backward 
value substantially. Since this polarizability contains the amplitude 04, which in turn is related 
to the forward polarizability 70, we expect that the forward value is more realistic. 

In summary, we obtain the following results : 



713 = 4.30 ± 0.02 , 

714 = -1.95 ±0.41 , 

70 = -0.91 ± 0.16 , (190) 

where the central value and the errors are derived by combining forward (last line of Table 3) 
and backward (first line of Table 2) DRs. 
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Fig. 19. The integrand for the t-channel contribution from S- waves (left panel) and D- waves (right panel) 

to the polarizability a^i — /?mi- Dashed curve: contribution from the vr+vr" channel, dotted curve: 
contribution from the -k^tt^ channel, solid curve: full result, sum of charged and neutral channels. 

In order to improve the convergence, we shall also consider hyperbolic DRs that are once 
subtracted at s = u = M^, 



ReAi{s,t,a)^Af{s,t,a) + [Ai{M\0,a) - Af{M\0,a 



+ 



ds' lmsAi{s', i, a) 



+ [Al-P°^\t) - Al-P^'^O)] + - 
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(191) 



In addition to the better convergence in the s- and t-channel integrals, the subtraction at s = 
u = allows us to pursue a similar strategy as in the case of fixed-t subtracted DRs. We note 
in fact that the subtraction constants in Eq. (191) are again related to the polarizabilities, i.e. 
Oj = [Ai{M'^, 0, a) - Af{M'^, 0, a)] independent of the value of a. 



3.9 Comparison of different dispersion relation approaches to RCS data 



In this subsection we compare the results from fixed-t and hyperbolic DRs, in both their 
subtracted and unsubtracted versions, with some selected experimental data. Fig. 20 shows the 
differential cross section in the low-energy region for various lab angles, obtained at fixed val- 
ues of ctBi, Pmi: and ^t^. The results from subtracted and unsubtracted fixed-f DRs (full and 
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dashed curves) are nearly identical except for extreme backward scattering. Hyperbolic DRs, 
on the other hand, can only be trusted in the backward hemisphere. The unsubtracted version 
(dashed-dotted curve) clearly fails at ^^lab = 107° and above pion threshold, and of course even 
more so at smaller scattering angles. However, it is extremely satisfying that in all other cases 
the four different approaches agree within the experimental error bars. We can therefore con- 
clude that the analysis of the low-energy data is well under control and that quite reliable values 
can be extracted for the polarizabilities, in particular oei, Pmi, and 7,^, which have a large 
influence on the low-energy cross sections. 
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Fig. 20. Differential cross section for Compton scattering off the proton as a function of the lab photon 
energy E-y and at different scattering angles. Full curves: results from fixed-t subtracted DRs, dashed 
curves: fixed-t unsubtracted DRs, dotted curves: hyperbolic subtracted DRs, dashed-dotted curves: hy- 
perbolic unsubtracted DRs. All result are shown for fixed values of a^i + Pmi = 13.8, oei — Pmi = 10, 
and 7^ = —37. The experimental data are from Ref. [87] (full circles), Ref. [84] (diamonds), Ref. [83] 
(triangles), and Ref. [110] (open circles). 

In Table 5, we show the results from subtracted fixed-t DRs for the fit of the polarizabilities 
Pmi, and 7,^ to the modem low-energy data. In particular, we analyzed both the set of re- 
cent data from the TAPS [87] experiment and the full set of low-energy dataofRefs. [83,84,87,110]. 
For the fitting procedure we used the standard minimization by following two different 



The program package MINUIT from the CERNlib was used. 
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strategies: 1) the spin polarizability, 7^, the polarizabiUty difference, a^i — /9mi, and the polar- 
izability sum, aEi + Pmi, were all used as independent free parameters; 2) 7,^ and aEi — Pmi 
were varied in the fit, while aEi+ Pmi was constrained by Baldin's sum rule, using as additional 
datum point the result from the recent re-evaluation of Ref. [87], asi + Pmi — 13.8 ± 0.4. We 
note that the two fit procedures give consistent results within the error bars. Only the value for 
the polarizability sum aEi + 13 mi obtained from the fit is slightly underestimated with respect 
to the expected value from the Baldin sum rule. However, we note that the set of the fitted data 
mainly covers the backward-angle region where cross sections are quite insensitive to the scalar 
polarizability sum and mainly depend on the correlated effects of 7^ and aEi — Pmi- 

Table 5 

The polarizabilities aEi — f^Mi and 7,^ as obtained by fitting the differential cross sections from different 
experiments with fixed-t subtracted DRs: "TAPS" refers to the data from Ref. [87] (65 data points) fitted 
by using the definition of in Eq. (192), and "global" denotes the values of the fit to the set of data from 
Refs. [83,84,87,1 10] (a total of 101 data points) with the x^-function of Eq. (193). The term fixed denotes 
that aEi + Pmi = 13.8 ib 0.4 is included as a constraint, while free indicates that this combination is 
also a free parameter. The first error band is statistical and the second one is systematic. 





OiEl + PmI 


13.8 ± 0.4 (fixed) 


12.6 ± 1.0 ± 1.0 (free) 


TAPS 


OiEl — PmI 


11.2 ±1.2 ±1.9 


11.4 ± 1.3 ± 1.7 




In 


-35.7 ±3.9 ±0.6 


-35.6 ±2.1 ±0.4 




-^red 


82.1/(66 -2) = 1.3 


80.6/(65 - 3) = 1.3 




aEi + Pmi 


13.8 ± 0.4 (fixed) 


13.2 ± 0.9 ± 0.7 (free) 


global 


OLE\ — PmI 


11.3 ±1.1 ±2.7 


11.1 ± 1.1 ±0.8 




In 


-35.9 ± 1.8 ±3.2 


-36.0 ± 1.8 ±3.2 




-^red 


116.0/(102 - 7) = 1.2 


115.7/(101 -8) = 1.2 



In the fit procedure to the TAPS data, we used the standard definition of x^, i.e. 



X 



E 



cxp 



Aa 



(192) 



where (Texp are the experimental and an^,,,, the calculated cross-sections, and Act are the exper- 
imental error bars. In Eq. (192), the experimental errors were estimated according to Ref. [87], 
by adding in quadrature the statistical errors and the "random systematic errors" which were 
estimated, from uncertainties in the experimental geometry and from the statistics of the simu- 
lation, to be equal to ±5% of the measured cross sections. These statistical errors, including the 
random systematic uncertainties, determine the first error bar in the values of the polarizabilities 
reported as TAPS -fit values in Table 5. The second error bar in these fitted values corresponds to 
systematic uncertainties. These systematic error bars were obtained by rescaling the differential 
cross section by ±3%, assuming that the systematic uncertainties in the data are mainly due to 
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Table 6 

The polarizabilities agi — Pmi and as obtained by the fit of Ref . [87] with fixed-i unsubtracted DR 
to the differential cross sections from different experiments : "TAPS" refers to the data from Ref. [87] 
(65 data points) fitted by using the definition of in Eq. (192), and "global" denotes the values of the fit 
to the set of data from Refs. [83,84,87,110] (a total of 101 data points) with the x^-function of Eq. (193). 
The first error band is statistical and the second one is systematic. The terms fixed and free are defined in 
Table 5. 





OiEl + PMl 


(free) 




f\ 7-11 


12 2 + 08^1 4 


TAPS 


/?M1 


0.8 ±0.9 ±0.5 




Ttt 


-35.9 ± 2.3 =F 0.4 




''^red 


80.6/(65 - 3) = 1.3 




OLE\ + /3m1 


(fixed) 




OLE\ 


12.4 ± 0.6 T 0.5 


global 


Pmi 


1.4 ±0.7 ±0.4 






-36.1 ± 2.1 T 0.4 






108.4/(102 - 7) = 1.1 



errors in the normalization of the measured cross sections. 

The global fit to the different data sets of Refs. [83,84,87, 110] was performed by using a differ- 
ent x^-function, namely 



X 



E 



(AT 



+ 



iV 

A(7, 



(193) 



sys 



where is a normalization parameter used to change the normalization for each data set within 
its systematic errors Ao-gys taken equal to ±3% of the measured cross sections. According to 
Refs. [87] and [83], the minimization of this extended x^-function was performed by taking 
the polarizabilities and the normalization constants for each data set as free parameters. The 
resulting uncertainties in the fitted values of the polarizabilities include contributions from both 
the statistical and the systematic errors. The purely statistical contribution to these error bars 
was obtained by fitting the data with fixed values of the normalization constants. On the other 
hand, the net systematic error bars were derived by assuming that the total uncertainty is the 
result of the sum in quadrature of the statistical and systematic contributions. These statistical 
and systematic errors are given by the first and second error bar, respectively, in the values of 
the polarizabilities quoted as global-fit values in Table. 5. 

In Table 6 we also show the results of Ref. [87] from a fit within fixed-t unsubtracted DRs. 
We note that all the analyses, fixed-t subtracted, fixed-t unsubtracted and subtracted hyperbolic 
DRs, are in quite good agreement, giving us confidence that the model dependence of polariz- 
ability extraction is well under control. 
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Fig. 21. Differential cross section for Compton scattering off the proton as a function of the lab photon 
energy and at different scattering angles. Full curves: results from hyperbolic unsubtracted DRs, 
dashed curves: fixed-t subtracted DRs, dotted curves: fixed-t unsubtracted DRs. All results are shown 
for fixed values of aEi + /?mi = 14.05, uei — Pmi = 10, and = —38. The experimental data are 
fromRef. [89]. 



The analysis is more model dependent when we move towards higher energies. As shown in 
Fig. 21, the results of subtracted fixed-t DRs have serious numerical problems at energies above 
the A(1232). The reason for this failure is, however, not the high energy region itself, which is 
strongly suppressed by the denominator in Eq. (180). Quite on the contrary, the denominator 
creates the problem near the lower limit uq of the integral, which extends into the region where 
the amplitude has to be constructed by a continuation of the partial wave series into the unphys- 
ical region (the area between the line s = (M + m-„Y and the s-channel region in Fig. 16). 
Moreover, there are some systematic differences between fixed-t and hyperbolic DRs at these 
higher energies. The fact that the data seem to favor fixed-t DRs is not surprising: The calcula- 
tions are performed with polarizabilities essentially derived by this method. Since these data are 
taken at backward angles, hyperbolic DRs should in fact be quite appropriate. As can be seen 
from the following Fig. 22, subtracted and unsubtracted hyperbolic DRs agree quite nicely, with 
the exception of the lowest scattering angle where the subtraction is necessary. 
Figure 23 shows the angular distribution at an energy somewhat below the A(1232), which turns 
out to be quite sensitive to the backward spin polarizability. Our calculations confirm the finding 
of Ref. [89]: The value of derived by the LEGS collaboration [86] is related to the fact that 
the LEGS and SAL data lie systematically above the recent MAMI results. This difference can 
be partly compensated by a small change of the Afi+ multipole, e.g., a 2 % increase of 
raises the cross section by nearly 10 %. However, the backward-forward asymmetry can not be 
changed that way but requires a strong variation of 77r, in addition of the effect of the (known) 
£'i+/Mi+ ratio. 
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Fig. 22. Differential cross section for Compton scattering off the proton as a function of the lab photon 
energy and at different scattering angles. Full curves: results from hyperbolic subtracted DRs, dashed 
curves: hyperbolic unsubtracted DRs. All results are shown for fixed values of a^i + Pmi = 14.05, 
ctEi — Pmi = 10, and 7,r = —38. The experimental data are from Ref. [89]. 

In order to get new and independent information on the spin polarizabilities, it will be neces- 
sary to perform double polarization experiments. Figure 24 shows the differential cross sections 
for circularly polarized photons and target polarized perpendicular or parallel to the photon 
beam. Both for parallel and perpendicular polarization, a spin-flip of the target proton changes 
the cross section by large factors. The sensitivity to the backward spin polarizability turns out 
to be largest at the higher energies and for circularly polarized photons hitting protons with 
polarizations perpendicular to the photon beam. It is also demonstrated in Fig. 24 that even an 
unreasonably large 20 % decrease of aEi — Pmi can only simulate a change in of about 2-3 
units, making this an ideally suited observable to access 77r. 

In the case of linearly polarized photons, one can access three additional independent observ- 
ables. In particular, we can classify these polarization observables by assuming the xz plane as 
the photon scattering plane, with the quantization axis along the direction of the incoming pho- 
ton momentum, and denoting with $ the angle between the polarization vector of the photon 
and the x axis. With respect to this frame, one can measure the cross sections with the target 
polarized along the x or z direction and the photon polarization at ^ — ±45°, and the cross 
sections with the target polarization perpendicular to the scattering plane and the photon polar- 
ization parallel ($ = 0°) or perpendicular ($ = 90°) to the scattering plane. The results from 
fixed-t subtracted DRs for these observables are displayed in Fig. 25, and show similar sensi- 
tivity to the and aei — Pmi polarizabilities as in the case of double polarization observables 
with circularly polarized photons. 




67 



A LARA 2000 

■ LEGS 97 

□ MAMI 97 

• SAL 93 



^ 300 

a 250 
^ 225 
200 
175 
150 
125 
100 



Fig. 23. The angular distribution at fixed photon lab energy = 285 MeV. The results are displayed 
for aEi + Pmi = 14.05, a^i — Pmi = 10 and different values of the backward spin polarizability. 
The dashed and soUd curves are the results from fixed-f subtracted DR and hyperboUc subtracted DR, 
respectively, for 7,^ = — 27 (pair of upper curves) and 7^ = — 38 (pair of lower curves). The results from 
hyperbohc DRs are shown at backward angles, ^cm > 100°. The experimental data are from Refs. [89] 
(triangles), [86] (full diamonds), [111] (open diamonds), and [85] (circles). 
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Fig. 24. Double polarization differential cross sections for Compton scattering off the proton, with 
circularly polarized photon and target proton polarized along the photon direction (upper panels) or per- 
pendicular to the photon direction and in the plane (lower panels). The thick and thin curves correspond 
to a proton polarization along the positive and negative directions, respectively. The results of the disper- 
sion calculation at fixed-i are for fixed aEi + (3mi = 13.8, fixed a^i — /9mi = 10, and 77r = —32 (full 
curves), = —27 (dashed curves), and 7,r = —37 (dashed-dotted curves). We also show the result for 
ctEi + Pmi = 13.8, aEi — Pmi = 8 and 7,r = —37 (dotted curves). 
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Fig. 25. Double polarization differential cross sections for Compton scattering off the proton, with lin- 
early polarized photon and target proton polarized parallel or perpendicular to the scattering plane. Upper 
panels: target polarized along the x direction in the scattering plane and photon with linear polarization 
at an angle $ = +45° (thick curves) and <1> = —45° (thin curves) with respect to the scattering plane. 
Middle panels: target polarized along the y direction perpendicular to the scattering plane and linearly 
polarized photon parallel (thick curves) and perpendicular (thin curves) to the scattering plane. Lower 
panels: target polarized along the z direction in the scattering plane and photon with linear polarization at 
an angle $ = +45° (thick curves) and $ = —45° (thin curves) with respect to the scattering plane. The 
results of the dispersion calculation at fixed-t are for fixed uei + Pmi = 13.8, fixed a^i — /3mi = 10, 
and = —32 (full curves), = —27 (dashed curves), and = —37 (dashed-dotted curves). We also 
show the result for aEi + Pmi = 13.8, uei — Pmi = 8 and = —37 (dotted curves). 
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3.10 Physics content of the nucleon polarizabilities 

The physical content of the polarizabihties can be visuaUzed best by effective multipole in- 
teractions for the coupling of the electric {E) and magnetic {H) fields of the photon with the 
internal structure of the nucleon [80,112], 

i^eff = -47r (^i? + ^e2"^'0 ' (194) 
n=l 

where the even and odd upper indices refer to scalar and vector polarizabilities, respectively. In 
particular, the lowest scalar polarizabilities are contained in 

-^iff — I OiEi,u ^ + I Pmi,u + ^ aE2 E^j + Y2 ^M2 Hfj . (195) 

The leading term contains the (static!) electric and magnetic dipole polarizabilities, a = oei 
and /3 = (3 mi- hi the subleading term there appear two derivatives of the fields with regard to 

— * — * -, 

either time or space, E = dfE and Eij = ^{ViEj + VjEi) respectively. Applied to a plane 
wave photon, the subleading term is therefore O{io'^) relative to the leading one. The terms in 
aEi,u and /3ei,u are, of course, retardation or dispersive corrections to the respective leading 
order dipole polarizabilities, while aE2 and Pe2 are the electric and magnetic quadrupole po- 
larizabilities. Combining the static dipole polarizabilities with all terms in the sum with time 
derivatives only, we obtain the "dynamical dipole polarizabilities" aEi{uj) and ^Mii^)- The 
terms involving the gradients build up higher polarizabilities, at fourth order the (static) electric 
{aE2) and magnetic {/3m2) quadrupole polarizabilities. In a similar notation the lowest vector 
or spin polarizabilities are defined by 

-1M1E2 Eij aiHj - 7ijiAf2 Hij aiEj , (196) 
^'yEiEi,u d-{E X E) + \'^MiMi,v a -{H X H) 
—lMiE2,v Eij aiHj — ^EiM2,v Hij aiEj 

—'^lE2E2 ^ijk CTiEjiEkl — 2'yM2M2 ^ijk CTiHjiHkl 

+37M2E3 EijkHjk — 37E2M3 HijkEjk , (197) 

\ {ViVjEk + ViVkEj + VjVkEi) - ^ {6i,^Ek + Sjk^Ei + 6ik^Ej) . (198) 

As in the spin-averaged case, four of the terms in the 0{tu'^) polarizabilities are simply dis- 
persive corrections to the 0{u^) expressions. All polarizabilities defined above can be related 



(3) 



(5) 



where 

Eijk 
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to the multipole expansions given in Ref. [80]. In terms of the standard notation of spherical 
tensors, the polarizabilities correspond to the following coupling of electromagnetic transition 
operators: 



dEL ~ [El X El^ , (3ml ~ [Ml x MlT , 
lELEL ~ [El X ^l]' , IMLML ~ [Ml X Mi]S (199) 

lM{L-l)EL ^ [Ml-1 X ElY , lE(L-l)ML ~ [El-1 X Ml^ ■ 

The higher order polarizabilities given above are uniquely defined by the quantities Oj, and 
tti^t of Eq. (161), as discussed in detail in Refs. [80] and [112]. In particular we find for the 
leading terms the relations 



(y-Ei = (ai + 03 + oe), 
otE2 = -f (ai,t + 03,* + oe.t), 



lEiEi — (^2 — ^4 + 2a5 + ae) , 7mimi — 



7M1E2 - 8^ («2 + 04 + as, , _ 

1E2E2 — ('^2,t " ^4,* + 3a5,t + 2a6,t) JM2M2 = 24^ (~'^2,t " ^4,* — SOs^j + 2a6,t), 

7M2B3 = ~T2^ ("2,t + 04,* + Qe,*) , 7B2M3 = ~12^ {~(^2,t + 0'A,t + ^e,*) ■ 

(200) 

where we neglected recoil contributions of 0{M~'^). For details see Ref. [112]. In terms of 
Ragusa's polarizabilities 7^ one has 



Pmi = +^ (oi - 03 - oe) , 
/3m2 = +f (ai,t — a3,t — a6,t) , 

^ ~8;Fm ("^2 + ot4 + 2a5 - ae) 
7E1M2 = 8^ (02 — a4 — ae) 



7B1B1 = -71 - 73 , 7M1M1 = 74 , 
7M1B2 = 72 + 74 , 7B1M2 = 73 ■ 



(201) 



With these definitions we can now complete the expansion of the forward scattering amplitudes, 
Eqs. (54) and (55), to the next order: 



2„2 



e e 



N 



AnM 
2,^2 



+ {cCElil^) + /3mi(J^)) + n i.<^E2{l^) + f3M2{l^)) + C»(l/6) , (202) 



N 



SnM^ 



(203) 



with aEL{y) = oiEL + ttEL,i/i^^ + O^u'^) and similarly for the magnetic terms. In the spin-flip 
amplitude we have defined 



lo{y) = -(lElEliy) + iMlMli^) + lMlE2{y) + 1e1M2{v)) , 
7o(i^) = -{lE2E2{y) + lM2M2{y) + |7M2i!;3('^) + |7E2M3(i^)) ■ 



(204) 
(205) 
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We repeat that all combinations of polarizabilities appearing in the forward direction can be 
evaluated safely on the basis of DRs. However, the individual polarizabilities suffer from the 
non-convergence of the unsubtracted DRs for Ai and A2, and the bad convergence of ^43. In the 
following section we shall compare the predictions of DRs, ChPT and some QCD-motivated 
models amongst each other and with the available experimental data. 
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Fig. 26. The real part of the proton polarizabilities aEi, Pmi (upper panels) and aE2, Pm2 (lower 
panels) as function of the photon cm. energy oj. Full curves: results from fixed-t subtracted dispersion 
relations. Dashed curves: predictions in leading order HBChPT from Ref. [114] for the isoscalar contri- 
bution to the dynamical polarizabilities up to = 170 MeV. The diamonds are the experimental values 
for the dipole static polarizabilities [87] which are used to fit low-energy constants. 



The imaginary parts of the dynamical polarizabilities are determined from the scattering am- 
plitudes of photoproduction on the nucleon by the unitarity relation. If we take into account 
only the contribution from one-pion intermediate states, the unitarity relations take the follow- 
ing simple form [80] : 

^ c 
^ c 

^ c 

Im/3M2(a;) = ^36E(3|Mg|2 + |m£)|2), (206) 

^ c 

where is the pion momentum, uj the photon cm. energy, and Ef^ and M^i^ are pion photo- 
production multipoles which are summed over the different isotopic or charge channels. 
The real part of these amplitudes, calculated both in dispersion theory and HBChPT [1 13,1 14], 
is displayed in Fig. 26. The dynamical polarizabilities allow for a very detailed study of the 
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internal degrees of freedom. For example, aEi and aE2 clearly show cusp effects due to the 
opening of the pion threshold, and Pmi exhibits the A-resonance structure, with the real part 
passing through zero at the resonance position. The HBChPT calculation nicely reproduces the 
results of DRs. 



3.11 DR predictions for nucleon polarizabilities and comparison with theory 



In a nonrelativistic model like the constituent quark model (CQM), the scalar dipole polariz- 
abilities can be expressed by 



aEi = 2aem2^ ^ _^ +^(Q|z^g'^i|0) ' (207) 

n^O n j 

/3mi = 2ae^ E - fi(OM^ + E dim > (208) 



where d = X] dj = ^ eiYi and = I^'i = '^^i ^^e the electric and magnetic dipole oper- 
ators in the cm. frame of the nucleon. For simplicity the quark masses are taken as rrii = |M, 
and the quark charges are in units of e. The terms 0{M~^) in Eqs. (207,208) are retardation 
or recoil terms, which are small corrections in atomic physics but actually quite sizeable for the 
quark dynamics of the nucleon. Clearly the first term on the rhs of both equations is positive, 
because the dipole matrix elements appear squared and the excitation energy En — Eq is pos- 
itive. The higher order terms 0{M~^), however, are positive for aEi but negative for Pmi- In 
the case of the magnetic polarizability, the leading term describes the paramagnetism which is 
essentially due to the spin-flip transition from the nucleon to the A (1232), while the subleading 
term represent Langevin's diamagnetism. The simple CQM with an oscillator potential connects 
the rms radius (r^)^/^ with the oscillator frequency, loq — 3/ (M(r^)), and yields [1 15] 

aEi = (2aem)/(Mc^o') + 0{M-^) . (209) 



Unfortunately, it is not possible to describe both size and excitation energy in this model. If we 
use the proper size, say the electric Sachs radius of the proton, (r^) = (r^)^, aEi is grossly 
overestimated with a value of about 40. On the other hand, the correct excitation energy for the 
dominant dipole mode N* (1520), leads to a value much too small, agi ~ 3.5. Concerning the 
magnetic polarizability, the magnetic dipole transition to the A (1232) yields a large paramag- 
netic value, f3^i 12, which is somewhat reduced by the diamagnetic terms. 

The fact that we underestimate uei if using the excitation energy of the N* (1520) is easily 
understood: The energy denominator in Eq. (207) has been taken to be nearly 600 MeV, while 
electric dipole absorption due to pion S-wave production already takes place at much smaller 
energies. The strong dependence of aEi on the size of the oscillator parameter can of course 
be used to get close to the experimental numbers, and indeed reasonable results were obtained 
using the MIT bag model [116,117]. However, it was also early recognized that no complete 
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picture of the nucleon can emerge without including the pion cloud. In fact, a detailed study of 
the polarizabilities in a chiral quark model showed that for a reasonable quark core radius of 
0.6 fm, the pion cloud contributions are clearly dominant [118]. 

Systematic calculations of pion cloud effects became possible with the development of chiral 
perturbation theory (ChPT), an expansion in the external momenta and the pion mass ("p expan- 
sion"). The first calculation of Compton scattering in that scheme was performed by Bernard, 
Kaiser and MeiBner in 1991 [119]. Keeping only the leading term in l/m^r, they found the 
following simple relation at order (one-loop calculation) 

aEi = IOPmi = /"T/^ = 12.2 , (210) 



in remarkable agreement with experiment. The calculation was later repeated in heavy baryon 
ChPT, which allows for a consistent chiral power counting, and extended to 0{p'^) yield- 
ing [120] 

all = 10.5 ± 2.0 , = 3.5 ± 3.6 . (211) 



The error bars for these values indicate that several low-energy constants were determined by 
resonance saturation, e.g., by putting in phenomenological information about the A (1232) 
resonance. Since this resonance lies close, it may not be justified to "freeze" the degrees of 
freedom of this near-by resonance. It is for this reason that the "small scale expansion" (SSE) 
was proposed which includes the excitation energy of the A (1232) as an additional expansion 
parameter ("e expansion"). Unfortunately, at 0{e^) the "dynamical" A [121,104] increases the 
polarizabilities to values far above the data [104], 

all = 16.4 and /3^i = 9.1 . (212) 



Since large loop corrections are expected at 0{e'^), a calculation to this order might remedy the 
situation. Elsewise, one would have to shift the problem to large contributions of counterterms, 
thus loosing the predictive power. 

The comparison between the predictions for the scalar polarizabilities in heavy baryon ChPT 
and in fixed-t DR is given in Table 7. The differences between the DR analyses of Ref. [112] 
(HDPV) and Ref. [80] (BGLMN) can be explained by different inputs for the one-pion multi- 
poles (in Ref. [80] the solution of SAID-SP97K was used) and different approximations for the 
multipion channels. In Ref. [80], in addition to the parametrization of the resonant contribution 
of the inelastic channels mentioned in Sec. 3.7, the nonresonant contribution to the two-pion 
photoproduction channel was modeled by calculating the OPE diagram of the tt A reac- 

tion. The difference between the data and the model for two-pion photoproduction consisting of 
resonant mechanism plus the OPE diagram for the nonresonant mechanism, was then fitted and 
attributed to a phenomenological, nonresonant 7iV ttA S-wave correction term. The effect 
of the multipion channels can be seen mainly in the sum agi + Pmi, which, within the BGLMN 
analysis, approximately reproduces the value of Baldin's sum rule as given in Eq. (59). 

Furthermore, in Table 7, we also show the results for uei and Pmi obtained in the dressed 
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K-matrix model of Ref. [122], which turn out to be quite close to the DR results. 
Table 7 

Theoretical predictions for scalar polarizabilities of the proton: to 0{p^) in HBChPT [119], to 0{e^) 
in the small scale expansion [104], in the fixed-t dispersion relation analyses of Ref. [112] (HDPV) and 
Ref. [80] (BGLMN), and in the dressed K-matrix model of Ref. [122] (KS). In the DR calculations aci - 
Pmi = 10.0 is used as input. The values are given in units of 10~^ fm^ for the dipole polarizabilities and 
in units of 10~^ fm^ for the quadrupole polarizabilities. 





0(/) 


0{s') 


HDPV 


BGLMN 


KS 


aEi 


13.6 


16.4 


11.0 


11.9 


12.1 


Pmi 


1.4 


9.1 


1.0 


1.9 


2.4 


aE2 


22.1 


26.2 


28.8 


27.5 




Pm2 


-9.5 


-12.3 


-23.7 


-22.4 





The spin polarizabilities were calculated within the HBChPT approach in Ref. [123]. Taking 
out a common factor C = oiemgA/ (47r^/^m^), Ragusa's polarizabilities at 0{p^) read 

r ^ f t , 9a , 9a I , 9a I 9a\ 

{71,72,73,74} = ^{-1 + y + 2 + 24' "2 " 24/ ' ^^^^^ 



where the first term for each 7^ is the contribution of the t-channel 7r° pole term, and the second 
one the dispersive contribution. Clearly, the pole term is the dominant feature except for the 
case of 72. Whether or not the pole term should be included or dropped in the definition of the 
spin polarizabilities is an open discussion, though from the standpoint of DRs the pole terms 
and dispersion integrals are clearly separated. 

In Table 8, we compare the proton results from heavy baryon ChPT, fixed-t DR and hyper- 
bolic DR analyses for the dispersive contribution to the spin polarizabilities of Eq. (200). The 
agreement between the different DR results is quite satisfactory in all cases, and the spread 
among the different DR values can be seen as the best possible error estimate of such calcula- 
tions to date. Let us also notice that the dressed K-matrix model of Ref. [122] also yields values 
which are quite close to the DR analysis, except for jEiAn which comes out much larger in 
absolute value in the dressed K-matrix model and is responsible for the too large and positive 
value obtained in that model for 70 compared to experiment. One also sees from Table 8 that 
the ChPT predictions disagree in some cases, both among each other and with the DR results. 
It is obvious that the reason for these problems deserves further study. 

In the following we shall discuss the forward spin polarizabilities 70. As is obvious from 
Eqs. (201) and (213), the 7r° pole term cancels in the forward direction, and in agreement with 
forward dispersion relations, Eq. (58), only excited intermediate states contribute. Two recent 
calculations of 70 at O(p^) yield the following result [124,125] 

70 = y C (1 - ^(21 + - 2«:„)) = 4.5 - 8.4 = -3.9 . (214) 
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In another independent investigation, Gellas et al. [126] arrived at a value 70 = —1, close to the 
experimental value 70 = (—1.01 ± 0.08 ± 0.10) of Eq. (60). However, the apparent discrepancy 
within the HBChPT calculations is not related with any differences concerning the observables, 
but merely a matter of definition of a polarizability [127,128]. For comparison, the SSE at 0{e^) 
predicts [104] 

70 = 4.6(7r7V) - 2.4(A) - 0.2(7rA) = 2.0 , (215) 

the individual contributions being due to tiN loops, A poles, and ttA loops. 
Table 8 

Theoretical predictions for the dispersive contribution to spin polarizabilities of the proton: to 0{p^) 
in HBChPT [104], to 0(p^) in HBChPT from the two derivations of Refs. [128] and [124], to 0{e^) 
in the small scale expansion [128], in fixed-t dispersion relation analysis of Ref. [112] (HDPV) and 
Ref. [80] (BGLMN), and in our calculation with hyperbolic dispersion relations (HYP. DR) at ^lab = 
180°. Furthermore, the column (KS) gives the results in the dressed K-matrix model of Ref. [122]. The 
values are given in units of 10~^ fm^ for the lower order polarizabilities and in units of 10~^ fm^ for the 
higher order polarizabilities. 

* : Ref. [129] has suggested that a contribution of +2.5 from the A-pole is still missing. 





0(p3) 


0(/) [128] 


0{e^) 


©(p^) [124] 


HDPV 


BGLMN 


HYP DR 


KS 


7E1B1 


-5.7 


-1.8 


-5.4 


-1.4 


-4.3 


-3.4 


-3.8 


-5.0 


7M1M1 


-1.1 


0.4* 


1.4 


3.3 


2.9 


2.7 


2.9 


3.4 




1.1 


0.7 


1.0 


0.2 


-0.01 


0.3 


0.5 


-1.8 


7M1E2 


1.1 


1.8 


1.0 


1.8 


2.1 


1.9 


1.6 


1.1 


70 


4.6 


-1.1 


2.0 


-3.9 


-0.7 


-1.5 


-1.1 


2.4 


disp 


4.6 


3.3 


6.8 


6.3 


9.3 


7.8 


7.8 


11.4 


1E2E2 


-0.4 


0.08 


-0.28 




-0.16 








7M2M2 


-0.03 


0.06 


-0.03 




-0.09 








7E2M3 


0.11 


0.03 


0.11 




0.08 








7M2£;3 


0.11 


0.03 


0.11 




0.06 









The corresponding results for the neutron scalar and spin polarizabilities are shown in Ta- 
bles 9 and 10. In the case of the scalar polarizabilities, the difference between proton and neutron 
dispersive results is quite small. This is in qualitative agreement with the ChPT calculations, in 
which the isovector effects appear only at the fourth order. At this order however, unknown 
low-energy constants enter the scalar polarizabilities. For the lower-order spin polarizabilities 
of the neutron, shown in Table 10, it is amusing to note that the 0{p^) HBChPT predictions of 
Ref. [124] are quite close to the hyperbolic DR results. The higher-order spin polarizabilities of 
the proton and neutron are predicted to be equal, because only the isoscalar contribution to the 
t-channel structure constants ai^t of Eq. (200) was taken into account. 
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Finally, we also like to mention that very recently the first lattice QCD calculations of hadron 
electric and magnetic polarizabilities were reported [130]. While these initial results look quite 
promising, it is too premature to make a quantitative comparison with experiment at this time, 
because the current lattice calculations were performed for rather large values of the quark mass 
(corresponding with pion masses > 500 MeV). However, in the near future such calculations 
can be envisaged for values of down to about 300 MeV. Furthermore, in this range one 
may make use of ChPT results which calculate the dependence of the polarizabilities on m^r. 
This opens up the prospect to extrapolate the lattice results downwards in m^, and bridge the 
gap in between the existing lattice calculations and the chiral limit. Such a study is very 
worthwhile to investigate in a future work. 

Table 9 

Same as Table 7, but for the scalar polarizabilities of the neutron. In the DR calculations aEi — Pmi = 
11.5 is used as input. 





0{p^) 




HDPV 


BGLMN 


KS 


aEi 


13.6 


16.4 


12.3 


13.3 


12.7 


Pmi 


1.4 


9.1 


0.8 


1.8 


1.8 




22.1 


26.2 


28.8 


27.2 




Pm2 


-9.5 


-12.3 


-23.7 


-23.5 





Table 10 

Same as Table 8, but for the spin polarizabilities of the neutron. 

* : Ref. [129] has suggested that a contribution of -h2.5 from the A-pole is still missing. 
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4 Dispersion relations in virtual Compton scattering (V CS) 

4.1 Introduction 



In this section, we discuss dispersion relations for the virtual Compton scattering (VCS) pro- 
cess. In this process, denoted as 7* + p ^ 7 + p, a spacelike virtual photon (7*) interacts with 
a nucleon (we consider a proton in all of the following as experiments are only performed for a 
proton target so far) and a real photon (7) is produced. At low energies, this real photon plays 
the role of an applied quasi-static electromagnetic field, and the VCS process measures the re- 
sponse of the nucleon to this applied field. In the real Compton scattering process discussed in 
Sec. 3, this response is characterized by global nucleon structure constants such as the nucleon 
dipole and higher order polarizabilities. In contrast, for the VCS process, the virtuality of the 
initial photon can be dialed so as to map out the spatial distribution of these nucleon polariz- 
abilities, giving access to so-called generalized polarizabilities (GPs). 

First unpolarized VCS observables have been obtained from experiments at the MAMI ac- 
celerator [131] at a virtuality = 0.33 GeV^, and recently at JLab [132] at higher virtualities, 
1 < <2 GeV^. Both experiments measured two combinations of GPs. Further experimen- 
tal programs are underway at the intermediate energy electron accelerators (MIT-Bates [133], 
MAMI [134], and JLab [135]) to measure both unpolarized and polarized VCS observables. 

VCS experiments at low outgoing photon energies can be analyzed in terms of low-energy 
expansions (LEXs), proposed in Ref. [136]. In the LEX, only the leading term (in the energy 
of the produced real photon) of the response to the quasi-constant electromagnetic field, due to 
the internal structure of the system, is taken into account. This leading term depends linearly on 
the GPs. As the sensitivity of the VCS cross sections to the GPs grows with the photon energy, 
it is advantageous to go to higher photon energies, provided one can keep the theoretical uncer- 
tainties under control when approaching and crossing the pion threshold. The situation can be 
compared to RCS as described in Sec. 3, where it was shown that one uses a dispersion relation 
formalism to extract the polarizabilities at energies above pion threshold, with generally larger 
effects on the observables. 

In this section, we describe the application of a dispersion relation formalism to the VCS 
reaction with the aim to extract GPs from VCS experiments over a larger energy range. We will 
also review the present status and future prospects of VCS experiments and describe the physics 
contained in the GPs. For more details on VCS, see also the reviews of Refs. [137,138]. 



4.2 Kinematics and invariant amplitudes 

The VCS process on the proton is accessed through the ep — > ep7 reaction. In this process, the 
final photon can be emitted either by the proton, which is referred to as the fully virtual Compton 
scattering (FVCS) process, or by the lepton, which is referred to as the Bethe-Heitler (BH) 
process. This is shown graphically in Fig. 27, leading to the amplitude T^^ ''' of the ep — > ep7 
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reaction as the coherent sum of the BH and the FVCS process : 

rpee'i _ rpBH _|_ rpFVCS (216) 

The BH amplitude T^^ is exactly calculable from QED if one knows the nucleon electromag- 
netic form factors. The FVCS amplitude T^"^^^ contains, in the one-photon exchange approx- 
imation, the VCS subprocess 7*p — > 7p. We refer to Ref. [137] where the explicit expression 
of the BH amplitude is given, and where the construction of the FVCS amplitude from the 
1*P ~^ IP process is discussed. 





Fig. 27. (a) FVCS process, (b) BH process. 

We characterize the four- vectors of the virtual (real) photon in the VCS process 7*p — > 7p 
by q iq') respectively, and the four-momenta of initial (final) nucleons by p ip') respectively. In 
the VCS process, the initial photon is spacelike and we denote its virtuality in the usual way by 

= - Q^. Besides Q^, the VCS process can be described by the Mandelstam invariants 

s^{q + p)\ t^{q-q'f, u^{q-p')\ (217) 



with the constraint 

s + i + ii = 2M2-g^ (218) 

where M denotes the nucleon mass. We furthermore introduce the variable v, which changes 
sign under s ^ u crossing : 

where E]!^'' is the virtual photon energy in the lab frame. In the following, we choose Q^, ^ and t 
as the independent variables to describe the VCS process. In Fig. 28, we show the Mandelstam 
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plane for the VCS process at a fixed value of Q"^ = 0.33 GeV^, at which the experiment of 
Ref. [131] was performed. 
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Fig. 28. The Mandelstam plane for virtual Compton scattering at = 0.33 GeV^. The boundaries of 
the physical s-channel region are ^ = 0° and 9 = 180° for > 0, the n-channel region is obtained by 
crossing, u —>■ —v. The curves for 6* = 0° and 9 = 180° intersect at i/ = 0, t = —Q^, which is the point 
where the generalized polarizabilities are defined. 

The VCS helicity amplitudes can be written as 

TvA^;AA^ = -e%(g,A)4*(g',A')M(p',A'^)A<'^'^M(p,A^) , (220) 

where the polarization four- vectors of the virtual (real) photons are denoted by e {e), and their 
helicities by A (A'), with A = 0, ±1 and A' = ±1. The nucleon helicities are Aat, A'^ = ±1/2, 
and n, u are the nucleon spinors. The VCS tensor M.^^ in Eq. (220) can be decomposed into a 
Born (B) and a non-Bom part (NB) : 

M^"" = M^B + -^aTb • (221) 

In the Born process, the virtual photon is absorbed on a nucleon and the intermediate state 
remains a nucleon, whereas the non-Bom process contains all nucleon excitations and meson- 
loop contributions. The separation between Born and non-Born parts is performed in the same 
way as described in Ref. [136], to which we refer for details (see also Ref. [139]). 
One can proceed by parametrizing the VCS tensor of Eq. (221) in terms of 12 independent am- 
plitudes. In Ref. [140], a gauge-invariant tensor basis was found so that the resulting non-Bom 
invariant amplitudes are free of kinematical singularities and constraints, which is an important 
property when setting up a dispersion relation formalism. This tensor takes the form 

M^'' = i: MQ\u,t)pr , (222) 
1=1 
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where the 12 independent tensors pf^ are given in App. B. The corresponding 12 amplitudes fi 
are expressed in terms of the invariants Q"^, v and t. The tensor basis pf^ is chosen such that 
the resulting invariant amplitudes /j are either even or odd under crossing, which leads to the 
following symmetry relations for the /j at the real photon point : 



/,(0,z/,t) = +/,(0,-z/,t), (^ = 1,2,6,11), 
/,(0,z/,t) = -/,(0,-z/,t), (2 = 4,7,9,10), 



(223) 



while the amplitudes /a, /s, /g, /12 do not contribute at this point, because the corresponding 
tensors vanish in the limit — > 0. 

Nucleon crossing combined with charge conjugation provides the more general constraints 
on the ji at arbitrary virtuality : 



1,2,5,6,7,9,11,12), 



fi (q', i^,t) = - fi (q\ -V, t),{% = 3, 4, 8, 10). 



(224) 



When using dispersion relations, it will be convenient to work with 12 amplitudes that are all 
even in v. This is achieved by introducing the amplitudes Fj (i = 1,...,12) as follows : 



Fi (q\ V, i) = /, (g^ V, t), (i = 1, 2, 5, 6, 7, 9, 11, 12) , 

Fi (g^ u,t)=- f\ (g^ u, t), {t = 3, 4, 8, 10) , 



(225) 



satisfying Fi (Q^, —u, t) — Fi (Q^, i', t) for i = 1,...,12. As the non-Bom invariant amplitudes 
fi^fs 10 ~ ^ fo'" u ^ 0, the definition of Eq. (225) ensures that also all the non-Born amplitudes 
Ff^^ (i = 1,...,12) are free from kinematical singularities. The results for the Born amplitudes 
Ff are fisted in App. B of Ref. [141]. 

From Eqs. (223) and (224), one furthermore finds that F7 and Fg vanish at the real photon 
point. Since 4 of the tensors also vanish in the limit — > 0, only the six amplitudes Fi, F2, 
Fi, Fq, FiQ and Fu enter in real Compton scattering (RCS). These 6 amplitudes are related to 
the RCS amplitudes of Eq. (158) by 



-e'^F, 



-e'Ffi 



-e^F 



10 



1 



2M2 



A3 + Ae 



^3 + ^^4 + ^6, 



A, 



4M2 



2M2 
1 

4M2 
1 



2M 



It - 4Af2 
~ 4M2 

[A -A] , 



^4 + A. 
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11 



4M 



A2 ^ A4 + At 



(226) 



where the charge factor — appears explicitly on the Ihs of Eq. (226), because this factor is 
included in the usual definition of the A . 



4.3 Definitions ofnucleon generalized polarizabilities 



The behavior of the non-Bom VCS tensor M^j^^ of Eq. (222) at low energy (q' = |qj^| 0) 
but at arbitrary three-momentum q = |qcm | of the virtual photon, can be parametrized by 6 
generalized polarizabilities (GPs), which will be denoted by piM' l',ml)S(^^ [136,140,142]. In 
this notation, M. {M') refers to the electric [E), magnetic (M) or longitudinal (L) nature of 
the initial (final) photon, L (L') represents the angular momentum of the initial (final) photon, 
and S differentiates between the spin-flip (S = 1) and non spin-flip (S = 0) character of the 
transition at the nucleon side. Assuming that the emitted real photons have low energies, we 
may use the dipole approximation (L' = 1). For a dipole transition in the final state, angular 
momentum and parity conservation lead to 10 GPs [136]. Furthermore, it has been shown [140] 
that nucleon crossing combined with charge conjugation symmetry of the VCS amplitudes pro- 
vide 4 additional constraints among the 10 GPS. A convenient choice for the 6 independent GPs 
appearing in that approximation has been proposed in Ref. [137] : 



p(Li,Li)0(q)^ p(Mi,Mi)0(q)^ (227) 
p(Li,Li)i(q)^ p(Mi,Mi)i(q)^ p(Mi,L2)i(q)^ p(Li,M2)i(q)_ (228) 

We note at this point that the difference between the transverse electric and longitudinal transi- 
tions is of higher order in q, which explains why the electric multipoles can be replaced by the 
longitudinal ones in the above equations. 

In the limit q — > one finds the following relations between the VCS and RCS polarizabili- 
ties [140] : 

p(Li,Li)0(0) = a , P(^i'^i)°(0) = /5 , 

p(Li,M2)i(0) = -If f 73 , P(^^'^2)i(0) = -fl^ (72 + 74) , (229) 
p(^W(0) = , P(^Wi(o)=0. 



In terms of invariants, the limit q' — > at finite three-momentum q of the virtual photon corre- 
sponds to u ^ and t — at finite Q^. One can therefore express the GPs in terms of the 
VCS invariant amplitudes Fj at the point u — 0,t — —Q'^ for finite Q'^, for which we introduce 
the shorthand : 

FiiQ') = F^^ (q\ iy^O,t^ -Q') . (230) 
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The relations between the GPs and the Fi{Q'^) can be found in Ref. [140]. 

Analogously to the sum rules which we discussed in Sec. 3 for the nucleon polarizabilities 
at = 0, we now turn to dispersion relations for the GPs. From the high-energy behavior 
of the amplitudes Fj, it was found in Ref. [141] that the unsubtracted DRs do not exist for the 
amplitudes Fi and F5, but can be written down for the other amplitudes. Therefore, unsubtracted 
DRs for the GPs will hold for those GPs which do not depend on the two amplitudes Fi and F5. 
However, the amplitude F5 can appear in the combination F5 + 4 Fn, because this combination 
has a high-energy behavior leading to a convergent integral [141,143]. Among the six GPs we 
find four combinations that do not depend on Fi and F5 : 



P^'''''^' + ip(*^^.^^«° = -| (^^) Mgo II A + (2 Fe + F,) - F^^ } , (23 1) 
P^''^''^' = ^ ^E^y" {(F5 + F, + 4Fn) +4MF,,}, (232) 

p(Ll,M2)l _ J_piMl,Ml)l = i ( E + M \'/^ Mqo 
V2qo 3\ E J q2 

X {(F5 + i"7 + 4Fii) +4M(2F6 + i"9)}, (233) 

p{Ll,M2)l , V^p(Ml,L2)l ^ 1 / F + M y/^ go 

2 6\ E J 

X {qo[F, + Fr + AFu)+8M^{2F, + Fg)}, (234) 



where F = y + denotes the initial proton cm. energy and qo = M — E the virtual photon 
cm. energy in the limit q' = 0. For small values of q, we observe the relation go ~ — q^/ (2M). 
Furthermore, in the limit q' = 0, the value of is always understood as being = "~ Q'o' 
which we denote by Q"^ for simplicity of the notation. 



4.4 Fixed-t dispersion relations 



With the choice of the tensor basis of App. B, and taking account of the crossing relation 
Eq. (224), the resulting non-Bom VCS invariant amplitudes Fj (z = 1,...,12) are free of kine- 
matical singularities and constraints, and even in u, i.e., FiiCp, z/, t) = Fi{Q'^, —u, t). Assum- 
ing further analyticity and an appropriate high-energy behavior, these amplitudes fulfill unsub- 
tracted dispersion relations with respect to the variable v at fixed t and fixed virtuality Q^, 



As a historical remark we note that dispersion relations have been considered for the first time for the 
virtual Compton scattering process in [144]. This work considered however a different set of amplitudes 
as discussed here. To avoid numerical artefacts due to kinematical singularities we will only consider 
here DRs in the amplitudes Fj which are free from such singularities and constraints. 
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71 J 



2^7^ , z/'Im,F,(Q2,j.',t) 



_ ^2 



I'D 



(235) 



where we explicitly indicate that the Ihs of Eq. (235) represents the non-Bom (NB) parts of the 
amplitudes. In Eq. (235), is defined as in the discussion following Eq. (221), whereas _p?°'^ 
represents the nucleon pole contribution (i.e. energy factors in the numerators are evaluated 
at the pole positionp"^. Furthermore in Eq. (235), Im^Fj are the discontinuities across the s- 
channel cuts of the VCS process, starting at the pion production threshold, which is the first 
inelastic channel, i.e., uq = + {ml + t/2 + Q'^/2)/{2M). 

Besides the absorptive singularities due to physical intermediate states which contribute to 
the rhs of dispersion integrals as Eq. (235), one might wonder if other singularities exist giving 
rise to imaginary parts. Such additional singularities could come from so-called anomalous 
thresholds [145,146], which arise when a hadron is a loosely bound system of other hadronic 
constituents which can go on-shell (such as is the case of a nucleus in terms of its nucleon 
constituents), leading to so-called triangular singularities. It was shown that in the case of strong 
confinement within QCD, the quark-gluon structure of hadrons does not give rise to additional 
anomalous thresholds [147,148], and the quark singularities are turned into hadron singularities 
described through an effective field theory. Therefore, the only anomalous thresholds arise for 
those hadrons which are loosely bound composite systems of other hadrons (e.g., the S particle 
in terms of A and tt). For the nucleon case, such anomalous thresholds are absent, and the 
imaginary parts entering the dispersion integrals of Eq. (235) are calculated from absorptive 
singularities due to nN, nnN, ... physical intermediate states. 

The assumption that unsubtracted dispersion relations as in Eq. (235) hold, requires that at 
high energies (z/ ^ oo at fixed t and fixed Q'^) the amplitudes lrasFi{Q'^,iy,t) {i = 1,...,12) 
drop fast enough so that the integrals of Eq. (235) are convergent and the contribution from 
the semi-circle at infinity can be neglected. The high-energy behavior of the amplitudes Fi 
was investigated in [141] by considering the Regge limit (z/ oo, at fixed t and fixed Q'^) 
of the VCS helicity amplitudes. As mentioned above, it follows from this analysis that for the 
amplitudes Fi and F5, an unsubtracted dispersion integral does not exist, whereas the other ten 
VCS amplitudes can be evaluated through unsubtracted dispersion integrals as in Eq. (235). 

Having specified the VCS invariant amplitudes and their high energy behavior, we are now 
ready to set up the DR formalism. The difference between Bom and pole terms in Eq. (235) 
vanishes for the four combinations of GPs on the Ihs of Eqs. (231 - 234). They can be directly 
evaluated by unsubtracted DRs through the following integrals for the corresponding Fi{Q'^) : 



vr J u' 



We will next discuss in Sec. 4.4.1 how the s-channel dispersion integrals of Eqs. (235) and 
(236) are evaluated. In particular, unitarity will allow us to express the imaginary parts of the 

Note that of the twelve VCS amplitudes Fj, only for the amplitudes Fi, and a combination of F5 and 
Fii there is a difference between the Born and pole parts. 



84 



VCS amplitudes in temis of -kN, ttttN,... intennediate states. Subsequently, we will show in 
Sec. 4.4.2 how to deal with the remaining two VCS invariant amplitudes for which one cannot 
write down unsubtracted DRs. 

4.4.1 s-channel dispersion integrals 

The imaginary parts of the amplitudes Fj in Eq. (235) are obtained through the imaginary part 
of the VCS helicity amplitudes Ta'a^iAAjv defined in Eqs. (220) and (222). The VCS helicity am- 
plitudes can be expressed by the Fi in a straightforward manner, even though the calculation 
is cumbersome. The main difficulty, however, is the inversion of the relation between the two 
sets of amplitudes, i.e., to express the twelve amplitudes Fi in terms of the twelve independent 
helicity amplitudes. This problem has been solved in Refs. [141,149] in two different ways. 
Firstly, the inversion was performed numerically by applying different algorithms. Secondly, 
an explicit analytical inversion was found as detailed in Ref. [149]. The two different methods 
allow us to cross-check the results. 

Having expressed the amplitudes F^ in terms of the helicity amplitudes, the latter are deter- 
mined by using unitarity. Denoting the VCS helicity amplitudes by Tfi, the unitarity relation 
takes the generic form 

2 Im, Tf, = ^{27Ty6\Px - Pi)TifTxi , (237) 

X 

where the sum runs over all possible intermediate states X. Here we are mainly interested in 
VCS through the A(1232)-resonance region. Therefore, we restrict ourselves to the dominant 
contribution by only taking account of the nN s-channel intermediate states. If one wants to 
extend the dispersion formalism to higher energies, the influence of additional channels, like 
the nnN intermediate states has to be addressed. The helicity amplitudes for nN intermediate 
states are expressed in terms of pion photo- and electroproduction multipoles as specified in 
App. C.4 of Ref [141]. The calculations are performed by use of the phenomenological MAID 
analysis [18], which contains both resonant and non-resonant pion production mechanisms. 
This state-of-the-art analysis is based on the existing pion photo- and electroproduction data. 
A direct evaluation of the rhs of Eq. (237) is not possible due to an incomplete coverage of 
the phase space with the present data sets. Therefore, a phenomenological analysis is needed to 
fully calculate the dispersive input. 

4.4.2 Asymptotic parts and dispersive contributions beyond ttN 

To evaluate the VCS amplitudes Fi and F5 in an unsubtracted DR framework, we proceed 
as in the case of RCS [75]. This amounts to perform the unsubtracted dispersion integrals of 
Eq. (235) for Fi and F5 along the real u-axis in the range —Umax < < +i^max', and to close 
the contour by a semi-circle with radius Uj^ax in the upper half of the complex i/-plane, with the 
result 

ReFnQ', ^, t) = Ff\Q\ V, t) - Ff{Q\ u, t) + i^-*(g^ u, t) + F^Q', u, t), (238) 
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for (i = 1,5), where the integral contributions F/"* (for i = 1,5) are. given by 

.,t) = iv T,i.' '-'■■"fw'-"'-'' . (239) 

n J i/'^ — i/^ 



and with the contributions of the semi-circle of radius Umax identified with the asymptotic con- 
tributions (Fl"^ F^'). 

Evidently, the separation between asymptotic and integral contributions in Eq. (238) is spec- 
ified by the value of Umax- The total result for Ff^ is formally independent of the specific 
value of Umax- In practice, however, Umax is chosen to be not too large so that one can eval- 
uate the dispersive integrals of Eq. (239) from threshold up to Umax sufficiently accurate. As 
we are mainly interested here in a description of VCS up to A(1232)-resonance energies, the 
dispersion integrals are saturated by their ttN contribution and we choose Umax = 1-5 GeV. 
In the following, we denote this contribution by F[^ . Furthermore, the remainder is estimated 
by an energy-independent function, which parametrizes the asymptotic contribution F"'' due to 
t-channel poles, and which contains all dispersive contributions beyond the value u^ax = 1-5 
GeV. We will next discuss the asymptotic contributions F^^ and F"''. 



• The asymptotic contribution F^ 



The asymptotic contribution to the amplitude F5 predominantly results from the t-channel 
-exchange |33 : 



Fr(Q^z.,t) ^ Ff{Q\t) = -AF-l{Q\t) = _A_JlEIlJkT^, (24O) 



For the Q^-dependence of F^-o^^ (Q^)^ one can use the interpolation formula proposed in 
[150] : 

F^o^^(O) 



where F^o^^(O) has been given in Eq. (179). Equation (241) provides a rather good parametriza- 
tion of the 7r°7*7 form factor data over the whole Q"^ range. 

When fixing the asymptotic contribution F^^ through its 7r°-pole contribution as in Eq. (240), 
one can determine one more GP of the nucleon, in addition to the 4 combinations of Eqs. (23 1 ) 
- (234). In particular, the GP p(^^i'^^i)i which contains F5 can be calculated through 

p(M,Mi)i ^ _V2 i^F + M^ V2 ^^^^^^^ ^ ^,F,m} . (242) 



As mentioned before, the vr^-pole only contributes to the amplitudes F5 and Fn, but drops out in the 
combination (F5 + 4Fii), which therefore has a different high-energy behavior. 
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• The asymptotic part and dispersive contributions beyond tiN to Fi 



We next turn to the high-energy contribution to Fi. As in the case of RCS, the asymptotic 
contribution to the amplitude Fi originates predominantly from the i-channel tttt intermediate 
states. In a phenomenological analysis, this continuum is parametrized through the exchange 
of a scalar-isoscalar particle in the t-channel, i.e. an effective "cr"-meson, as suggested in [75] 
and discussed in Sec. 3.6 for the RCS case. In this spirit, the difference between F^^ and its 
T^N contribution can be parametrized at = by the energy-independent function : 

Fr(0, V, t) - F-^(0, V, t) ^ [Fr(0, 0, 0) - F-^(0, 0, 0)1 ] ^ , (243) 

where F[^ is evaluated through a dispersive integral as discussed in Eq. (236), and the a- 
meson mass m^. is a free parameter as in the fixed-i unsubtracted RCS dispersion analysis. A 
fit to the t-dependence of RCS data results in nia- ^ 0.6 GeV [75]. The value Ff^{0, 0, 0) is 
then considered as a remaining global fit parameter to be extracted from experiment. It can 
be expressed physically in terms of the magnetic dipole polarizability P : 

4-7r 

F,^^(0,0,0) = —f3. (244) 

The term Fi^{0, 0, 0) in Eq. (243), can be calculated through a dispersion integral and results 
in the value : 

ae^Fi-^(0,0,0) = Z?'^^ = 9.1, (245) 

in units of 10~^ fm^. From the nN contribution of Eq. (245), and the phenomenological 
value P of Eq. (168), one obtains the difference 

p-p^^ = -7.5 , (246) 

which enters in the rhs of Eq. (243). As discussed before, the small total value of the magnetic 
polarizability P comes about by a near cancellation between a large (positive) paramagnetic 
contribution (P'^'^) and a large (negative) diamagnetic contribution (P — P'^^), i.e., the asymp- 
totic part of Fi parametrizes the diamagnetism. 

Turning next to the Q'^ dependence of the asymptotic contribution to Fi, it has been pro- 
posed in [141] to parametrize this part of the non-Bom term F^^{Q'^, v, t) beyond its nN 
dispersive contribution, by an energy independent i-channel pole of the form : 

(Q^ ^, t) - Fr{Q\ V. t) ^ • (247) 

The function f{Q'^) in Eq. (247) can be obtained by evaluating the Ihs of Eq. (247) at the 
point where the GPs are defined, i.e., i/ = and t = —Q^, at finite Q'^. This leads to : 

fm = [A(g^) - Fr(Q')] {l+QVml) , (248) 



87 



where the shorthand Fi(Q^) is defined in Eq. (230). Equations (247) and (248) then lead to 
the following expression for the VCS amplitude Ff^ : 



1 + QVml 



FriQ',u,t)^FriQ^iy,t) + F,m - F^m : (249) 

where the vrA^ contributions F^^(Q'^, u,t) and F[^(Q'^) are calculated through dispersion 
integrals as given by Eqs. (236) and (239) respectively. Consequently, the only unknown 
quantity on the rhs of Eq. (249) is Fi{Q'^), which can be directly used as a fit parameter 
at finite Q^. The quantity Fi{Q'^) can be expressed in terms of the generalized magnetic 
polarizability p(A^i,A/i)o (227) as [140] : 

where (3{Q'^) is the generalized magnetic polarizability, which reduces to the polarizability 
(3 of RCS at = 0. The parametrization of Eq. (249) for Fi then permits to directly ex- 
tract P{Q^) from VCS observables at a fixed Q^. In the following, we consider a convenient 
parametrization of the dependence of P{Q^) in order to provide predictions for VCS ob- 
servables. For this purpose it was proposed in [141] to use a dipole parametrization for the 
difference /3{Q^) - /^^^(Q^), which enters in the rhs of Eq. (249) via Eq. (250), as Qf] : 

- r^(Q') = 7^ {2 , (251) 

(1 + QVA|) 

where the RCS value (/? — P""^) on the rhs is given by Eq. (246). The mass scale A^j in 
Eq. (25 1) determines the dependence, and hence gives us the information how the diamag- 
netism is spatially distributed in the nucleon. Using the dipole parametrization of Eq. (251), 
one can extract A/3 from a fit to VCS data at different values, and check the parametriza- 
tion of Eq. (251) for the asymptotic contribution to I3{Q'^). 

To have some educated guess on the physical value of A/3, we next discuss two microscopic 
calculations of the diamagnetic contribution to the GP P{Q'^). The diamagnetism of the nu- 
cleon is dominated by the pion cloud surrounding the nucleon. This diamagnetic contribution 
has been estimated in Ref. [141] through a DR calculation of the t-channel nn intermediate 
state contribution to Fi. Such a dispersive estimate has been discussed before for RCS in Sec. 
3.7, where it was shown that the asymptotic part of Fi (or equivalently Ai) can be related 
to the 77 ^ TTTT ^ NN process. The dominant contribution is due to the tttt intermediate 
state with spin and isospin zero (/ = J = 0). The generalization to VCS leads then to the 
identification of with the following unsubtracted DR in t at fixed energy u = : 

™.,^2^ 1 7 ,.Im,Fi(O^0.«') 



Frm--Jdt' ■ (252) 



4mJ 



The dipole form displays the l/Q^ behavior at large as expected from perturbative QCD. 
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The imaginary part on the rhs of Eq. (252) has been evaluated in Ref. [141] through the sub- 
processes 7*7 — i> TTTT and tttt — NN. To describe the dependence of the 7*7 tttt 
amplitude, which is dominated by the unitarized Bom amplitude (on the pion), the pion elec- 
tromagnetic form factor was included. The result for this dispersive estimate of F"* through 
^-channel tttt intermediate states is shown in Fig. 29, and compared with the corresponding 
evaluation of Ref. [151] in the linear cr-model (LSM). The LSM calculation overestimates 
the value of Ff''(0) (or equivalently Pas) by about 30% at any realistic value of rria, which is 
a free parameter in this calculation. However, as for the dispersive calculation, it also shows 
a steep dependence. 
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°: -fi 
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Fig. 29. Theoretical estimates of the asymptotic contribution Ff^ : DR calculation [141] of the 
7*7 —s- TTTT ^ NN process (solid curve); linear a-model (LSM) calculation [151] with = 0.5 GeV 
(dotted curve) and m^r = 0.7 GeV (dashed-dotted curve). The dashed curves are dipole parametrizations 
according to Eq. (251), which are fixed to the phenomenological value at = and are shown for two 
values of the mass-scale, = 0.4 GeV (upper dashed curve, nearly coinciding with solid curve) and Ap 
= 0.6 GeV (lower dashed curve). 

Furthermore in Fig. 29, the two model calculations discussed above are compared with the 
dipole parametrization of Eq. (251) for the two values = 0.4 GeV and A3 = 0.6 GeV. It 
is seen that these values are compatible with the microscopic estimates discussed before. In 
particular, the result for A^ = 0.4 GeV is nearly equivalent to the dispersive estimate of tttt 
exchange in the t-channel. The value of the mass scale is small compared to the typical 
scale of Ad ^ 0.84 GeV appearing in the nucleon magnetic (dipole) form factor. This reflects 
the fact that diamagnetism has its physical origin in the pionic degrees of freedom, i.e., the 
diamagnetism is situated in the surface and intermediate region of the nucleon. 

• Dispersive contributions beyond nN to F2 

Though we can write down unsubtracted DRs for all invariant amplitudes (or combinations 
of invariant amplitudes) except for Fi and F5, one might wonder about the quality of our 
approximation to saturate the unsubtracted dispersion integrals by nN intermediate states 
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only. This question is particularly relevant for the amplitude F2, for which we next investigate 
the size of dispersive contributions beyond the vrA^ channel. We start with the case of RCS, 
where one can quantify the higher dispersive corrections to F2, because the value of F^^ 
at the real photon point can be expressed exactly (see Eqs. (229) and (231)) in terms of the 
scalar polarizability (a + /3) as 

Att 1 

FnO,0,0) = --^(a + /3). (253) 

The nN dispersive contribution provides the value 

{a + (3y^ = 11.6 , (254) 

which falls short by about 15 % compared to the Baldin's sum rule value of Eq. (59). The 
remaining part originates from higher dispersive contributions (ttttN, ...) to F2. These higher 
dispersive contributions could be calculated through unitarity, by use of Eq. (237), similarly 
to the ttN contribution. However, the present data for the production of those intermediate 
states (e.g., 7*iV — > ttttN) are still too scarce to evaluate the imaginary parts of the VCS 
amplitude F2 directly. Therefore, we estimate the dispersive contributions beyond nN by an 
energy-independent constant, which is fixed to its phenomenological value at u — t — 0. 
This yields at = : 



An 1 



Fi'''{0,u,t)^F^''{0,u,t) - [(« + /?)- {a + Pr''\, (255) 

which is an exact relation at z/ = t = 0, the point where the polarizabilities are defined. 

The approximation of Eq. (255) to replace the dispersive contributions beyond nN by a 
constant can only be valid if one stays below the thresholds for those higher contributions. 
Since the next threshold beyond nN is nnN, the approximation of Eq. (255) restricts us in 
practice to energies below the A(1232)-resonance. 

We next consider the extension to VCS, and focus our efforts to describe VCS into the 
A(1232)-resonance region. Analogously to Eq. (255) for RCS, the dispersive contributions 
beyond nN are approximated by an energy-independent constant. This constant is fixed at 
arbitrary Q^, u — 0, and t — —Q^, which is the point where the GPs are defined. One thus 
obtains for F^^ [141] : 

F^^{Q\v,i) « Fll^'iQ^v.i) + \F2{Q^) - n''{Q^)\ , (256) 

where F2{Q'^) is defined as in Eq. (230), and can be expressed in terms of GPs through 
relations as given by Eqs. (231) - (234). We saturate the 3 combinations of spin GPs in 
Eqs. (232) - (234) by their nN contribution, and include for the fourth spin GP of Eq. (242) 

also the 7r°-pole contribution. Therefore, we only consider dispersive contributions beyond 
the nN intermediate states for the two scalar GPs, which are therefore the two fit quantities in 
the present DR formalism for VCS. In this way, one can use Eq. (231), to write the difference 
F2(g2) _ F^^{Q^) entering in the rhs of Eq. (256) as : 
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e2 \E + MJ q2 2M 
X { [a(g2) - a-^(g^)] + [/3(g^) - P^'^iQ')] } , (257) 

in terms of the generalized magnetic polarizability (3{Q'^) of Eq. (250), and the generalized 
electric polarizability a{Q'^), which is related to the GP P^^^'^^^^{Q'^) by Eq. (229). 

We stress that Eqs. (249) and (257) are intended to extract the two GPs a{Q'^) and (3{Q'^) 
from VCS observables minimizing the model dependence as much as possible. As discussed 
before for I3{Q'^), we next consider a convenient parametrization of the Q"^ dependence of 
a{Q'^) in order to provide predictions for VCS observables. For this purpose, a dipole form 
has been proposed in Ref. [141] for the difference a{Q'^) — a^^ {Q"^) which enters in the rhs 
of Eq. (257), 

f - 

q;(Q') - a'^^(Q') = ^ ^ , (258) 

where the dependence is governed by the mass scale A^, the second free parameter of the 
DR formaUsm. In Eq. (258), the RCS value 

(a - a^^) = 9.6 , (259) 

is obtained from the phenomenological value of Eq. (168) for a, and from the calculated nN 
contribution, a^^ = 2.5. Using the dipole parametrization of Eq. (258), one can extract the 
free parameter A^ from a fit to VCS data at different Q"^ values. 



4.5 VCS data for the proton and extraction of generalized polarizabilities 



Having set up the dispersion formalism for VCS, we now show the predictions for the dif- 
ferent ep ep7 observables for energies up to the A(1232)-resonance region. The aim of the 
experiments is to extract the 6 GPs of Eqs. (227) and (228) from both unpolarized and polarized 
observables. We will compare the DR results, which take account of the full dependence of the 
ep — > ep7 observables on the energy (q') of the emitted photon, with a low-energy expansion 
(LEX) in q'. In the LEX of observables, only the first three terms of a Taylor expansion in q' are 
taken into account. In such an expansion in q', the experimentally extracted VCS unpolarized 
squared amplitude A^'^^'p takes the form [136] : 

M-P = ilzl + + Mr + O(q') . (260) 

q q 



Due to the low energy theorem (LET), the threshold coefficients 7V/f'!^| and TW^f are known [136], 
and are fully determined from the Bethe-Heitler + Born (BH + B) amplitudes. The information 
on the GPs is contained in , which contains a part originating from the BH+B amplitudes 
and another one which is a linear combination of the GPs, with coefficients determined by the 
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kinematics. The unpolarized observable Ai^"^ can be expressed in terms of 3 structure functions 
PLiiq), PMq), and PiT(q) by [136] : 



Mr - -Mo''^^ = 2i^2|^^i [sPLLiq) - PTT(q)] + [v2 - ^vsj ^2s{l + s)PLT{q) |, (261) 

where K2 is a kinematical factor, e is the virtual photon polarization (in the standard notation 
used in electron scattering), and vi,V2, V3 are kinematical quantities depending on e and q as 
well as on the cm polar and azimuthal angles (62^ and 0, respectively) of the produced real 
photon (for details see Ref. [137]). The 3 unpolarized observables of Eq. (261) can be expressed 
in terms of the 6 GPs as [136,137] : 



Pll = -2V6MGeP^^^'^^^^ , (262) 
P^^ ^ -SGm- f P(^i'^^)i -V2qo P(^i'^2) A ^ (263) 

-xH^Ge + I^Gm P^'''''^\ (264) 

V 2 Q 2 qo 

where Ge and Gm stand for the electric and magnetic nucleon form factors Ge{Q'^) and 
Gm{Q'^), respectively. 

The first VCS experiment was performed at MAMI [131] and the response functions Plt and 
Pll - Ptt/^ were extracted at Q'^ = 0.33 GeV^ by performing a LEX to these VCS data, accord- 
ing to Eq. (261). To test the validity of such a LEX, we show in Fig. 30 the DR predictions for 
the full energy dependence of the non-Bom part of the ep ep-f cross section in the kinemat- 
ics of the MAMI experiment [131]. This energy dependence is compared with the LEX, which 
predicts a linear dependence in q' for the difference between the experimentally measured cross 
section and its BH -1- B contribution. The result of a best fit to the data in the framework of the 
LEX is indicated by the horizontal bands in Fig. 30 for the quantity (d^a - d^(7^"+^°'°) /$q', 
where $ is a phase space factor defined in Ref. [136]. The fivefold differential cross section d^a 
is differential with respect to the electron lab energy and lab angles and the proton cm. angles, 
and stands in all of the following for da / dkf^^ c^^Lfc di^c.m.- It is seen from Fig. 30 that the DR 
results predict only a modest additional energy dependence up to q' ~ 0. 1 GeV and for most of 
the photon angles involved, and therefore support the LEX analysis of [131]. Only for forward 
angles, ^ 0, which is the angular range from which the value of Plt is extracted, the DR 
calculation predicts a stronger energy dependence in the range up to q' ~ 0.1 GeV, as compared 
to the LEX. 

In Fig. 3 1 , we display the response functions Pll — Ptt / £ and Plt at - 0.33 GeV^, which 
have been extracted from the cross section data of Fig. 30 [131], and compare them with the 
corresponding DR calculations. For the electromagnetic form factors in Eqs. (262)-(264) we use 
the Hohler parametrization [ 1 52] as in the analysis of the MAMI experiment [ 1 3 1 ] . In the lower 
panel of Fig. 31, the Q^-dependence of the VCS response function Plt is displayed, which re- 
duces to the magnetic polarizability /3 at the real photon point. At finite Q^,it contains both the 
scalar GP /^(Q^) and the spin GP p(^i'^i)i, as seen from Eq. (264). It is obvious from Fig. 31 
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The data and the shaded bands, representing the best fit to the data within the LEX formalism, are 

from [131]. The solid curves are the DR results taking into account the full q' dependence of the non-Born 
contribution to the cross section. The asymptotic contributions are calculated according to Eqs. (251) and 
(258), with = 0.6 GeV and A^ = 1 GeV, respectively. 

that the structure function P^t results from a large dispersive vrA^ (paramagnetic) contribution, 
which is dominated by A (1232) resonance excitation, and a large asymptotic (diamagnetic) 
contribution to /3 with opposite sign, leading to a relatively small net result. The asymptotic 
contribution is shown in Fig. 31 with the parametrization of Eq. (251) for the values = 0.4 
and Af3 = 0.6 GeV, which were also displayed in Fig. 29. Due to the large cancellation in Plt, its 
dependence is a very sensitive observable to study the interplay of the two mechanisms. In 
particular, one expects a faster fall-off of the asymptotic contribution with in comparison to 
the ttN dispersive contribution, as discussed before. This is highlighted by the measured value 
of Plt at = 0.33 GeV^ [131], which is comparable to the value of Plt at = 0. As seen 
from Fig. 31, this points to an interesting structure in the Q'^ region around 0.05 - 0.1 GeV^, 
where forthcoming data are expected from an experiment at MIT-Bates [133]. 

In the upper panel of Fig. 31, we show the Q^-dependence of the VCS response function 
Pll - Ptt/^, which reduces at the real photon point to the electric polarizability a. At non- 
zero Q'^, Pll is directly proportional to the scalar GP a{Q'^), as seen from Eq. (262), and the 
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response function Ptt of Eq. (263) contains only spin GPs. As is shown by Fig. 31, the tyN 
dispersive contribution to a and to the spin GPs in Ptt are smaller than the asymptotic contri- 
bution to a. At Q'^ - 0, the ttN dispersive and asymptotic contributions to a have the same sign 
and lead to a large value of a, in contrast to P where both contributions have opposite sign and 
largely cancel each other in their sum. 
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Fig. 31. Results for the unpolarized structure functions Pll — Ptt/^ (upper panels), and Plt (lower 
panels), for e = 0.62. Upper left panel : dispersive vrA^ contribution of the GP a (solid curve, SO), disper- 
sive ttN contribution of the spin-flip GPs (dashed curve, SI), and the asymptotic contribution (AS) of a 
according to Eq. (258) with = 1 GeV (dotted curve). Upper right panel : total result for Pll — Ptt/s 
(sum of the three contributions on the upper left panel) for = 1 GeV (solid curve) and Aq, = 1.4 
GeV (dashed curve). Lower left panel : dispersive vrA^ contribution of the GP (3 (solid curve, SO), con- 
tribution of the spin-flip GPs (dashed curve, SI), and the asymptotic contribution (AS) of /3 according to 
Eq. (251) with = 0.6 GeV (dotted curve). Lower right panel : total result for Plt, for = 0.7 GeV 
(dotted curve), A^ = 0.6 GeV (solid curve), and A^ = 0.4 GeV (dashed curve). The RCS data are from 
Ref. [87], and the VCS data at = 0.33 GeV^ from Ref. [131]. 



Increasing the energy, we show in Fig. 32 the DR predictions for photon energies in the 

A(1232)-resonance region. It is seen that the ep ep7 cross section rises strongly when cross- 
ing the pion threshold. In the dispersion relation formalism, which is based on unitarity and 
analyticity, the rise of the cross section with q' below pion threshold, due to virtual nN in- 
termediate states, is connected to the strong rise of the cross section with q' when a real ttN 
intermediate state can be produced. It is furthermore seen from Fig. 32 (lower panel) that the 
region between pion threshold and the A-resonance peak displays an enhanced sensitivity to the 
GPs through the interference with the rising Compton amplitude due to A-resonance excitation. 
For example, at q' ~ 0.2 GeV, the predictions for Plt in the lower right panel of Fig. 31 for 
A^ = 0.4 GeV and A^ = 0.6 GeV give a difference of about 20 % in the non-Born squared am- 
plitude. In contrast, the LEX prescription results in a relative effect for the same two values of 
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Plt of about 10% or less. This is similar to the situation discussed in Sec. 3.4 for RCS, where 
the region between pion threshold and the A-resonance position also provides an enhanced sen- 
sitivity to the polarizabilities and is used to extract those polarizabilities from data using a DR 
formalism. Therefore, the energy region between pion threshold and the A-resonance seems 
promising to measure VCS observables with an increased sensitivity to the GPs. Such an exper- 
iment has been proposed at MAMI and is underway [134]. 




Fig. 32. Upper panel: The differential cross section for the reaction ep epj as function of the out- 
going-photon energy q' in MAMI kinematics : e = 0.62, q = 0.6 GeV, and for 62m = 0°, in plane 
((j) = 0°). The BH + B contribution is given by the dashed-dotted curve. Lower panel: Results for 
{d^a — (^s^BH+Born-j ^ ^^^/ function of q' . The total DR results are obtained with the asymptotic parts 
of Eqs. (251) and (258), using a fixed value of Aq = 1 GeV and for the same three values of as dis- 
played in the lower right plot of Fig. 31, i.e. A^ = 0.7 GeV (dotted curve), A^ = 0.6 GeV (soUd curve), 
and A/3 = 0.4 GeV (dashed curve). In the lower panel, the DR calculations taking into account the full 
energy dependence of the non-Bom contribution (thick curves) are compared to the corresponding results 
within the LEX formalism (thin horizontal curves). The data are from Ref. [131]. 

When crossing the pion threshold, the VCS amplitude acquires an imaginary part due to the 
coupling to the ttN channel. Therefore, single polarization observables become non-zero above 
pion threshold. A particularly relevant observable is the electron single spin asymmetry (SSA), 
which is obtained by flipping the electron beam helicity [137]. For VCS, this observable is 
mainly due to the interference of the real BH -i- B amplitude with the imaginary part of the VCS 
amplitude. As the SSA vanishes in-plane, its measurement requires an out-of-plane experiment. 
Such experiments have been proposed both at MAMI [134] and at MIT-Bates [153]. In Fig. 33, 
the SSA is shown for a kinematics in the A(1232) region, corresponding with W 1.2 GeV. 
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The DR calculation firstly shows that the SSA is quite sizable in the A(1232) region. The SSA, 
which is mainly sensitive to the imaginary part of the VCS amplitude, displays only a rather 
weak dependence on the magnetic GPs (3{Q'^), and shows a modest dependence on a{Q'^). 
Therefore, it provides an excellent cross-check of the dispersive input in the DR formalism for 
VCS, in particular by comparing at the same time the pion and photon electroproduction chan- 
nels through the A region. 



E, = 883 MeV q = 600 MeV q' = 230 MeV (0 = -135° 




■ 10 20 30 40 50 60 70 80 90 
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Fig. 33. Electron single spin asymmetry (SSA) for VCS in MAMI kinematics as function of the photon 
scattering angle. The full dispersion results are shown for the values: Aq = 1 GeV, = 0.6 GeV (solid 
curve), Aa = 1 GeV, A^ = 0.4 GeV (dashed curve), Aq, = 1 GeV, A^ = 0.7 GeV (dotted curve), and 
A« = 1.4 GeV, A^ = 0.4 GeV (dashed-dotted curve). 

Going to higher Q^, the VCS process has also been measured at JLab and data have been 
obtained both below pion threshold at = 1 GeV^ [154], at Q"^ = \.9 GeV^ [155], as well as 
in the resonance region around = I GeV^ [156] (see Ref. [132] for a short review of these 
JLab data). 

In Fig. 34, we show the results for the ep epj reaction in the resonance region at = 
1 GeV^ and at a backward angle. These are the first VCS measurements ever performed in 
the resonance region. We also display the DR calculations of [141] for the cross section. The 
data clearly show the excitation of the A (1232) resonance, and display a second and third res- 
onance region, mainly due to the excitations of the Dis{1520) and Fi5(1680) resonances. The 
DR calculations reproduce well the A(1232) region. Due to scarce information for the disper- 
sive input above the A (1232) resonance, the DR calculations cannot be extended at present into 
the second and third resonance regions. Between pion threshold and the A(1232) resonance, 
the calculations show a sizable sensitivity to the GPs, in particular to Pll in this backward an- 
gle kinematics, and seem very promising to extract information on the electric polarizability. 
The precise extraction of GPs from VCS data at these higher values of Q'^, requires an accu- 
rate knowledge of the nucleon electromagnetic form factors (FFs) in this region. For the proton 
electromagnetic form factors, we use the new empirical fit of [157], which includes the recent 
high accuracy measurements performed at JLab for the ratio of proton electric FF to the 
magnetic FF Gm in the Q"^ range 0.4 - 5.6 GeV^ [70,71]. From Fig. 34, one sees that a good 
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description of the JLab data is obtained by the values = 1.0 GeV and hp = 0.45 GeV. 




Fig. 34. The differential cross sections for the ep — > ep-y reaction as function of the cm. energy W in 
JLab kinematics : = 4.032 GeV, = l.o GeV^, and for fixed scattering angle 92^ = -167.2°, 
for different out-of-plane angles 0. The BH + B contribution is given by the dashed curve. The total 
DR result is shown by the solid curve (limited to W < 1.25 GeV) for the values = 1.0 GeV and 
Al3 = 0.45 GeV. The data are from Ref [132]. 

Besides the measurement in the resonance region, the ep — > ep7 reaction has also been mea- 
sured at JLab below pion threshold for three values of the outgoing photon energy at = 1 
GeV^ [154], and at = 1.9 GeV^ [155]. For those kinematics, we show in Fig. 35 the dif- 
ferential cross sections as well as the non-Bom effect relative to the BH -i- B cross section. It 
is seen from Fig. 35 that the sensitivity to the GPs is largest where the BH -i- B cross section 
becomes small, in particular in the angular region between 0° and 50°. In Fig. 35, we show 
the non-Bom effect for different values of the GPs. From the JLab data below pion threshold, 
the two unpolarized structure functions Pn — Ptt/^ and Plt have been extracted at = 1 
GeV^ and at Q'^ = 1.9 GeV^ [132]. For this extraction below pion threshold, both the LEX and 
the DR formalisms can be used. A nice agreement between the results of both methods for the 
structure functions was found in Ref. [132]. The preliminary results at = 1 GeV^ and at Q'^ = 
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Fig. 35. The differential cross section for the reaction ep — > epj as function of the photon scattering 
angle, at different values of the outgoing-photon energy in JLab kinematics: Q"^ = 1 GeV^ and e = 0.95 
(left panels) and = 1.9 GeV^ and e = 0.88 (right panels). The BH + B cross sections are shown 
by the dashed-dotted curves. The DR results are displayed with the asymptotic terms parametrized as in 
Eqs. (258) and (251), using the values : = 1 GeV and = 0.6 GeV (solid curves), A^ = 1 GeV and 
A^ = 0.4 GeV (dashed curves), A^ = 1.4 GeV and A^ = 0.6 GeV (dotted curves). 

1.9 GeV^ for Pll Q ^lt are displayed in Fig. 36, alongside the RCS point and the results 
at = 0.33 GeV^. By dividing out the form factor Ge, one sees from Eq. (262) that Pll is 
proportional to the electric GP a{Q'^), whereas Plt is proportional to the magnetic GP 13 {Q"^) 
plus some correction due to the spin flip GP p(^i'^i)i vvhich turns out to be small in the DR for- 
malism as discussed further on. One sees from Fig. 36 that the best fit value for ~ 0.92 GeV 
yields an electric polarizability which is dominated by the asymptotic contribution and has a 
similar behavior as the dipole form factor. However, the best fit value for ^ 0.66 GeV 
is substantially lower, indicating that the diamagnetism, which is related to pionic degrees of 
freedom, drops faster with Q"^. One nicely sees that the data confirm the interplay between para- 
and dia-magnetism in [3 as function of Q^. 
Until now, we discussed only unpolarized YCS observables. An unpolarized YCS experi- 



The present experiments, which are performed at a fixed value of e only measure the combination 
Pll — Ptt/^- To extract Pll from these data, we calculate the relatively small (spin-flip) contribution 
Ptt (shown by the curves labeled SI on the left panel of Fig. 31) in the DR formalism and subtract it 
from the measured value. 
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Fig. 36. Results for the unpolarized VCS structure functions Pll (left panel) and Plt (right panel) di- 
vided by the proton electric form factor. Dashed hnes: dispersive ttN contributions. Dotted lines: asymp- 
totic contributions calculated according to Eqs. (251) and (258) with Aq, = 0.92 GeV (left panel) and 
= 0.66 GeV (right panel). Solid curves: total results, sum of the dispersive and asymptotic contribu- 
tions. The RCS data are from Ref. [87], the VCS MAMI data at = 0.33 GeV^ are from Ref. [131], 
and the preliminary VCS JLab data at = i GeV^ and = 1.9 GeV^ from Ref. [132] (inner error 
bars are statistical errors only, outer error bars include systematical errors). The values for P^i^ at > 
were extracted by use of a dispersive estimate for the not yet separated Ptt contribution. 

ment gives access to only 3 combinations of the 6 GPs, as given by Eqs. (262) - (264). It was 
shown in Ref. [158] that VCS double polarization observables with polarized lepton beam and 
polarized target (or recoil) nucleon, will allow us to measure three more combinations of GPs. 
Therefore a measurement of unpolarized VCS observables (at different values of e) and of 3 
double-polarization observables will give the possibility to disentangle all 6 GPs. The VCS 
double polarization observables, which are denoted by AA4{h, i) for an electron of helicity h, 
are defined as the difference of the squared amplitudes for recoil (or target) proton spin orien- 
tation in the direction and opposite to the axis i (i = x, y, z), where the z-direction is chosen 
along the virtual photon momentum (see Ref. [158] for details). In a LEX, this polarized squared 
amplitude yields : 



AM 



exp 



V2 



+ AMrP + 0(q'). 



(265) 



Analogous to the unpolarized squared amplitude (260), the threshold coefficients ATW^I, AAit 
are known due to the LET. It was found in Ref. [158] that the polarized squared amplitude 
AM^"^ can be expressed in terms of three new structure functions Plrp{q), P^^(q), and P^-^(q). 
These new structure functions are related to the spin GPs according to Refs. [158,137] : 
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= ^C^pC^l,^!)! _ ^GgpCiV^l-Ml)!^ (266) 

2% Q 

P'lt- ~QGm + Ge P(^^'^^)\ (267) 

^ (^^^''^'^' " yf ?oP(^^'^^)^] . (268) 



While P£y and can be accessed by in-plane kinematics ((f) — 0°), the measurement of P'fj^ 
requires an out-of-plane experiment. 

In Fig. 37, we show the dispersion results for the double polarization observables, with po- 
larized electron and by measuring the recoil proton polarization either along the virtual photon 
direction (^-direction) or parallel to the reaction plane and perpendicular to the virtual pho- 
ton (x-direction). The double polarization asymmetries are quite large (due to a non- vanishing 
asymmetry for the BH + B mechanism), but the DR calculations show only small relative ef- 
fects due to the spin GPs below pion threshold. However, a heavy-baryon chiral perturbation 
theory (HBChPT) calculation to 0{p^) [159] shows a significantly larger effect due to larger 
values of the spin GPs in this calculation, as will be discussed in the next section. Although 
these double polarization observables are tough to measure, a first test experiment is already 
planned at MAMI [134]. 



q=600MeV q'=111.5 MeV e=0.62 (|)=0° 




Fig. 37. VCS double-polarization asymmetry (polarized electron, recoil proton polarization along either 
the z- or x- directions in the cm. frame) in MAMI kinematics as function of the photon scattering angle. 
The dotted curves correspond to the BH+B contribution. The solid curves show the total DR results for 
the values of Aq, = 1 GeV, = 0.6 GeV. The dashed curves are the HBChPT predictions from [159]. 
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4.6 Physics content of the nucleon generalized polarizabilities 



Having discussed the present status of the VCS experiments and the combinations of GPs 
which have been extracted from such experiments so far, we now turn in some more detail to 
the physics content of these GPs and compare different model predictions. 

We start our discussion with the VCS unpolarized structure functions as shown in Fig. 38, 
and compare the DR results of Ref. [141] with the 0{p^) HBChPT calculations [159,160]. The 
DR results have been shown before in Fig. 31, where we discussed the different mass scales 
parametrizing the asymptotic parts in the GPs a and (3. In Fig. 38, we show in addition the 
effect of the spin GPs on these response functions and compare them with the corresponding 
calculation in HBChPT. One notices that the effect of the spin GPs is much smaller in the DR 
calculation than in 0{p^) HBChPT, in particular for the spin GPs entering Ptt- The good agree- 
ment with the data found in the 0{p^) HBChPT calculation is for an important part due to the 
larger size of the spin GPs in this calculation. 



Pll-Ptt/e (GeV-^) 



Pll-Ptt/e (GeV-^) 




(GeV^) 



(GeV^) 



Fig. 38. Comparison between the VCS unpolarized structure functions calculated within the DR for- 
malism [141] and the 0{p^) HBChPT [159,160]. Upper panels : results for the Pll - Ptt/£ structure 
function for e = 0.62 predicted from DR (left) and 0{p^) HBChPT (right). For the DR predictions, the 
contribution from the electric GP a{Q'^) for A„ = 1.4 GeV (dashed curve), is compared with the sum of 
the contributions from the scalar and spin-flip GPs (solid curve). Lower panels : results for Plt within 
DR (left) and 0{p^) HBChPT (right). In the left panel, the contribution from /^(Q^) is shown for the 
values = 0.6 GeV (dashed curve) and = 0.4 GeV (dotted curve). The total results, sum of the 
contributions from scalar and spin-flip GPs, are shown for A^ = 0.6 GeV (solid curve) and for A/-j = 0.4 
GeV (dashed-dotted curves). In the right panel, the predictions from HBChPT are shown for the con- 
tribution from f3{Q^) alone (dashed curve), and for the total result (soUd curve), which includes the GP 
p(Li,Li)i Tjjg f^Qjj^ [37]^ ^jjg Y(3s q2 ^ Q 33 QgY2 ^^^^ j^^f ^3^^ 
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The comparison between the spin independent GPs in both calculations is shown in Fig. 39. 
From this figure, we see a qualitative agreement between both the DR and the 0{p^) HBChPT 
results. In particular, we see that in both calculations the Q"^ dependence of the electric and mag- 
netic GPs is quite different. The electric GP shows a rather smooth Q'^ behavior, much as the 
nucleon electric form factor, whereas the magnetic GP has a characteristic structure at small Q^. 
In the DR calculation, this results due to a cancellation between a large paramagnetic A contri- 
bution and a diamagnetic contribution (due to t-channel tttt exchange) which have a different 
behavior, as was already noticed in the early effective Lagrangian calculation of Ref. [161]. In 
the 0{p^) HBChPT, this structure in /^(Q^), at low Q"^, results from nN loop effects. By Fourier 
transforming the GPs ol{Q'^) [^and PiQ"^) in the Breit frame, it was argued in [162] that one 
obtains a spatial distribution of the induced electric polarization a(r) and magnetization /3(r) 
of the nucleon. The picture which then emerges from the nN loop contribution in the HBChPT 
calculation is as expected from a classical interpretation of diamagnetism. Due to a change in the 
external magnetic field, pionic currents start circulating around the nucleon, and give rise to an 
induced magnetization, opposite to the applied field. This diamagnetic effect leads at distances 
r > l/rriTr to an negative value for /3(r), whereas for distances r < l/rriT, the paramagnetism 
dominates and /3(r) is positive. Therefore, as the momentum transfer increases, the nega- 
tive long-distance contribution to the magnetic GP due to the pion cloud, no longer contributes 
and hence P{Q'^) increases. This nicely explains the positive slope of P{Q'^) at = and the 
characteristic turn-over at low in the HBChPT calculation as is seen in Fig. 39. Hence it 
will be interesting to see the results of a measurement around ~ 0.05 — 0.1 GeV^ for Plt 
performed at MIT-Bates [133] to reveal the nature of the diamagnetism in the nucleon. 



a(Q') ^(Q^) (10-^fm') 




0.2 0.4 0.2 0.4 

(GeV^) (GeV^) 

Fig. 39. Left panel : comparison between the results for the electric GP a{Q'^) predicted by the DR 
formalism for = 1.4 GeV (solid curve) and 0{p^) HBChPT (dashed curve). Right panel : comparison 
between the results for the magnetic GP /3{Q'^) predicted by the DR formalism for A/3 = 0.6 GeV (full 
line) and = 0.4 (dashed-dotted line), and by 0{p^) HBChPT (dashed line). 

We next discuss the spin-flip GPs. In Fig. 40, we show the dispersive and 7r°-pole contribu- 
tions to the 4 spin GPs as well as their sum, according to the calculations of [141,143]. For the 
presentation, we multiply in Fig. 40 the GPs p(^i'A^2)i p(Mi,L2)i ^^^j^ -^^ ^^^^^ better 

compare the dependence when including the 7r°-pole contribution, which itself drops very 

In the notation of Ref. [162], a{Q^) is denoted as the so-called longitudinal electric GP a^, to distin- 
guish it from a higher order (in the outgoing photon energy) so-called transverse electric GP ax- 
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fast with Q^. The vr^-pole does not contribute to the GP p(^i'^i)i^ but is seen to dominate the 
other three spin GPs. It is however possible to find, besides the GP p(^i:^i)i^ the two combi- 
nations given by Eqs. (233) and (234) of the remaining three spin GPs, for which the 7r°-pole 
contribution drops out [143]. 






0.25 0.5 0.75 1 

(GeV^) 

Q P^^'-"«' (10-'fm^) 



0.25 0.5 0.75 1 

(GeV^) 




0.25 0.5 0.75 1 



(GeV^) 



0.25 0.5 0.75 1 



(GeV^) 



Fig. 40. Q^-dependence of the spin-flip GPs as calculated in Refs. [141,143]. The dashed curves corre- 
spond to the dispersive vrA^ contribution, the dotted curves show the vr'^-pole contribution, and the solid 
curves are the sum of the dispersive and 7r°-pole contributions. For comparison, we also show the 7r°-pole 
contribution when setting the 7r°7*7 form factor equal to 1 (dashed-dotted curves). Note that p(^i'^i)i 
has no Tr'^-pole contribution. 



In Fig. 41 we show the results of the dispersive contribution to the four spin GPs, and com- 
pare them to the results of the nonrelativistic constituent quark model [163], the HBChPT to 
[159,160], the recent HBChPT calculation to C(/) [164], and the linear a-model [151]. 
The constituent quark model (CQM) calculation gives negligibly small contributions for the 
GPs P(^i'^i)i and pi^^^L2)i^ whereas the GPs P(a^i>a^i)i and p(li,m2)i j-eceive their dominant 
contribution from the excitation of the A(1232) (Ml -> Ml transition) and the A^* (1520) 
{El M2 or LI M2 transitions), respectively. The smallness of p(^i'^i)i and p(^i.i2)i 
in the CQM can be understood by noting that those two GPs can be expressed in terms of a GP 
which involves a transition from LO (Coulomb monopole) to Ml, through a crossing symmetry 
relation [140] as : 

2 2 
^q_p(Ll,Ll)l ^ p(Ml,LO)l _^ _q_p(Ml,L2)l ^-259) 

The GPs on the rhs of Eq. (269) encode the response to a static magnetic dipole field (Ml) of the 
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nucleon charge density (LO) or the electric quadrupole density (L2). In a non-relativistic CQM 
calculation [136,163], the only response to an applied static magnetic field is the alignment of 
the quark spins, whereas the charge density or electric quadrupole density remain unchanged. 
Therefore, both GPs p{^i'^o)i and p(^i'^2)i vanishingly small in the quark model, as well 
as p(^i'^i)i through Eq. (269). Consequently, p(^^'^^)^ and p(A^i.^2)i promising observ- 
ables to study the effects of the pion cloud surrounding the nucleon. A large contribution of 
pionic effects for these GPs is indeed observed in the HBChPT and in the linear ci-model calcu- 
lations. One furthermore notices from Fig. 41 that the 0{p^) HBChPT predicts a rather strong 
increase with for the GPs p(^i'^i)i and p(A^i,Mi)i p^j. p(Li,Li)i ^^^^ j-gg^n confirmed by 
the calculation [164]. For the GP pm,Mi)i^ ^ ^^s found in Ref. [164] that the 0{p^) 

calculation gives a large reduction compared to the 0{p'^) result, and calls the convergence of 
the HBChPT result for this observable into question. The linear cr-model, which takes account 
of part of the higher order terms of a consistent chiral expansion, in general results in smaller 
values for the GPs p(^i'^^)i and p(^i.^i)i compared with the corresponding calculations to 
leading order in HBChPT. For the GP p(^i.A^2)i^ ^^^^^ q2 ^ which is related to the spin 
polarizability 73 through Eq. (229), was reported in Sec. 3.11. In particular, we can notice that 
the 0{p^) HBChPT result yields a relatively large correction, bringing it in better agreement 
with the DR result. From Fig. 41 one notices that the dependence of the 0{p^) HBChPT 
calculation for the GP p(^i'^2)i rather weak [164], and results in a near constant reduction 
for this observable compared to the 0{p^) calculation. 

The comparison in Fig. 41 clearly indicates that a satisfying theoretical description of the 
spin-dependent GPs over a larger range in is still a challenging task. This calls for VCS 
experiments which are sensitive to the spin-dependent GPs. Two types of experiments can be 
envisaged in this regard. Firstly, one notices from Eq. (261) that an unpolarized VCS experi- 
ment at different values of e (by varying the beam energy) allows one to disentangle the re- 
sponse functions Pll and Ptt- The latter contains the combination of the spin GPs p(^i'^i)i 
and p{^i'Af2)i gjyej^ \jy £q (263). In Fig. 42, we show the response function Ptt and compare 
the DR predictions [143,141] with the 0{p^) HBChPT result [160] and the 0{p^) HBChPT 
result [164]. One notices large corrections at O(p^) to the HBChPT result. Therefore, the main 
difference between the DR result and the 0{p^) HBChPT result for the measured response 
function P^ — Ptt/^, as shown in the upper panels of Fig. 38, is largely reduced by the 0{p^) 
HBChPT calculations. It will be very worthwhile to directly measure the response function 
Ptt which will provide an interesting check on our understanding of the spin densities of the 
nucleon, and allow to extract the electric polarizability a{Q'^) unambiguously from the mea- 
surement of Pll — Ptt/s- 

To access the other spin GPs, which do not appear in Ptt it was discussed before that one has 
to resort to double polarization observables. It was shown in Fig. 37 that such observables are 
particularly sensitive to the different predictions for spin GPs, and are very promising to mea- 
sure in the near future [134], so as to sharpen our understanding of the spin dependent response 
of the nucleon to an applied electromagnetic field. 
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Fig. 41. Results for the spin-flip GPs excluding the 7r"-pole contribution in different model calculations. 
The solid curves correspond to the dispersive ttN contribution [141,143]. The thin dashed curves show 
the results of 0{p^) HBChPT [159], whereas the thick dashed curves for p(LhLi)i^ p(Mi,Mi)i^ ^j^^ 
p(Li,M2)i jjjg 0{p'^) HBChPT results [164]. The dashed-dotted curves correspond to the predictions 
of the linear cr-model [151], and the dotted curves are the results of the nonrelativistic constituent quark 
model [163]. Note that the constituent quark model results (CQM) for p{^i.^i)i and p(A^i'^2)i 
multiplied (for visibility) by a factor 100. 
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Fig. 42. Results for the VCS structure function Ptt- Dotted curve : 0{p^) HBChPT [160]; dashed curve 
: 0{p^) HBChPT [164]; soUd curve : dispersive evaluation [143,141]. 
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5 Conclusions and perspectives 



In this review, we have applied dispersion relations to real and virtual Compton scattering 
processes off a nucleon as a powerful tool to connect different observables and to extract nu- 

cleon structure quantities. 

For forward real Compton scattering, sum rules directly connect low energy quantities to the 
polarized or unpolarized total absorption cross sections. We discussed in some detail the recent 
evaluations of the Baldin sum rule and the status of the GDH sum rule. The latter involves an 
integral over the helicity difference cross section (Ji/2 — 0-3/2, for photon and proton helicities 
having the same or opposite signs. This helicity difference cross section for the proton has now 
been measured at MAMl and ELSA through the resonance region, up to W < 2.5 GeV. It dis- 
plays a region around pion threshold which is dominated by S-wave pion production, for which 
(7i/2 dominates. Furthermore, these data clearly exhibit three resonance regions with dominance 
of as/2. By performing the GDH integral up to 1^ < 2 GeV, one overestimates the sum rule 
value for the proton by about 15 %, indicating that the anomalous magnetic moment is mostly 
related to the low-lying degrees of freedom. A measurement of ai/2 — 0^12 at energies up to 
< 9 GeV will be performed in the near future at SLAC, in order to find out whether the 
present "oversaturation" of the sum rule will be removed by high-energy contributions. Such 
an experiment will be quite important, because it will test both our understanding of soft Regge 
physics in the spin-dependent forward Compton amplitude and the validity of high-energy ex- 
trapolations of DIS data at large to the real photon point. 

For the neutron, the convergence of the GDH sum rule is less clear at the moment because of 
a lack of data. Theoretical estimates based on our present knowledge of pion photoproduction 
multipoles yield only about 85 % of the sum rule value. This may point to systematic defi- 
ciencies in these multipoles, which have mostly been obtained from experiments on a deuteron 
target, or to large contributions from higher intermediate states, such as two pions. It is there- 
fore of extreme interest to see the outcome of dedicated experiments on the neutron which are 
planned in the near future at MAMI and GRAAL. 

Extending the sum rules to forward scattering of spacelike virtual photons, we have shown 
how to relate nucleon structure quantities to the inclusive electroproduction cross sections. The 
unpolarized cross section (weighted with leads to a generalization of Baldin's sum rule, 
whereas the polarized cross sections (weighted with 1 / v^) lead to two nucleon spin polarizabil- 
ities. We estimated these quantities at low and intermediate by a phenomenological model 
(MAID), and at large by the corresponding moments of DIS structure functions. As a result, 
we find that a transition occurs around ^ 1 — 2 GeV^ from a resonance-dominated descrip- 
tion at lower to a partonic description at larger . Furthermore, we also studied the general- 
ized GDH integrals, using very recent experimental results at intermediate values measured 
at SLAC, HERMES, JLab/CLAS, and JLab/HallA. In particular, the JLab/CLAS data for the 
helicity difference cross section (T1/2 — 0-3/2 of the proton in the range Q'^ ~ 0.15 — 1.2 GeV^, 
clearly demonstrate a sign change from a large negative value at low Q^, where 0-3/2 dominates 
due to resonance excitation, to the positive DIS value at larger Q^, where (T1/2 survives. We have 
shown that this transition can also be nicely understood in a quantitative way. For the proton- 
neutron difference, where isospin 3/2 resonances such as the A drop out, the validity of chiral 
perturbation theory (ChPT) extends towards somewhat larger Q^, and there is hope to bridge 
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the gap between ChPT and perturbative QCD, which eventually leads to the well established 
Bjorken sum rule at large . 

In Sec. 3, we extended the dispersion formalism for forward scattering to real Compton scat- 
tering (RCS) on the nucleon for all angles. At low photon energies, this process has a well- 
known low-energy limit, the Thomson term, which is determined by the total mass and electric 
charge of the system. Moving to larger photon energies, one can identify the higher order terms 
in a low energy expansion (LEX) with the response of the nucleon to an external electromag- 
netic field, parametrized by dipole and higher order nucleon polarizabilities. However, such a 
LEX is only valid up to about 80 MeV photon energy, and a direct experiment in this energy 
range would have to be extremely precise to disentangle the small effects due to nucleon polar- 
izabilities. In practice, one also has to include experiments at higher energies, up to and above 
the pion threshold, and to rely on dispersion relations to analyze the experiments. We have re- 
viewed and compared several such dispersion relation formalisms for RCS. 

In the literature, most of the recent experiments have been analyzed using unsubtracted fixed- 
t dispersion relations for the six RCS amplitudes. In such an approach, one has to estimate the 
asymptotic contributions for two of the six RCS amplitudes for which the unsubtracted disper- 
sion integrals do not converge. These asymptotic contributions can be parametrized as energy 
independent t-channel poles. In such parametrizations, the most important contributions are 
the 7r°-pole for one of the spin-dependent amplitudes and a 'V'-pole for a spin-independent 
amplitude. This procedure is relatively safe for the 7r°-pole which is well established both ex- 
perimentally and theoretically. However, since the a-meson mass and coupling constants enter 
as free parameters in such a formalism, the "cr"-pole introduces a considerable model depen- 
dence. Instead we replace the a-meson by existing physical information on the / = J = 
part of the two-pion spectrum, within the formalism of fixed-t dispersion relations. This has 
been achieved by subtracting the fixed-t dispersion relations (at = 0) and by evaluating the 
subtraction functions through a dispersion relation in the variable t. The absorptive parts enter- 
ing the t-channel dispersion integrals can be saturated by tttt intermediate states in the reaction 
77 TTTT — iViV, constructed by means of the phenomenological information on the 77 nn 
and TTTT — > NN subprocesses. In this way we found that a consistent description for Compton 
scattering data at low energy can be achieved in both formalisms. Going to higher energies 
and backward scattering angles, a large part of the integration range lies outside the physical 
region, and the full amplitude has to be constructed by an analytical continuation of the par- 
tial wave expansion. Since this expansion converges only in a limited range, the application 
of fixed-t dispersion relations is restricted in practice to energies up to the A-resonance and 
to forward angles. In order to overcome this shortcoming, we also studied fixed-angle disper- 
sion relations, in which case the integration range of the s-channel contribution falls into the 
physical region. The t-channel dispersion integrals can be reconstructed from a partial wave ex- 
pansion which converges for angles > 100". Furthermore, such fixed-angle dispersion relations 
can quantitatively explain the large value for the difference of the electric and magnetic dipole 
polarizabilities, a — (3, without invoking a cr -meson contribution. Evaluated at ^ = 180°, the 
predictions are a - /3 = (10.7 ± 0.2) • 10-^fm^ and 7^ = (-38.8 ± 1.8) • lO^^fm^ for the 
backward spin polarizability. 

In conclusion, fixed-t and fixed angle subtracted dispersion relations nicely complement each 
other, the former being convergent at small scattering angles and the latter at large scattering 
angles. We applied this combined formalism to all existing data. Below pion threshold, we 
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found that all methods nicely agree. This comparison proves that the polarizabilities can indeed 
be extracted with a minimum of model dependence for the energy range below the A reso- 
nance. However, subtracted dispersion relations also provide a quantitative description of the 
data through the A resonance. 

We have furthermore shown that the sensitivity to the backward spin polarizability can 
be substantially increased by an experiment with polarized photons hitting a polarized proton 
target. Such an experiment, although challenging to perform, could become feasible in the near 
future, and can teach us more about the spin response of the nucleon to a static electromagnetic 
field. 

In Sec. 4, we have extended the dispersion relation formalism to virtual Compton scatter- 
ing (VCS) off a proton target, as a tool to extract generalized polarizabilities (GPs) from VCS 
observables over a large energy range. The way we evaluated our dispersive integrals using 
'kN intermediate states, allows us to apply the present formalism for VCS observables through 
the A (1232) -resonance region. The presented dispersion relation framework, when applied at a 
fixed value of Q^, involves two free parameters, which can be expressed in terms of the electric 
and magnetic GPs, and which are to be extracted from a fit to VCS data. 

We confronted our dispersive calculations with existing VCS data taken at MAMI and JLab 
below pion threshold. Our dispersion relation formalism yields results consistent with the low- 
energy expansion analysis for photon energies up to about 100 MeV. When increasing the pho- 
ton energy, the dispersive calculations show that the region between pion threshold and the 
A-resonance peak displays an enhanced sensitivity to the GPs. We also compared our disper- 
sion relation calculations to JLab data taken at higher photon energies, through the A(1232)- 
resonance region, and found a good agreement. The extraction of GPs from the preliminary 
JLab data below and above pion threshold yields consistent results. These data indicate a Q"^ 
dependence of the electric GP similar to a dipole form factor, whereas the magnetic GP follows 
a more complicated behavior. As was already shown for RCS, the magnetic dipole transi- 
tion involves a strong cancellation between a diamagnetic mechanism due to pion cloud effects 
and a paramagnetic contribution due to nucleon resonance excitation. Since the cloud effects 
have a considerably longer range in space than the resonance structures, the Q"^ behavior of the 
magnetic GP is able to disentangle both physical mechanisms, which is already displayed in 
the existing data. Given this initial success, future experiments to measure VCS observables in 
the A-energy region hold the promise to extract GPs with an enhanced precision, within the 
dispersion relation formalism presented here. 

Besides the VCS experiments without polarization degrees of freedom, which give access to 
a combination of only 3 of 6 GPs, we investigated the potential of double polarization VCS ob- 
servables. In fact, a first double polarization experiment is now underway at MAMI. Although 
such investigations will be challenging indeed, they are prerequisite to access and quantify the 
full set of scalar and spin GPs of the nucleon. 

In conclusion we find that dispersion relations are indeed a powerful tool to analyze real and 
virtual Compton scattering processes, linking low-energy structure quantities to the excitation 
spectrum of the nucleon. Though the experiments with virtual photons have only become fea- 
sible very recently, they have opened up a new and systematic way to map out, in quantitative 
detail, the transition from hadronic degrees of freedom at low virtuality to partonic degrees of 
freedom at large virtuality. We are looking forward to increasing theoretical and experimental 
activities in the fields of both real and virtual Compton scattering, and hope that the present 
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review will be useful to stimulate and analyze such further work. 
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A ^-channel tttt exchange 



We express the invariant amplitudes Ai{i',t) (i = 1,.",6) in terms of the t-channel helicity 
amplitudes Tl^^ x (i^, for which we have found the expressions 
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We decompose the i-channel helicity amplitudes for 77 — > A/^A?^ into a partial wave series. 



(A.2) 



where dj^^j^^ are Wigner ci-functions and is the scattering angle in the t-channel, which is 



related to the invariants v and t by cos 6*4 = 4M u/^/iy/t — AM'^. We calculate the imagi- 
nary parts of the i-channel helicity amplitudes Tl^)^_ )^, ^^{u, t) through the unitarity equation 

by inserting tttt intermediate states, which should give the dominant contribution below KK 
threshold. 
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Combining the partial wave expansion for 77 — > tttt 
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■ Pp{cos9^^), (A.4) 
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and the partial wave expansion for tttt — > NN, 



2J+1 



-, J (tttt ^ AfAT) 
Aw 



We can now construct the imaginary parts of the Compton t-channel partial waves, 



2 IinT^^^^^ ^/^^^ (^) 



1 Ptt 



(Stt) 



T, 



/ (77— » tttt) 



A. 



T, 



J (tttt -» AT AT) 



A^ 



it) 



(A.6) 



The partial wave amplitudes Qf Eq_ (A.5) are related to the amplitudes f±{t) of 

Frazer and Fulco [165] by the relations 



J (tttt ^ NN) 



Ajv=0 



it) 



16% 



Pn 



Pn 



(A.7) 



with Pat and p^r the cm. momenta of nucleon and pion respectively (p^ = Jt/A — AP and 



= ^t/A — ml). For the reaction 77 tttt, we will use the partial wave amplitudes Fj a^, (t), 
which are related to those of Eq. (A.4) by 



T, 



J (77 tttt) 



it) 



V2J + 1 



(A.8) 



Inserting the partial-wave expansion of Eq. (A. 2) into Eq. (A.l), we can finally express the 2n 
t-channel contributions lm.tAi{v, t)'^^ by the partial wave amplitudes for the reactions 77 — > tttt 
and TTTT — > NN, 



2^ 1 



Fj=oAMt)fi='*it) 



tViPN 

-SU + M^)F,=2A,=o(t) fi='*{t) - JyM^.^p^F,=2A,=o(t) f^='*{t) 



'15^ 
2t''Vt 

2 tViPN 



AMu'Fj=,j,-o{t)fl='%t), 



Fj=2A,=2{t) fi='*{t) - MFj=2A,=2{t) fi=''*{t) 



ImtA5(i/,i)'^ = - 



'15 M 
~2u/t 



piFj=2AMt) f-="*it) , 
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We note that the s-wave {J — 0) component of the 2n intermediate states contributes only to 
Ai and that only waves with J > 4 contribute to the amplitude A4. 



B Tensor basis 

In writing down a gauge-invariant tensor basis for VCS, we use the combinations of the four- 
momenta given in Eq. (157), 



The 12 independent tensors pf' entering the VCS amplitude of Eq. (222) and introduced in 
Ref. [140], are given by : 



p = - (p + p') , 



K=l{q + q') . 



(B.l) 



pf = -2MuQ^g'"' - 2Mu q^q" + 2Q^P''q'>' + 2q-q' P^^q", 

pf = P''q^$- ^ (P'^7'^ - P'^7^) - Mv q^Y " ^ 75 e^'^'^^ K^^p, 



p^" = -Sq.q'P^'P'' + 4Mi/ [P^q" + P^q'^) + 4Mq-q' (P^'Y + P^l") 

-4MV (g'/^7'' + ^'^7'^) + i4Mt/ (q^a'^'Ko, - q'a'^'K^ + q-q'a'"') 

+i4M5.?'75£'^'^"^i^a7/3, 



= (p/^g- - P'^q'^^^- q.q' (^P^'Y - P"Y) + Mu (q'^Y - (fl"), 




pf = 2Mv (P^'q" - P^q'^) - 2Mq-q' (P^Y - P^Y) + '^M^^ {q'^Y - Q^Y) 
+i 2q-q' (P'^cr'^^X^ + P'^a^'^Ka) - i 2Mv (q'^^a'^'^Ka + q'^a^'^Ka), 
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-2i (g''^(7^"Xa - gV^^Xa + q-q'f^""), 

= 2g2p'^P^ + 2Mi/ P^g^ - 2MQ^P^'Y - SM^i/ + i 2Mi/ q''a'"^K^ 

+i Q2 {P^'a^^'K^ + P'^a^'^'K^ - Mu a^"") - i MQ^ 75 e'"''"^Ka^/3, (B.2) 

where we follow the conventions of Bjorken and Drell [145], in particular a^'' — i/2 [7^*, Y] 
and £0123 = +1- 
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